Technical  Report  CERC-95-4 
May  1995 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Directional  Wind  Wave  Characteristics 
at  Harvest  Platform 

by  Charles  E.  Long 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


19950628  025 

DTIf  QUALTIT  mSPECTED  8 


Prepared  for  Headquarters,  U.S.  Army  Corps  of  Engineers 


The  contents  of  this  report  are  not  to  be  used  for  advertising, 
publication,  or  promotional  purposes.  Citation  of  trade  names 
does  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 


PRINTED  ON  RECYCLED  PAPER 


Technical  Report  CERC-95-4 
May  1995 


Directional  Wind  Wave  Characteristics 
at  Harvest  Platform 

by  Charles  E.  Long 

U.S.  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Accesion  For 

NTIS  CRA&I 
OTIC  TAB 
Unar'.nounced 
Justification 


By . 

Distribution  / 


□ 

□ 


Availability  Codes 


DIst 

Avail  and/or 

Special 

Final  report 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


Under 


Work  Unit  32484 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


HEADQUARTERS 

BULDilG 


FOR  INFORMATION  CO^f^ACT : 

PUBLIC  AFFAIRS  OFFICE 
U.  S.  ARMY  ENGINEER 
WATERWAYS  EXPERIMENT  STATION 
3909  HALLS  FERRY  ROAD 
VICKSBURG,  MISSISSIPPI  39100-6199 
PHONE :  (601)634-2502 


AREA  OF  RESERVATION  -  2.7  sq  km 


Waterways  Experiment  Station  Cataioging-in-Publication  Data 

Long,  Charles  E. 

Directional  wind  wave  characteristics  at  Harvest  Platform  /  by  Charles 
E.  Long  ;  prepared  for  U.S.  Army  Corps  of  Engineers. 

182  p. :  ill. ;  28  cm.  —  (Technical  report ;  CERC-95-4) 

Includes  bibliographic  references. 

1 .  Ocean  waves  —  Measurement  —  Instruments.  2.  Wind  waves  — 
Measurement.  3.  Water  waves  —  Measurement.  4.  Frequency  spectra. 
I.  United  States.  Army.  Corps  of  Engineers.  II.  U.S.  Army  Engineer 
Waterways  Experiment  Station.  III.  Coastal  Engineering  Research 
Center  (U.S.)  IV.  Title.  V.  Series:  Technical  report  (U.S.  Army  Engineer 
Waterways  Experiment  Station) ;  CERC-95-4. 

TA7  W34  no.CERC-95-4 


Contents 


Preface  . . . . .  v 

1 —  Introduction  . 1 

2—  Directional  Gauge .  4 

Gauge  Layout  .  4 

Pressure  Gauges  and  Data  Collection  Scheme .  5 

3—  Primary  Data  Analysis .  7 

Error  Checking  .  7 

Frequency-Direction  Spectra  .  12 

4—  Characterizing  Parameters .  16 

Bulk  Parameters .  16 

Frequency-Dependent  Parameters  .  21 

Summary  of  Parameters  .  24 

5—  Bulk  Parametric  Behavior .  25 

Time  Series  of  Bulk  Parameters  .  25 

Bulk  Parameter  Probability  Densities .  32 

Parameter  Correlations  with  Wave  Height .  40 

Summary  of  Bulk  Parameters .  47 

6—  Frequency-Dependent  Parameters  .  48 

Spectral  Density  and  Peak  Direction  .  48 

Circular  Moment  Parameters .  52 

Quartile  Spread  and  Asymmetry .  56 

Modal  Analysis  .  60 

Summary  of  Frequency-Dependent  Parameters .  62 

7—  Examples  of  Spectra  .  63 


8—  Characteristic  Directional  Distribution  Functions .  67 

Classification  Algorithms  .  67 

Resulting  Distributions . 69 

9—  Accuracy  of  Directional  Distribution  Estimates .  72 

Test  Method  .  72 

Test  Distributions .  75 

Test  Results  .  76 

10—  Accuracy  of  Directional  Parameters .  80 

Peak  Direction .  80 

Circular  Moment  Parameters . 83 

Quartile  Parameters .  89 

Summary  of  Parameter  Accuracy  .  91 

11—  Summary  .  94 

Gauge  and  Basic  Data  Set .  94 

Bulk  Parameters .  95 

Frequency-Dependent  Parameters  .  96 

Representative  Results . 97 

Gauge  Accuracy .  98 

Closing  Comment .  98 

References  .  100 

Appendix  A:  Time  Series  Graphs  of  Bulk  Parameters  .  A1 

Appendix  B:  Graphs  of  Example  Spectra .  B1 

Appendix  C:  Graphs  of  Characteristic  Directional  Distributions  ....  Cl 

Appendix  D:  Graphs  of  Test  Results .  D1 

Appendix  E:  Notation .  El 

SF  298 


IV 


Preface 


This  report  describes  measurement  and  results  of  analysis  of  wind  wave 
directional  distributions  at  Texaco  Oil  Company’s  Harvest  Platform,  located 
near  the  200-m  (656-ft)  isobath  about  20  km  (10.8  n.m.)  west  of  Point 
Conception,  California.  The  work  was  motivated  by  a  need  to  quantify  and 
develop  a  climatology  of  the  variability  in  natural  wind  wave  energy  distribu¬ 
tions  at  a  deepwater  site  on  the  exposed  California  coast.  A  six-element, 
high-resolution  directional  wave  gauge  was  used  for  this  purpose.  This  effort 
was  authorized  by  Headquarters,  U.S.  Army  Corps  of  Engineers 
(HQUSACE),  under  Civil  Works  Coastal  Flooding  and  Storm  Protection  Pro¬ 
gram  Research  Work  Unit  32484,  "Directionality  of  Waves  in  Shallow  Wa¬ 
ter.  "  Funds  were  provided  through  the  Coastal  Engineering  Research  Center 
(CERC),  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES),  under 
the  program  management  of  Ms.  Carolyn  M.  Holmes,  CERC.  Messrs.  John 
H.  Lockhart,  Jr.,  Charles  Chesnutt,  Barry  W.  Holliday,  and  John  F.  C.  Sanda 
were  HQUSACE  Technical  Monitors. 

This  rq)ort  was  prepared  by  Dr.  Charles  E.  Long  at  WES’s  Field  Research 
Facility  (FRF)  in  Duck,  NC,  under  the  direct  supervision  of  Mr.  William  A. 
Birkemeier,  Chief,  FRF,  and  Mr.  Thomas  W.  Richardson,  Chief,  Engineering 
Development  Division  (EDD),  CERC.  General  supervision  was  provided  by 
Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Director  and  Assis¬ 
tant  Director,  CERC,  respectively. 

Mr.  David  D.  McGehee,  Prototype  Measurement  and  Analysis  Branch, 
EDD,  CERC,  was  instrumental  in  coordinating  the  efforts  of  CERC  and  the 
State  of  California  in  gauge  installation.  Drs.  Robert  T.  Guza  and  Thomas  H. 
C.  Herbers  of  the  Center  for  Coastal  Studies,  Scripps  Institution  of  Oceanog¬ 
raphy  (SIO),  and  Dr.  Joan  M.  Oltman-Shay  of  Quest  Integrated,  Inc.,  Kent, 
WasWngton,  assisted  in  gauge  array  design.  Data  transfer  between  SIO  and 
the  FRF  was  coordinated  through  SIO’s  Coastal  Data  Information  Program, 
under  the  direction  of  Dr.  Richard  J.  Seymour,  with  particularly  helpful  assis¬ 
tance  from  Ms.  Julianna  Thomas.  The  contributions  of  all  of  these  individuals 
are  gratefully  acknowledged. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert 
W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


1  Introduction 


Long-term  observations  of  well-resolved,  deep-ocean,  wind  wave 
frequency-direction  energy  spectra  are  essential  for  advancing  knowledge  in 
ocean  surface  dynamics.  This  knowledge  is  required  for  such  things  as  pro¬ 
viding  ground  truth  to  help  interpret  satellite  imagery  of  the  ocean  surface, 
testing  evolution  and  propagation  models  of  open-ocean  wind  waves,  and 
establishing  seaward  boundary  conditions  for  models  of  wave  propagation  and 
transformation  from  deep  water  to  coastal  regions.  The  latter  problem  is  of 
considerable  interest  in  coastal  engineering  because  wind  waves  generated 
offshore  are  among  the  dominant  forcing  mechanisms  for  coastal  dynamic 
processes. 

This  problem  is  particularly  pronounced  in  the  region  of  the  Southern 
California  Bight  (Figure  1),  where  waves  generated  in  the  deep  ocean  must 
propagate  past  a  number  of  offshore  islands  to  reach  the  coast.  Refractive, 
diffractive,  and  dissipative  effects  of  the  bathymetry  around  the  islands  have 
strong  local  effects  on  the  amount  of  wave  energy  reaching  the  shore,  and 
these  effects  are  quite  sensitive  to  the  directions  of  incident  wave  propagation 
(O’Reilly  and  Guza  1991,  1993).  To  gain  insight  into  the  nature  of  direction- 
ally  distributed  wind  wave  energy  from  the  open  ocean  in  the  vicinity  of  the 
Southern  California  Bight,  a  high-resolution  directional  wave  gauge  has  been 
installed  on  a  deepwater  oil  production  platform  near  the  north  end  of  the 
bight.  Called  the  Harvest  Platform,  the  structure  is  owned  by  Texaco  Oil 
Company,  which  allowed  the  platform  to  be  used  for  scientific  purposes 
through  an  agreement  with  the  State  of  California. 

As  indicated  in  Figure  1,  Harvest  Platform  is  located  about  20  km 
(10.8  n.m.)  west  of  Point  Conception,  California,  in  water  with  a  mean  depth 
of  202  m  (663  ft).  Waves  originating  in  the  greater  Pacific  Ocean  can  reach 
the  platform  via  relatively  unobstructed  paths  from  the  north,  west,  and  south, 
and  thus  represent  conditions  on  the  seaward  side  of  the  offshore  island  chain. 
Though  the  water  depth  is  considerable,  it  is  expected  that  some  bathymetric- 
ally  induced  transformation  occurs  for  low-frequency  wind  waves  in  the  vicin¬ 
ity  of  the  platform.  With  a  nominal  frequency  range  of  0.04  to  0.3  Hz,  wind 
waves  obeying  linear  wave  theory  have  wavelengths  in  the  range  18  to  870  m 
(59  to  2,854  ft)  in  water  of  about  200-m  (656-ft)  depth.  Waves  with  lengths 
of  about  twice  the  water  depth  and  longer  will  be  influenced  by  the  ocean 
bottom.  Thus,  waves  with  lengths  of  400  to  870  m  (1,312  to  2,854  ft),  or 
frequencies  in  the  range  0.04  to  0.06  Hz,  are  expected  to  be  affected.  The 


Chapter  1  Introduction 


Figure  1 .  Southern  California  Bight  and  location  of  Harvest  Platform 

primary  effect  is  a  modest  turning  of  the  deepwater  wave  propagation  direc¬ 
tion  towards  a  bath5mietry-normal  azimuth.  This  process  can  be  modeled  with 
reasonable  simplicity  so  that  true  deepwater  conditions  can  be  estimated  from 
observations  at  Harvest  Platform. 

Because  the  platform  structure  is  stationary,  reasonably  rigid,  and  of  con¬ 
siderable  horizontal  dimensions,  it  is  nearly  ideal  for  deploying  a  spatial  array 
of  subsurface  pressure  sensors  to  establish  the  geometry  of  a  directional  wave 
gauge.  When  coupled  with  a  high-resolution  mathematical  analysis  technique, 
the  system  allows  routine  collection  of  well-resolved  wind  wave  frequency- 
direction  spectra  that  can  be  acciunulated  in  a  database  for  use  in  climatologi¬ 
cal  studies  as  well  as  the  objectives  mentioned  above. 

The  primary  purpose  of  this  report  is  to  describe  some  of  the  characteris¬ 
tics  of  directionally  distributed  wave  energy  observed  at  Harvest  Platform  in 
the  course  of  1  year.  The  directional  gauge  was  installed  and  operating  as  of 
mid-November  1992.  It  is  convenient  to  report  observations  beginning  on  1 
January  1993  and  to  include  all  observations  during  calendar  year  1993.  This 
report  additionally  includes  several  parts  related  to  the  design  and  expected 
performance  of  the  directional  gauge.  Chapter  2  describes  the  gauge  geome¬ 
try,  including  some  of  its  design  considerations  and  constraints,  and  the  data 
collection  scheme.  Chapter  3  shows  the  basic  directional  estimation  algorithm 
along  with  routine  error-checking  procedures.  Chapter  4  defines  some  param¬ 
eters  used  to  characterize  the  observations.  Chapters  5  and  6  discuss  how 
these  parameters  behave  for  bulk  spectra  and  at  individual  frequencies,  respec¬ 
tively.  Time  series  graphs  of  bulk  parameters  are  presented  in  Appendix  A. 
Representative  frequency-direction  spectra  are  discussed  in  Chapter  7,  and 
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illustrated  in  Appendix  B.  Chapter  8  discusses,  and  Appendix  C  illustrates, 
some  characteristic  directional  distribution  functions  derived  from  classifying, 
grouping,  and  averaging  observations  from  the  whole  data  set.  Chapters  9 
and  10  contain  analyses  of  directional  resolution  and  parametric  resolution, 
respectively,  based  on  the  characteristic  distributions  described  in  Chapter  8. 
Graphs  of  test  results  discussed  in  Chapter  9  are  presented  in  Appendix  D. 
Main  results  of  this  study  are  summarized  in  Chapter  11. 
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Directional  Gauge 


Gauge  Layout 

The  fundamental  detector  in  a  directional  gauge  of  the  type  used  here  is  a 
set  or  array  of  subsurface  pressure  gauges.  A  number  of  design  consider¬ 
ations  and  constraints  had  to  be  satisfied  in  constructing  the  Harvest  Platform 
array.  These  conditions  resulted  in  some  compromises  in  both  wave  frequen¬ 
cy  and  direction  resolving  ability  of  the  final  design.  For  the  spacings  (called 
lags)  between  pairs  of  gauges  to  be  tuned  to  the  wind  waves  of  interest,  one 
would  like  long  lags  for  long  waves  and  shorter  lags  for  short  waves  within 
the  wind  wave  frequency  band.  It  is  desired  to  have  the  gauges  near  the 
ocean  surface  so  that  signals  from  high-frequency  waves  are  not  overly  attenu¬ 
ated  by  depth,  and,  at  the  same  time,  have  the  gauges  deep  enough  so  that 
they  do  not  become  exposed  to  air  in  the  troughs  of  extremely  large  waves. 

Conditions  for  reasonable  wave  detection  had  to  be  balanced  against  the 
geometrical  and  logistical  constraints  of  an  operating  oil  platform.  In  plan 
view.  Harvest  Platform  has  a  rectangular  outline  of  nominal  dimensions  23  m 
(75  ft)  by  65  m  (213  ft).  These  dimensions  represent  the  largest  possible 
horizontal  gauge  spacings  in  two  orthogonal  directions.  Additional  limitations 
were  imposed  on  these  dimensions  by  the  availability  of  structural  members 
on  which  to  mount  the  gauges.  A  safe  gauge  depth  was  determined  to  be  in 
excess  of  9  m  (29.5  ft),  based  on  observations  from  prior  gauging  of  an  ex¬ 
treme  wave  height  of  about  18  m  (59  ft).'  A  depth  that  satisfied  this  con¬ 
straint  and  that  also  had  platform  structure  for  gauge  mounting  was  about 
16  m  (52.5  ft). 

In  initial  array  design,  it  was  determined  that  nine  gauges  would  provide 
the  best  wave  response  by  establishing  lag  spacings  that  increased  in  reason¬ 
ably  uniform  increments  within  the  spatial  constraints  of  the  platform.  Addi¬ 
tionally,  it  was  desired  to  establish  redundant  lags  to  ensure  that  wave  data 
were  not  being  compromised  by  passing  around  and  through  the  open  struc¬ 
ture  of  Harvest  Platform  itself.  Funding  limitations  reduced  to  six  the  number 


'  Personal  communication,  1991,  Dr.  R.  J.  Seymour,  Coastal  Data  Information  Program 
(CDIP),  Scripps  Institution  of  Oceanography  (SIO). 
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of  gauges  that  could  be  purchased,  installed,  and  maintained  on  a  regular 
basis. 

The  array  design  that  evolved  from  these  constraints,  and  the  one  installed 
on  Harvest  Platform,  is  shown  in  plan  view  in  Figure  2.  This  figure  also 
shows  the  gauge  numbering  scheme,  the  orientation  of  the  array  relative  to 
geophysical  coordinates,  and  the  direction  convention  used  in  this  report.  The 
outside  dimensions  of  the  array  are  22.7  m  (74.5  ft)  by  59.4  m  (194.9  ft), 
which  utilizes  almost  the  whole  of  the  platform  dimensions.  Because  the 
gauges  are  mounted  on  a  fixed  structure,  their  positions  are  known  to  within  a 
few  centimeters.  Note  that  the  regular  spacing  of  the  gauges  provides  a  num¬ 
ber  of  redundant  lags  for  testing  platform  interference.  Specifically,  cross 
spectra  of  wave  signals  between  gauge  pairs  1-4,  2-3,  and  5-6  should  be  iden¬ 
tical,  as  should  cross  spectra  between  gauge  pairs  1-2  and  4-3,  as  well  as  1-5 
and  4-6.  An  example  of  such  cross  spectra  computed  from  measured  wave 
signals  is  shown  in  Chapter  3. 

All  six  gauges  were  mounted  at  15.72  m  (51.57  ft)  below  mean  sea  level, 
which  satisfies  the  safe  depth  constraint.  Surface  correction  of  wave  signals 
from  this  depth  using  the  linear  wave  pressure  response  function  {Shore  Pro¬ 
tection  Manual  1984)  is  quite  conservative  in  that  the  signal-to-noise  ratio  for 
these  gauges  is  quite  high  (t)q)ically  order  100)  at  the  highest  resolvable  wave 
frequency.  This  highest  frequency  is  established  not  by  frequency  response 
limitations,  but  by  the  shortest  lag  spacing  of  the  array.  To  avoid  aliasing  in 
directional  estimation,  the  shortest  wavelength  that  can  be  resolved  is  twice 
the  shortest  lag  spacing,  or  45.4  m  (149.0  ft)  for  this  array.  This  wavelength 
corresponds  to  a  wave  frequency  of  about  0.18  Hz.  To  be  conservatively 
clear  of  aliasing  effects,  signal  analysis  used  in  this  report  was  truncated  at 
0.16  Hz  for  directional  estimation. 

Though  the  final  installed  array  was  less  than  theoretically  ideal,  its  re¬ 
sponse  and  direction-resolving  ability  are  quite  reasonable.  Quantification  of 
directional  resolution  for  this  array  is  presented  in  Chapters  9  and  10  of  this 
report. 


Pressure  Gauges  and  Data  Collection  Scheme 


Individual  sensors  were  Model  TIE  absolute  pressure  sensors  manufactured 
by  Sensotec  Transducer  Company  with  operating  ranges  of  0  to  100  psia  (0  to 
689.5  kPa),  and  a  manufacturer’s  stated  accuracy  of  ±0.1  percent  of  full 
scale.  The  set  of  six  gauges  on  Harvest  Platform  were  sampled  simul¬ 
taneously  at  1  Hz,  digitized,  then  fed  to  a  concentrator  where  the  set  of  sam¬ 
ples  was  buffered.  Buffered  signals  were  periodically  transmitted  to  shore 
through  a  telephone  connection,  and  ultimately  stored  on  the  main  computer 
of  the  CDIP,  Ocean  Engineering  Research  Group,  Marine  Research  Division, 
SIO,  La  Jolla,  CA. 

Sample  sets  known  as  collections  consist  of  time  series  of  data  from  each 
of  the  six  gauges.  Each  time  series  is  8,192  sec  (2  hr  16  min  32  sec)  in^ 
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Figure  2.  Dimensions  and  orientation  of  the  Harvest  Platform  array 

length.  Collections  were  made  eight  times  daily,  at  approximately  3-hr  inter¬ 
vals.  From  1  January  to  2  July  1993,  nominal  collection  start  times  were 
0100,  0400,  0700,  1000,  1300,  1600,  1900,  and  2200  Greenwich  Mean  Time 
(GMT).  From  3  July  to  31  December  1993,  nominal  collection  start  times 
were  0200,  0500,  0800,  1100,  1400,  1700,  2000,  and  2300  GMT.  Actual 
collection  start  times  varied  by  several  minutes  on  either  side  of  a  given  nomi¬ 
nal  start  time  because  the  amount  of  time  required  to  establish  a  phone  link 
varied  from  collection  to  collection. 

Data  processing  for  results  presented  in  this  report  was  not  performed  at 
the  CDIP  site,  and  so  is  independent  of  the  processing  done  and  published  by 
that  group  (CDIP  Monthly  Reports).  Data  collections  were  transferred  to  the 
Field  Research  Facility  (FRF)  of  the  U.S.  Army  Engineer  Waterways  Experi¬ 
ment  Station’s  Coastal  Engineering  Research  Center  for  processing  by  high- 
resolution  techniques  that  are  different  from  those  used  by  CDIP.  Initially, 
data  transfer  was  achieved  by  way  of  magnetic  tape,  but  later  was  accom¬ 
plished  over  an  electronic  network. 

All  possible  collections  are  not  represented  in  the  data  set  used  for  this 
report.  A  number  of  collections  were  lost  because  of  the  inability  to  establish 
or  maintain  electrically  clean  phone  links  to  the  concentrator  on  Harvest  Plat¬ 
form.  An  additional  number  of  collections  were  not  processed  because  data 
did  not  satisfy  error-checking  constraints  established  at  the  FRF  analysis  site. 
These  error-checking  methods  are  discussed  in  Chapter  3  of  this  report.  Of 
the  maximum  possible  2,920  collections  during  calendar  year  1993,  a  total  of 
2,339  collections  were  acquired  and  processed. 


6 


Chapter  2  Directional  Gauge 


3  Primary  Data  Analysis 


For  logistically  successful  collections,  initial  data  processing  was  done  by 
checking  data  quality  through  a  series  of  gauge  intercomparisons,  and,  for 
data  of  sufficient  quality,  computing  frequency-direction  spectra  as  primary 
results.  All  of  these  steps  rely  on  Fourier  analysis  of  pressure  gauge  time 
series  data,  and  subsequent  computation  of  cross-spectral  densities.  A  discus¬ 
sion  of  error-checking  procedures  then  leads  logically  to  the  subsequent  steps 
involved  in  frequency-direction  spectral  computation. 


Error  Checking 

The  first  step  in  data  processing  is  computation  of  discrete  estimates  of 
frequency  autospectra  of  pressure  signals,  and  surface-corrected  cross-spectral 
densities  of  signals  from  all  pairs  of  gauges.  Cross  spectra  are  denoted 
-  iQg(f„),  where  C..(/„)  is  the  coincident  spectrum,  is  the 

quadrature  spectrum,  i  and  j  are  indices  ranging  in  value  from  1  to  6  that 
refer  to  the  gauge  numbers  shown  in  Figure  2,  and  is  the  n*  of  a  set  of  N 
discrete  frequencies.*  Frequency  autospectra  are  denoted  5(/„),  and,  if  sur¬ 
face  corrected,  are  identically  equal  to  .  All  spectra  are  computed 

using  Welch’s  method  (Welch  1967)  with  standard  Fourier  analysis  techniques 
(Bendat  and  Piersol  1971). 

In  a  collection,  the  8,192-sec  time  series  from  each  gauge  is  analyzed  in  15 
half-lapped  segments  of  1,024  sec  duration.  Each  segment  is  demeaned,  ta¬ 
pered  with  a  variance-preserving  window,  and  converted  to  the  frequency 
domain  with  a  discrete  Fourier  transform.  At  this  point,  the  analysis  is  split 
into  two  parts:  estimates  of  pressure  autospectra  from  each  gauge  at  depth, 
and  estimates  of  surface-corrected  cross  spectra  of  sea  surface  displacement. 
For  cross-spectral  estimates,  each  transformed  segment  of  pressure  data  is 
corrected  to  a  sea  surface  displacement  transform  using  the  linear  wave  pres¬ 
sure  response  function  (Dean  and  Dalrymple  1984).  Raw  cross-spectral 
estimates  are  formed  from  temporally  corresponding  transform  segments 
among  all  pairs  of  gauges.  Raw  autospectral  estimates  are  formed  from  each 
of  the  15  transform  segments  for  each  individual  gauge. 


'  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  E). 
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For  both  autospectra  and  cross  spectra,  smooth  estimates  are  formed  by 
averaging  raw  estimates  over  all  15  segments,  and  further  smoothing  is 
achieved  by  averaging  results  over  10  adjacent  frequency  bands.  Final  reso¬ 
lution  frequency  bandwidth  df=  0.00977  Hz,  and,  because  of  the  spatial 
Nyquist  limits  discussed  in  Chapter  2,  the  wind  wave  pass  band  of  frequencies 
ranges  from  0.044  to  0.162  Hz,  which  corresponds  to  (A^=)  13  discrete  fre¬ 
quency  bands.  Degrees  of  freedom  for  spectral  estimates  range  from  160  to 
about  200,  depending  on  the  extent  to  which  the  second  halves  of  time  series 
segments  are  correlated  with  the  first  halves  (Welch  1967). 


Autospectral  intercomparisons 

One  part  of  error  checking  is  graphic  intercomparisons  among  the  six 
pressure  gauges  of  signal  means  and  autospectra,  an  example  of  which  is 
shown  in  the  lower  left  graph  of  Figure  3.  Frequency  autospectral  estimates 
of  data  from  all  six  pressure  gauges  are  plotted  on  the  same  set  of  axes  from 
the  first  resolvable  frequency  band  out  to  the  temporal  Nyquist  frequency 
(0.5  Hz,  half  the  1-Hz  sampling  frequency).  As  can  be  seen,  the  most  ener¬ 
getic  parts  of  the  spectra  are  in  the  analysis  pass  band  of  frequencies.  At  the 
highest  frequencies  shown  in  this  graph  are  the  spectra  of  what  is  dominantly 
system  and  data  transmission  noise,  called  the  noise  floor.  If  a  pressure  gauge 
is  malfunctioning,  its  autospectrum  will  deviate  from  the  main  group  of 
curves,  and  become  immediately  obvious.  A  deviate  noise  floor  from  one  of 
the  gauges  indicates  spiky  or  jumpy  data,  another  cause  for  concern. 

The  small  inset  graph  in  the  lower  left  graph  of  Figure  3  is  an  analysis  of 
signal  means.  The  closely  packed,  nearly  constant  group  of  symbols  repre¬ 
sents  the  deviations  of  the  segment  means  from  the  median  of  the  set  of  seg¬ 
ment  means  for  each  of  the  15  segments.  If  a  gauge  develops  significant  drift 
problems,  it  will  become  immediately  obvious  as  a  symbol  that  deviates  from 
the  main  group  of  symbols.  As  shown.  Figure  3  is  somewhat  difficult  to 
read,  but  it  has  been  reduced  in  this  report  for  publication  purposes,  and  is 
much  easier  to  read  in  its  original  size.  The  triangle  symbols  in  the  small 
inset  graph  show  the  deviation  of  the  indicated  water  surface  from  mean  sea 
level  (gauge  height  off  the  bottom  plus  median  of  gauge  mean  depths  for  each 
segment  minus  the  total  long-term  mean  ocean  depth  of  202  m),  and  is  there¬ 
fore  an  indication  of  tide  stage  at  Harvest  Platform  for  each  of  the  15  seg¬ 
ments  in  a  collection. 


Coherence  and  phase  comparisons 

The  next  step  in  general  analysis,  and  one  useful  for  error  checking,  is  the 
computation  of  a  dimensionless  cross  spectrum  which  can  be  ex¬ 

pressed  as  a  complex,  Hermitian  matrix  in  terms  of  dimensional  cross  spectra 
by 
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Figure  3.  Autospectral  intercomparison  and  frequency-direction  spectral 
estimate 


Equation  1  is  used  in  directional  spectral  estimation,  but  its  use  in  error 
checking  is  in  its  relation  to  coherence  and  phase  estimates.  Coherence  of 
signals  from  gauges  i  and  j  at  discrete  frequency  is  defined  in  terms  of 
the  dimensionless  cross  spectrum  as 
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^Wn)  = 


(2) 


10 


Signal  phase  difference  of  gauge  i  relative  to  gauge  j  at  frequency  /„  is 


=  tan-' 


(3) 


where  Re[  ]  and  Im[  ]  are  the  real  and  imaginary  parts,  respectively,  of  the 
entity  contained  in  square  brackets. 


As  mentioned  in  Chapter  2,  if  multiple  pairs  of  gauges  are  arranged  so  that 
each  pair  is  separated  by  the  same  distance  and  has  the  same  directional  orien¬ 
tation  from  one  gauge  to  the  second  gauge,  then  in  a  uniform  wave  field, 
cross  spectra  of  the  pairs  of  gauges  should  be  identical.  In  the  Harvest  Plat¬ 
form  array  there  are  several  such  sets  of  pairs  as  can  be  seen  in  Figure  2. 
Specifically,  in  terms  of  coherences  and  phases  defined  by  Equations  2  and  3, 
one  would  expect 


r„(/,)  -  r„(/.)  -  r„tt) 


=  •f'M)  =  <“> 


as  well  as 


r,,!/.)  -  r„(/,)  ■I’M)  •  ■I’M.) 


and 


=  K(fn) 


(5) 


(6) 


Figure  4  is  an  example  of  coherence  and  phase  error  checking,  showing 
graphs  of  the  seven  coherence  and  phase  functions  named  in  Equations  4, 

5,  and  6  (upper,  middle,  and  lower  sets  of  graphs  in  Figure  4,  respectively). 
For  any  set  of  wave  approach  directions,  some  of  the  gauges  are  on  the  up- 
wave  side  of  the  platform,  and  others  are  on  the  down-wave  side.  If  the 
platform  legs,  risers,  and  other  structural  members  caused  significant  distor¬ 
tion  of  the  wave  field,  there  would  be  clear  disagreement  among  the  various 
curves  illustrated  in  Figure  4.  No  such  disagreement  is  obvious  in  Figure  4, 
which  is  exemplary  of  acceptable  collections,  so  it  is  concluded  that  Harvest 
Platform  does  not  interfere  in  a  significant  way  with  wind  waves  as  observed 
by  these  gauges. 

An  occasional  problem  occurs  in  the  data  collection  procedure  that  is  also 
detectable  with  coherence  and  phase  plots.  Sometimes  a  data  point  is  dropped 
during  telephone  transmission  from  the  data  buffer.  When  this  happens,  the 
time  series  of  data  for  the  affected  gauge  is  effectively  shifted  by  one  time 
step  relative  to  data  from  the  other  gauges.  If  not  caught,  a  lost  data  point 
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Harvest  Platform 
1553  20  Jun  93  Gages:  123456 

coherence  phase  (deg) 


Figure  4.  Sample  coherence  and  phase  function  comparisons 


can  have  devastating  effects  on  directional  spectral  estimation.  However,  such 
a  data  loss  is  readily  apparent  in  a  graphic  display  like  Figure  4.  At  a  1-Hz 
sample  rate,  a  single  lost  data  point  early  in  a  collection  induces  a  phase  shift 
of  about  15  deg  for  0.04-Hz  waves  and  about  60  deg  for  0.16-Hz  waves.  If 
the  phase  function  curves  like  those  in  Figure  4  deviate  from  each  other  by 
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such  phase  shifts,  it  is  quite  apparent.  If  a  data  point  is  lost  later  in  a  collec¬ 
tion,  the  evident  phase  shifts  are  relatively  less  because  only  part  of  a  time 
series  is  then  out  of  phase  with  other  gauge  records.  If  a  data  point  is  lost 
late  in  a  collection,  it  is  difficult  to  detect  in  graphs  like  those  of  Figure  4,  but 
the  effect  of  the  lost  point  is  greatly  reduced  because  most  of  the  time  series 
are  in  phase. 

The  combined  effects  of  intercomparing  frequency  autospectra  and  coher¬ 
ence  and  phase  functions  for  the  pressure  gauge  array  on  Harvest  Platform 
provide  clear  indications  of  faulty  or  suspect  data.  When  such  conditions  are 
detected  in  a  collection,  frequency-direction  spectra  are  not  computed.  Such 
rigorous  examination  of  the  data  ensures  that  only  high-quality  time  series  are 
used  in  directional  estimation. 


Frequency-Direction  Spectra 


Estimates  of  frequency-direction  spectra  are  made  using  the  iterative  maxi¬ 
mum  likelihood  estimator  (IMLE)  developed  by  Pawka  (1983).  The  first  step 
in  this  method  is  a  computation  of  the  maximum  likelihood  estimate  (MLE) 
described  by  Davis  and  Regier  (1977).  All  estimates  determine  directional 
distribution  functions  on  a  frequency-by-frequency  basis.  At  a  given  frequen¬ 
cy,  a  directional  distribution  function  is  related  to  its  corresponding 

frequency-direction  spectrum  S(f^,6^)  by 


S(fn) 


(7) 


where  6^  is  a  discrete  angle  indicating  the  direction  from  which  wave  energy 
arrives,  measured  counterclockwise  from  true  north,  as  indicated  in  Figure  2. 
The  direction  index  m  ranges  from  /M=lt0  7n  =  Af=181,  while  direction 
ranges  from  =  -180  deg  to  0,g,  =  180  deg  in  steps  of  dd  =  2  deg.  The 
normalizing  term  in  the  denominator  on  the  right-hand  side  of  Equation  7  is 
the  sea  surface  displacement  variance  spectral  density.  Because  the  integral 
with  respect  to  direction  of  S(f^,0^)  is  identically  S(f^),  the  directional  distri¬ 
bution  function  has  the  property 

i:D(f„,ejde  =  I  (8) 


for  all/„. 
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The  lowest  order  estimate  is  the  MLE,  which  takes  the  form 


Do(fM  =  — r-7 - ^ -  (9) 

(=1  >1 

where  is  a  factor  of  order  1  that  is  used  to  ensure  Equation  8  is  satisfied,  I 
is  the  number  of  gauges,  the  My^(f^)  are  elements  of  the  inverse  of  the  di¬ 
mensionless  cross-spectral  matrix  defined  by  Equation  1,  k^(6^)  is  the  wave 
number  vector,  and  x.  and  Xj  are  coordinate  position  vectors  of  gauges  i  and 

j  ,  respectively.  The  wave  number  vector  k^{d^)  is  a  function  of  both  fre¬ 
quency  and  direction,  being  defined  as 

X„(0„)  =  k^cosd  e  +  k  s'md  e 

m'  n  m  x  n  my 

where  and  are  spatial  coordinate  unit  vectors  in  the  x-  and  y -directions, 
respectively,  and  is  wave  number  vector  magnitude,  which  is  related  with 
gravitational  acceleration  g  to  frequency  and  water  depth  d  through  the 
linear  wave  dispersion  relation 

=  gkjdiihk^d 


The  need  for  the  iterative  part  of  the  estimator  arises  because  the  zeroth 
order  estimate  of  Equation  9  can  be  shown  to  produce  directional  distributions 
that  are  too  directionally  broad  relative  to  true  directional  distributions.  An 
IMLE  result  is  achieved  in  several  computational  steps.  At  the  r*  iteration 
step,  an  estimate  of  the  observed  cross-spectral  matrix  is 

computed  from  the  previous  directional  distribution  function  estimate 

by 


A  new  intermediate  directional  distribution  function  estimate  D'  (f^,  6^)  is 
computed  using  the  cross-spectral  matrix  of  Equation  12  by 
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where  a^,  like  a^,  is  adjusted  so  that  Equation  8  is  satisfied  for  £)/  A 

correction  is  computed  for  Z)/  by  first  computing 


\(fn>V  =  1  - 


DoifnA) 


(14) 


and  then  finding  a  new  directional  distribution  function  estimate 
from 


DrifnA)  -Di  if M 


1  + 


y\(fnA) 


(15) 


The  parameters  /3  and  7  in  Equation  15  control  the  rate  of  convergence  of  the 
estimator.  As  used  by  Pawka  (1983),  the  values  j8  =  1  and  7  =  5  were  used 
for  all  estimates  discussed  in  this  report. 

An  objective  of  the  iterative  technique  is  to  try  to  make  estimates  of  the 
cross  spectrum  computed  with  Equation  12  as  close  as  possible  to  the  mea¬ 
sured  cross  spectrum  of  Equation  1 .  In  each  iterative  loop,  a  convergence 
check  is  computed  as  the  sum  of  the  squared  magnitudes  of  the  differences 
of  the  two  cross  spectra.  This  takes  the  form 


>1 


(16) 


Iteration  continues  as  long  as  decreases  between  successive  iterations,  or 
until  an  upper  limit  R  of  iterations  has  been  completed.  In  computations 
reported  herein,  jR  =  30  because  tests  with  synthetic  data  have  indicated  that 
very  little  increase  in  accuracy  is  generally  obtained  if  a  higher  iteration  limit 
is  used. 

Equations  9  to  16  form  the  basis  of  the  IMLE  technique.  For  a  given 
frequency,  and  for  the  iteration  r  that  satisfies  the  convergence  check,  the 
dimensional  directional  distribution  of  energy  is  formed  from 

S(fM  -  (1' 


When  evaluated  for  all  frequencies  in  the  wind  wave  pass  band,  the  resulting 
entity  is  the  frequency-direction  spectrum  for  a  given  collection.  An  example 
of  such  a  spectrum  is  illustrated  in  Figure  3.  The  upper  graph  is  a  three- 
dimensional  plot  of  S(f^,d^),  and  the  lower  right  graph  is  a  contour  plot  of 
the  spectrum.  The  right  panel  in  the  three-dimensional  plot  is  a  linear  graph 
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of  the  discrete  frequency  spectrum  S(f^),  which  is  related  to  the  frequency- 
direction  spectrum  by 

s(f„)  =  E  sifjjde 

m-l 

The  left  panel  in  the  three-dimensional  plot  is  a  linear  graph  of  the  direction 
spectrum  S{dJ ,  which  is  the  directional  analog  of  the  frequency  spectrum, 
defined  by 

n=l 

Figure  3  illustrates  a  case  having  a  moderate  amount  of  energy,  with  a 
low-frequency  peak  of  wave  energy  arriving  from  slightly  West  of  south,  and 
a  second  mode  of  energy  having  greater  spread  in  both  frequency  and  direc¬ 
tion  arriving  from  the  northwest  quadrant.  To  classify  this  spectrum,  as  well 
as  all  the  spectra  in  the  1993  year  of  observations  at  Harvest  Platform,  a  set 
of  characterizing  parameters  is  defined  in  Chapter  4.  The  behavior  and  rela¬ 
tionships  among  those  parameters,  discussed  in  subsequent  chapters,  then 
reveal  some  of  the  directional  characteristics  at  this  site. 
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4  Characterizing  Parameters 


Two  groups  of  parameters  are  used  to  characterize  the  members  of  the  set 
of  frequency-direction  spectra  used  in  this  study.  One  set,  called  bulk  param¬ 
eters,  atten:q)ts  to  characterize  whole  spectra  (like  the  three-dimensional  plot  in 
Figure  3)  so  that  common  features  and  correlations  among  various  properties 
can  be  identified.  The  other  set  of  parameters  characterizes  components  of 
frequency-direction  spectra  on  a  frequency-by-frequency  basis.  Called 
frequency-dependent  parameters,  they  characterize  directional  properties  that 
can  be  identified  with  low-,  medium-,  and  high-frequency  parts  of  the  ob¬ 
served  spectra.  This  chapter  defines  the  parameters  used.  Parametric  behav¬ 
ior  is  described  in  Chapters  5  and  6. 


Bulk  Parameters 

A  set  of  nine  parameters,  along  with  total  water  depth  and  modal  analyses 
of  both  frequency  and  direction  spectra,  are  defined  to  help  characterize  com¬ 
plete  frequency-direction  spectra.  The  nine  parameters  are  characteristic  wave 
height,  peak  frequency,  two  measures  of  characteristic  direction,  two  mea¬ 
sures  of  directional  spread,  two  measures  of  as)mimetry  of  directionally  dis¬ 
tributed  wave  energy,  and  a  measure  of  kurtosis  of  directional  distributions. 
This  section  establishes  the  mathematical  definitions  of  these  parameters. 


Wave  height 


In  this  report,  the  conventional  definition  of  characteristic  wave  height 
based  on  4  times  sea  surface  displacement  variance  is  used.  Denoted  ,  it 
can  be  defined  in  terms  of  the  full  frequency-direction  spectrum,  the  frequen¬ 
cy  spectrum  defined  by  Equation  18,  or  the  direction  spectrum  defined  by 
Equation  19.  A  definition  that  relates  all  of  these  entities  is 


Hi 

16 


M  N 


M 


■  E  E  -  E  w.)4r-  E 


ms=l  /i=l 


n-l 


(20) 
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Peak  frequency 


Peak  frequency  is  defined  as  the  discrete  frequency  at  which  the  fre¬ 
quency  spectrum  5(/„)  is  maximum.  This  definition  is  also  conventional,  in 
that  it  is  the  usual  characteristic  firequency  defined  for  nondirectional  gauges. 
It  represents  the  frequency  that  has  the  most  overall  wave  energy,  and  is 
subject  to  less  statistical  estimation  uncertainty  than  the  most  logical  alterna¬ 
tive  peak  firequency  definition,  which  is  the  frequency  at  which  S(f^,  6^)  is 
maximum. 


Peak  direction 

There  are  several  ways  to  designate  a  characteristic  direction  for  a  distribu¬ 
tion  of  energy.  Two  ways  are  used  in  this  report.  One  is  based  on  a  circular 
moment  of  the  distribution,  and  is  defined  in  the  following  section.  A  simpler 
one,  but  one  that  is  subject  to  considerable  statistical  imcertainty,  is  a  point  of 
maximum  variance  density  in  a  directional  distribution.  In  this  report,  a  peak 
direction  is  defined  as  the  direction  of  maximum  variance  density  in  the 
directional  distribution  associated  with  the  peak  frequency.  In  symbols,  6^  is 
the  direction  at  which  S(f^,d^)  is  a  maximum.  It  is  interpreted  as  the  direc¬ 
tion  of  the  most  energetic  waves  at  the  frequency  containing  the  greatest  over¬ 
all  energy. 


Parameters  from  circular  moments 

A  set  of  parameters  that  define  a  characteristic  mean  direction,  as  well  as 
measures  of  directional  spread,  skewness,  and  kurtosis  was  proposed  by  Kuik, 
van  Vledder,  and  Holthuijsen  (1988)  based  on  theoretical  and  analogical  argu¬ 
ments  by  Borgman  (1969)  and  Mardia  (1972)  involving  circular  moments. 
Though  derived  for  directional  distributions  at  individual  frequencies,  the 
definitions  can  be  applied  to  any  directional  distribution  function.  For  the 
purposes  of  characterizing  a  frequency-direction  spectrum  as  a  whole,  the 
directional  spectrum  S(d^),  as  defined  by  Equation  19,  can  be  used  because  it 
represents  total  wave  energy  in  any  given  direction  arc. 

To  define  a  proper  directional  distribution  function  (one  that  integrates  to 
unit  area),  S(d^)  must  be  normalized  by  its  own  area.  By  Equation  20,  this 

area  is  identically  —H^,  so  the  appropriate  directional  distribution  function  is 
16 

Diej  =  -^S(ej  m  =  l,2,...,M  (21) 

Hmo 


Discrete  approximations  of  the  circular  moments  defined  by  Kuik, 
van  Vledder,  and  Holthuijsen  (1988)  in  terms  of  D(d^)  take  the  forms 


Chapter  4  Characterizing  Parameters 


m,  =  E  cos(«.-«.)D(«.)<i9 


Bij  =  E  cos(29„-29n)Z)(9„)ti9 


n,  =  £  sm(29,-29,)D(«.)<(9 


where  is  a  characteristic  direction  called  the  mean  Erection  defined  by 
requiring  n,  =  0 .  With  this  constraint.  Equation  23  can  be  solved  for  mean 
direction  with  the  result 


0o  =  tan-‘  - 

E  cosd^Didjde 


With  6^  determined  by  Equation  26,  the  remaining  circular  moments  m, ,  m^, 
and  Wj  can  be  computed  from  Equations  22,  24,  and  25,  respectively. 

Kuik,  van  Vledder,  and  Holthuijsen  (1988)  then  define  a  measure  of  direc¬ 
tional  spread  called  the  circular  width  a  as 

a  =  (2  -  2m,  <27; 

as  well  as  a  measure  of  asymmetry  of  a  directional  distribution  called  the 
circular  skewness  7,  given  by 


(1  - 


and  a  measure  of  the  flatness  of  a  directional  distribution  called  circular  kur- 
tosis  d ,  which  is 


18 


Chapter  4  Characterizing  Parameters 


(29) 


6  -  8fft,  +  2m^ 

(2  -  2m,  f 


The  circular  parameters  6^,  a,  y,  and  5  are  intended  to  be  directional  ana¬ 
logs  of  parameters  derived  from  moments  of  linear  functions  like  the  frequen¬ 
cy  spectrum.  Kuik,  van  Vledder,  and  Holthuijsen  (1988)  found  the  analogy  to 
be  quite  good  for  narrow  directional  distributions  of  wave  energy,  and  reason¬ 
able  for  a  variety  of  broader  distributions. 


Quartile  parameters 

Two  parameters  that  are  modestly  more  intuitive  than  the  corresponding 
circular  parameters,  and  are  also  useful  for  characterizing  spread  and  asymme¬ 
try  in  directional  distribution  functions  are  the  quartile  spread  A6  and  quartile 
asymmetry  A  used  by  Long  and  Oltman-Shay  (1991).  The  concept  is  based 
on  the  fact  that  any  directional  distribution  function  integrates  across  all  direc¬ 
tions  to  unity,  and,  for  the  Harvest  Platform  array,  spans  a  full  360-deg  view¬ 
ing  arc.  An  integral  of  the  discrete  distribution  function  from  the  direction  of 
minimum  energy  6^  (where  m^^  is  the  discrete  direction  index  at  which  the 
energy  minimum  occurs)  to  any  arbitrary  angle  creates  a  function  6^^) 

that  increases  monotonically  from  zero  to  an  upper  limit  of  unity.  The  direc¬ 
tions  at  which  this  integral  (interpolated  as  necessary  from  discrete  data)  has 

the  values  i,  and  -  are  the  first  quartile,  median,  and  third  quartile  direc- 

4  2  4 

tions  of  the  directional  distribution,  respectively.  Differences  among  these 
directions  then  provide  information  about  the  spread  and  asymmetry  of  the 
distribution. 

For  bulk  parameterization,  the  direction  spectrum  D(d^)  is  used  as  a 
representative  distribution  function.  The  integral  function  is  defined  by 


'  E  0(9, )<« 


(30) 


and  the  quartile  directions  satisfy  the  equations 
«9»,  -9..)  -  0.25 

((9»,-e.J  =0.50 


(31) 


(32) 


and 
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where,  if  necessary,  a  distribution  is  treated  as  cyclic  because  the  direction 
axis  covers  a  complete  circle. 


A  measure  of  directional  spread  A6  is 


Because  the  range  of  directions  in  Equation  34  spans  the  two  middle  quartiles, 
this  definition  of  directional  spread  has  the  specific  interpretation  that  it  is  the 
arc  subtending  the  central  50  percent  of  the  energy  distribution. 

A  measure  of  asymmetry  of  a  distribution  is  the  ratio  of  the  directional 
width  of  the  third  quartile  to  that  of  the  second  quartile.  By  taking  the  natural 
logarithm  of  this  ratio,  a  symmetric  distribution  has  an  asymmetry  parameter 
A  near  zero,  and  that  for  a  skewed  distribution  acquires  a  positive  or  negative 
sign  if  the  skewness  is  toward  larger  or  smaller  angles,  respectively.  The 
asymmetry  parameter  is  thus  defined  as 


A 


In 


^75%  ^50% 

^50*  ~  ^25* 


(35) 


Modal  analysis 

The  number  and  intensity  of  energy  peaks,  or  modes,  indicate  the  number 
and  relative  importance  of  various  wave  source  processes.  For  bulk  charac¬ 
terization  of  observed  spectra,  a  modal  analysis  is  done  both  for  directional 
modes,  using  D(6^)  as  a  characterizing  function,  and  for  fi-equency  modes, 

using  the  normalized  frequency  spectrum  l6S(fJIHlo  as  a  characterizing 
function. 

Usually,  a  mode  is  distinguished  as  a  local  maximum  in  a  distribution  with 
modal  boundaries  identified  as  the  adjacent  minimal  points.  However,  be¬ 
cause  the  IMLE  method  can  generate  modestly  spurious  maxima,  especially  in 
the  estimation  of  broad,  flat  true  distributions,  a  governing  condition  must  be 
imposed  to  ensure  the  distinction  of  true  modes.  In  this  report,  a  modal  edge 
is  defined  as  a  minimum  point  for  which  a  directional  distribution  has  a  value 
that  is  no  greater  than  half  the  value  of  the  lesser  of  the  adjacent  maximum 
points.  For  consistency,  the  same  definition  is  used  to  analyze  normalized 
frequency  spectra.  In  both  types  of  analysis,  the  number  of  modes  is  one  less 
than  the  number  of  modal  edges  found  in  a  distribution. 
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The  area  of  a  distribution  between  modal  edges  is  a  fraction  less  than  or 
equal  to  unity,  which,  when  multiplied  by  100  percent,  indicates  the  percent 
of  total  energy  represented  by  that  mode.  By  recording  the  percentage  of 
energy  in  the  modes  of  a  distribution,  it  can  be  determined  if  secondary  and 
higher  modes  contain  significant  amoimts  of  energy,  or  are  just  relatively 
minor  fluctuations  in  the  low-energy  tails  of  a  distribution.  In  both  direction 
and  frequency  modal  analyses,  percentages  of  energy  in  all  modes  are  comput¬ 
ed. 


Frequency-Dependent  Parameters 

Seven  parameters  are  used  to  characterize  directional  distribution  functions 
at  each  frequency  for  each  collection.  The  parametric  definitions  are  essen¬ 
tially  the  same  as  those  used  for  bulk  analysis  except  directional  distribution 
functions  at  each  frequency  are  used  as  basis  functions.  In  addi¬ 

tion,  a  modal  analysis  is  done  for  each  distribution.  This  section  introduces 
notation  and  provides  definitions  of  these  parameters.  Note  that  because  the 
basic  analysis  results  in  13  discrete  frequency  bands,  there  are  13  frequency- 
dependent  parameters  derived  from  each  collection. 


Peak  direction 

Using  n  in  subscript  form  to  denote  the  index  of  a  discrete  frequency  /^, 
the  peak  direction  6^^  of  a  directional  distribution  function  at  that  frequency 
D(f^,0^)  is  defined  as  the  direction  at  which  D(f^,6^)  is  a  maximum. 

Parameters  from  circular  moments 

As  in  bulk  characterization,  there  are  parameters  defining  mean  direction, 
circular  width,  circular  skewness,  and  circular  kurtosis  for  each  of  the  13 
directional  distributions  in  a  collection.  The  forms  of  the  definitions  of  these 
parameters  are  the  same  as  for  bulk  parameters,  but  the  basis  distributions  are 
different,  and  the  notation  is  adapted  to  the  fact  that  there  is  a  set  of  parame¬ 
ters  for  each  frequency. 

The  four  basic  circular  moments  at  the  n*  frequency  are 


M 


<3«) 


m=l 


M 


(37) 


m=l 
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n^.£sin(29.-29^)D(/„9,)i9  O' 

where  6^  is  the  mean  direction  at  frequency  defined  by  requiring  «,„  =  0 . 
Solving  for  6^^  using  Equation  37  leads  to 


E  sm9.D(/.,«.)rf9 

E  <^KD(f.,K)d» 


which  can  be  used  in  Equations  36,  38,  and  39  to  solve  for  m^,  and 
respectively. 

Again  following  the  work  by  Kuik,  van  Vledder,  and  Holthuijsen  (1988),  a 
frequency-dependent  circular  width  is 

-  (2  -  2'",.)“ 

a  circular  skewness  is 


(i  - 


and  a  circular  kurtosis  5_  is 


^  _  6-Sm^„  +  2m^ 
"■  (2-2/n,f 


It  is  noted  that  the  parameters  6^,  y^,  and  derived  from  directional 

distributions  at  individual  frequencies  are  the  types  discussed  by  Kuik,  van 
Vledder,  and  Holthuijsen  (1988)  in  their  work,  so  these  frequency-dependent 
parameters  may  be  more  comparable  to  their  results  than  the  corresponding 
bulk  parameters  defined  in  the  previous  section. 
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Quartile  parameters 

As  with  bulk  parameters,  frequency-dependent  parameters  for  quartile 
spread  and  quartile  asymmetry  are  useful  in  characterizing  the  direc¬ 
tional  distributions  D(/„,  0^) .  The  forms  of  their  definitions  follow  Equations 
30  to  35,  and  their  meaning  and  interpretation  are  the  same.  For  each  fre¬ 
quency,  an  integral  function  )  can  be  defined  as 


(44) 


where  m  .  is  understood  to  be  the  direction  index  of  the  minimum  of 

nun 

D(f^,dm)  at  the  n"'  frequency.  Quartile  directions  satisfy 


•^/i(  ^25  55,(1  ^m. 


)  =  0.25 


4«>5055.,.- V)  =  0.50 


(45) 


(46) 


^i|(^75%,/i  ^m. 


)  =  0.75 


(47) 


The  frequency-dependent  quartile  directions  6^^^,  and  deter¬ 

mined  from  Equations  45,  46,  and  47,  respectively,  are  used  to  define  direc¬ 
tional  spread 

~  ^7556,(1  “  ^2555,(1 


and  asymmetry 


A 


n 


In 


^75%.(i  ^5055,(1 

^5055,11  “  ^2555,(1 


for  each  frequency. 


(49) 


Modal  analysis 

Modal  analysis  is  also  performed  on  D(f^,6^)  for  each  fi'equency  to  deter¬ 
mine  the  number  and  relative  intensities  of  directional  peaks.  The  reasoning 
and  method  are  the  same  for  the  frequency-dependent  analysis  as  was  dis¬ 
cussed  for  bulk  modal  analysis.  The  intent  is  to  gain  a  more  detailed  descrip¬ 
tion  of  modal  structure  as  it  may  vary  with  frequency,  and  to  get  a  more 
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quantitative  understanding  of  modal  character  based  on  individual  frequencies 
rather  than  that  obtained  from  the  somewhat  smoothed  integrated  direction 
spectrum  D{d^ . 


Summary  of  Parameters 

Frequency-direction  spectra  observed  at  Harvest  Platform  are  characterized 
using  a  set  of  bulk  parameters  to  describe  whole  spectra,  and  a  set  of  frequen¬ 
cy-dependent  parameters  to  obtain  more  detail  of  directional  energy  distribu¬ 
tions  at  individual  frequencies.  In  addition  to  numbers  and  relative  intensities 
of  directional  and  frequency  modes  for  the  spectra  as  whole  units,  nine  bulk 
parameters  f^,  0^,  0^,  a,  y,  5,  A0,  A)  yield  considerable  insight  for 
general  wind  wave  energy  distribution.  A  more  refined  examination  of  these 
spectra  is  characterization  of  directional  distributions  frequency  by  frequency. 
Seven  parameters  (0^,  0^^,  y^,  A.0^,  A^),  as  well  as  an  analysis  of 

directional  modes,  are  used  in  this  study. 

Distributions,  behavior,  and  interactions  of  the  bulk  spectral  characteristics 
are  shown  and  discussed  in  Chapter  5.  Detailed,  frequency-by-frequency 
characteristics  are  examined  in  Chapter  6.  Some  insight  into  the  accuracy  of 
the  parameters  derived  from  the  present  data  set  is  given  in  Chapter  10. 
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5  Bulk  Parametric  Behavior 


A  way  to  examine  the  1-year  set  of  directional  spectra  observed  at  Harvest 
Platform  is  through  the  behavior  of  the  characterizing  bulk  parameters.  In 
this  chapter,  a  wave  climatology  is  discussed  briefly  based  on  time  series  of 
the  bulk  parameters  illustrated  in  Appendix  A.  Statistical  distributions  of 
these  parameters  are  also  shown  and  discussed  in  this  chapter  to  quantify 
ranges,  typical  values,  and  most  common  conditions.  In  the  third  part  of  this 
chapter,  correlations  among  various  bulk  parameters  are  illustrated  to  see  if, 
for  example,  there  is  a  typical  parametric  behavior  under  high-  or  low-energy 
conditions. 


Time  Series  of  Bulk  Parameters 

To  show  general  spectral  behavior  in  a  climatological  way,  time  series 
graphs  of  characterizing  parameters  are  shown  in  Appendix  A.  Parameters 
are  plotted  in  two  groups  for  each  month  in  the  year  of  observations,  with  one 
group  on  each  of  two  facing  pages.  Left-hand  pages  illustrate  f^,  d^, 

6^,  Ad,  and  a.  Right-hand  pages  show  d,  A,  y,  6,  percent  of  wave  energy 
in  each  of  up  to  three  directional  modes,  and  percent  of  wave  energy  in  each 
of  up  to  three  frequency  modes.  As  is  shown  in  the  following  section,  three 
is  the  maximum  number  of  modes  found  in  any  of  the  S(d^)  or  S(fJ  exam¬ 
ined  in  this  report. 


Wave  height 

Characteristic  wave  height  ,  being  derived  from  total  sea  surface  dis¬ 
placement  variance,  indicates,  among  other  conditions,  the  occurrence  of  high- 
energy  events  that  are  of  concern  to  those  concerned  with  shore  protection. 
Data  from  Harvest  Platform  indicate  a  persistent,  somewhat  energetic  wave 
climate  with  heights  commonly  in  the  range  1-2  m  (3. 3-6. 6  ft).  High-energy 
events  occur  in  late  fall,  winter,  and  spring.  Numerous  events  with  wave 
heights  between  3  and  4  m  (9.8  and  13.1  ft)  occurred  in  January,  February, 
March,  and  April,  with  the  annual  cycle  beginning  again  in  November  and 
December.  The  extreme  event  of  the  year  occurred  on  2  April  with  character¬ 
istic  wave  heights  in  excess  of  5  m  (16.4  ft).  The  lowest  energy  condition 
occurred  on  1  October,  where  wave  heights  fell  to  about  0.6  m  (2.0  ft). 
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Peak  frequency 


Frequencies  representing  peak  energy  levels  in  the  observed  frequency 
spectra  cover  the  full  range  of  discrete  frequencies  used  to  define  the  spectra. 
High  peak  frequencies  occur  in  the  initial  stages  of  storm  events.  Very  low 
peak  frequencies  tend  to  occur  when  the  spectrum  is  dominated  by  waves 
from  distant  sources.  Intermediate  frequencies  occur  for  all  other  wave  condi¬ 
tions,  with  common  values  of  0.06  to  0.10  Hz.  This  behavior  is  consistent 
with  conventional  observations  of  wind  sea  and  swell  (Kinsman  1984),  and  so 
is  not  surprising. 


Characteristic  directions 

Two  parameters,  peak  direction  and  mean  direction  6^  have  been  de¬ 
fined  to  yield  single  characteristic  directions  for  observed  frequency-direction 
spectra.  As  shown  in  Appendix  A,  6^  tends  mostly  to  lie  in  the  range  50  to 
70  deg  west  of  north  (following  the  direction  convention  shown  in  Figure  2), 
indicating  that  a  common  source  of  energetic  waves  is  west-northwest  of 
Harvest  Platform.  This  region  is  the  open  North  Pacific  Ocean  to  which 
Harvest  Platform  is  well-exposed.  Variations  in  6^  during  January,  February, 
and  March  appear  to  be  associated  with  some  high-energy  events  that  are  also 
characterized  by  high  peak  frequencies,  suggesting  the  presence  of  local 
storms.  In  the  summer  months  of  July,  August,  and  September,  6^  occasion¬ 
ally  shifts  abruptly  to  directions  near  ±180  deg,  and  is  associated  with  low 
peak  frequencies.  This  behavior  suggests  dominance  of  the  wave  field  by 
southern  swell,  which  consists  of  waves  generated  in  the  South  Pacific  Ocean 
during  the  southern  winter. 

Mean  direction  6^  generally  follows  the  pattern  of  6^,  but  has  a  smaller 
magnitude  of  large  fluctuations,  and  tends  to  vary  more  smoothly  than 
The  reason  for  these  differences  is  that  6^  is  defined  as  the  peak  of  the  direc¬ 
tional  distribution  at  whereas  is  based  on  integrals  of  D{B^)  (Equa¬ 
tion  26),  which  tends  to  smooth  the  result,  and  T>(0„)  is  itself  based  on  an 
integral  of  (Equations  19  and  21),  which  provides  further  smoothing. 

As  a  consequence,  behavior  of  these  two  parameters  must  be  interpreted  care¬ 
fully.  For  example,  in  a  spectrum  with  both  local  wind  sea  and  a  dominant 
southern  swell,  6^  will  (within  the  accuracy  of  the  IMLE  method)  be  the 
southern  swell  direction.  Because  of  the  smoothing  effect  of  the  integrals 
defining  mean  direction,  6^  will  be  neither  the  southern  swell  direction  nor  the 
local  wind  sea  direction  (unless  they  are  codirectional),  but  some  direction  in 
between.  The  two  characteristic  directions  are  more  likely  to  agree  if  only 
southern  swell  is  present,  or  if  only  local  wind  sea  is  present.  There  is  also 
likely  to  be  greater  variation  in  estimates  of  6^  because  it  is  based  on  a  single 
point  in  S(f^,d^),  and  so  has  fewer  degrees  of  freedom  than  an  integrated 
parameter  like  6^ . 
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To  illustrate  how  they  are  related,  Figure  5  shows  the  correlation  of  6^  and  6^ 
based  on  the  data  shown  in  Appendix  A.  For  this  figure,  data  were  recast 
along  axes  ranging  from  0  to  270  deg  instead  of  -180  to  180  deg  by  adding 
360  deg  to  any  negative  parameter.  This  procedure  clarified  the  relationship 
among  parameters  near  ±180  deg,  and  did  not  change  the  meaning  of  any 
data.  Figure  5  shows  only  three  quadrants  because  none  of  the  characteristic 
directions  were  detected  in  the  northeast  quadrant,  which  is  a  range  of  direc¬ 
tions  with  very  short  fetches  from  the  continental  land  mass,  and  from  which 
little  wave  energy  is  expected. 


Figure  5.  Correlation  of  characteristic  direction  parameters 

Considerable  scatter  is  seen  in  Figure  5,  as  well  as  the  fact  that  the  overall 
range  of  6^  is  less  than  that  of  6^.  An  extreme  example  of  scatter  occurs  for 
0Q  near  90  deg;  6^  ranges  from  about  20  deg  to  more  than  180  deg.  The 
interpretation  of  this  behavior  is  that  numerous  directional  spectra  have  cen¬ 
troids  represented  by  in  a  small  range  near  90  deg,  but  extreme  values  that 
cover  a  broad  range  because  some  distributions  are  broad  or  multimodal,  and 

is  a  direction  well  displaced  from  the  distribution  centroid.  Some  scatter  is 
also  accounted  for  by  less  statistical  reliability  in  estimates  of  6^,  and  a  slower 
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response  of  6^  to  changing  spectral  conditions.  A  user  of  these  data  must  use 
some  caution  in  choice  of  either  of  these  characteristic  directions.  The  two 
directions  have  different  meanings,  and,  in  complex  wave  conditions,  there 
may  not  be  a  single  direction  that  represents  any  spectrum  well. 


Directional  spread  parameters 

Two  parameters  have  been  defined  in  this  report  that  characterize  angular 
spread  of  observations  of  D(6^).  Time  series  of  circular  width  a  and  quar- 
tile  spread  A6  plotted  in  Appendix  A  show  quite  similar  behaviors,  though 
their  numerical  values  differ.  Using  Afl  as  a  representative  descriptor,  direc¬ 
tional  spreads  appear  to  have  a  persistent  background  range  of  values  between 
20  and  60  deg  with  occasional  excursions  to  a  range  of  rather  large  values 
between  90  and  120  deg.  Excursion  events  appear  to  be  associated  with  some 
increases  in  wave  energy  as  shown  for  the  months  of  January,  February,  and 
March,  and  also  during  southern  swell  or  storm  events  scattered  from  late  July 
to  late  November.  At  times  (e.g.,  in  July  and  August),  there  are  increases  in 
spread  parameters  without  corresponding  shifts  in  6^  or  6^ .  This  condition 
simply  suggests  that  there  may  be  multiple  modes  of  energy,  but  the  storm  or 
southern  swell  event  does  not  dominate. 

Because  a  and  AO  appear  to  correspond  well  in  the  graphs  in  Appendix  A, 
a  scatter  diagram  has  been  constructed  to  indicate  their  correlation.  As  shown 
in  Figure  6,  the  two  parameters  seem  well  correlated  for  high  spreads  (A0 
greater  than  90  deg),  and  less  correlated  at  low  spreads  (A6  less  than  about 
60  deg).  Data  points  cluster  in  these  two  spread  ranges,  which  crudely  quan¬ 
tifies  the  dual  regimes  discussed  in  the  previous  paragraph.  Because  the 
correlation  is  reasonably  good  in  Figure  6,  it  appears  that  one  can  estimate 
reasonably  well  one  of  these  parameters  given  the  other,  especially  at  higher 
spreads. 

It  is  noted  that  Kuik,  van  Vledder,  and  Holthuijsen  (1988)  multiply  their  o 
by  180/7r  to  estimate  a  spreading  angle  in  degrees,  using  an  analogy  with 
linear  (vice  circular)  moments.  Such  an  estimate  is  more  exact  for  small 
spread,  and  appears  reasonable  for  a  =  0.5 ,  for  which  their  spread  estimate  is 
approximately  29  deg,  and  Ad  is  in  the  range  20  to  40  deg  based  on  Figure  6. 
However,  for  a  =  1 ,  their  spread  estimate  is  about  57  deg,  but  Figure  6  indi¬ 
cates  that  A0  is  in  the  range  120  to  140  deg,  which  is  quite  different.  Though 
the  numerical  values  are  different,  the  reasonable  correlation  suggested  by 
Figure  6  indicates  that  a  and  AO  detect  similar  properties  of  a  distribution. 
This  subject  should  be  pursued  further  because  a  meaningful,  universally 
understandable  measure  of  directional  spread  has  yet  to  be  defined. 
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Figure  6.  Correlation  of  directional  spread  parameters 

Depth 

Time  series  plots  of  water  depth  d  are  included  in  Appendix  A  for  general 
information.  Total  water  depth  does  not  vary  enough  to  have  a  significant 
effect  on  directional-spectral  estimates.  However,  even  though  strongly  fil¬ 
tered  by  sample  averaging  of  2-hr  16-min  records,  the  figures  in  Appendix  A 
give  an  indication  of  tide  stage,  tide  range,  and  the  occurrence  of  spring  and 
neap  tides.  Such  information  may  be  of  use  to  modelers  and  climatologists 
interested  in  this  study  site. 


Asymmetry  parameters 

There  are  two  parameters  defined  in  this  report  that  indicate  skewness  or 
asymmetry  of  a  directional  distribution  of  wave  energy.  They  are  circular 
skewness  7  and  quartile  asymmetry  A .  They  are  important  parameters  be¬ 
cause  they  indicate  the  presence  of  significant  side  lobes  of  energy  in  a  distri¬ 
bution  that  may  not  be  evident  in  mean  direction  or  directional  spread  parame¬ 
ters.  For  a  symmetric  distribution  in  the  sense  of  these  parameters,  both  A 
and  7  are  zero.  The  general  background  levels  of  these  parameters  as  shown 
in  Appendix  A  tend  not  to  be  zero,  but  are  slightly  positive  with  A  varying 
typically  firom  0  to  1,  and  7  ranging  from  0  to  2.  Large  changes  in 
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asymmetry  do  not  appear  to  be  highly  correlated  with  energetic  events,  as 
evidenced  by  the  persistent  (with  some  exceptions)  background  levels  in  the 
energetic  months  of  January  through  April. 

Asymmetry  changes  markedly  during  events  associated  with  peak  directions 
indicating  strong  waves  from  the  south,  as  seen  most  clearly  in  time  series 
plots  for  August  and  September.  Due  either  to  storms  or  southern  swell,  such 
events  impose  a  lobe  of  energy  to  one  side  of  a  pre-existing  energy  distribu¬ 
tion,  leading  to  a  strong  asymmetry  indication.  This  is  especially  true  for  A 
if  the  side  lobe  is  at  a  large  azimuthal  separation  from  the  main  lobe.  Such 
behavior  suggests  there  may  be  two  main  regimes  for  as3Tnmetry  parameters: 
one  for  conditions  where  all  waves  are  from  a  single  source  leading  to  a  nar¬ 
row  range  of  relatively  small  asymmetry  scales,  and  another  where  multiple 
wave  sources  are  involved  leading  to  another  range  of  large  asymmetry  pa¬ 
rameters. 

The  two  parameters  A  and  7  show  evidence  of  some  correlation  in  the 
time  series  plots,  but  there  are  frequent  deviations.  Figure  7  is  a  scatter  plot 
of  the  two  parameters,  and  shows  that  there  is  a  coarse  relationship  between 
the  two.  When  one  is  large,  the  other  tends  to  be  large,  but  there  is  no  clear 
relationship.  The  dark  cluster  of  points  near  7  =  1  and  A  =  0.5  illustrates  the 
persistent  background  conditions  noted  previously.  A  second  banding  of 
points  for  which  A  ranges  between  1  and  2  where  7  =  1  is  indicative  of  bimo¬ 
dality  or  secondary  lobe  conditions.  Because  the  correlation  in  general  is  not 
strong,  the  two  parameters  must  be  sensitive  to  different  aspects  of  skewed 
energy  distributions,  and  so  should  be  interpreted  separately. 


Kurtosis 

Time  series  of  circular  kurtosis  h  derived  from  D(6^),  and  shown  in 
Appendix  A,  range  from  near  1  to  about  14.  Values  in  the  upper  part  of  this 
range  appear  to  be  common.  Values  in  the  lower  part  of  this  range  occur 
when  there  is  large  Ad.  This  behavior  is  qualitatively  consistent  with  the 
usual  interpretation  of  kurtosis.  A  normal  or  Gaussian  distribution  has  5  =  3. 
Distributions  with  5  >  3  are  considered  to  be  more  peaked  than  a  Gaussian 
curve,  and  those  with  5  <  3  tend  to  be  flatter  than  a  Gaussian  curve.  The 
tendency  for  5  to  be  commonly  above  3  indicates  that  many  directional  distri¬ 
butions  have  sharper  peaks  than  a  normal  distribution.  The  occurrence  of  low 
6  with  large  Ad,  which  also  corresponds  generally  with  d^  near  ±180  deg,  is 
consistent  with  multiple  sources  of  wave  energy  from  divergent  directions. 
This  condition  leads  to  flatter  (in  the  kurtosis  sense)  distributions  than  a  nor¬ 
mal  distribution  because  there  may  be  several  widely  separated  peaks  of  about 
the  same  amplitude  in  D(0^) .  If  the  peaks  or  modes  were  analyzed  separate¬ 
ly,  each  would  likely  have  a  higher  kurtosis. 
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Figure  7.  Correlation  of  asymmetry  parameters 

Modal  analysis 

The  final  sets  of  parameters  shown  in  Appendix  A  are  time  series  of  per¬ 
centages  of  wave  energy  in  the  modes  of  the  direction  spectra  D(6^)  and  the 

normalized  frequency  spectra  Such  curves  give  an  indication  of 

the  number  and  relative  intensity  of  the  energy  modes  in  these  distributions. 
The  maximum  number  of  modes  found  in  any  of  these  distributions  was  three. 
One  or  two  modes  was  far  more  common  than  three.  As  modes  were  isolated 
in  analysis  described  in  the  previous  chapter,  they  were  ordered  by  size,  with 
the  most  energetic  mode  (the  mode  having  the  highest  fraction  of  total  energy) 
being  called  the  primary  mode,  and  subsequent  modes  being  called  secondary 
and  tertiary,  respectively. 

Among  the  directional  modal  curves  labelled  0 -modal  energy  in  Appen¬ 
dix  A,  there  were  long  periods  when  most  of  the  directional  distributions  were 
unimodal,  i.e.,  having  100  percent  of  wave  energy  in  a  single  mode.  These 
periods  occurred  in  January,  February,  March,  April,  and  toward  the  end  of 
the  year  in  December.  In  all  the  other  months,  there  were  firequent  occur¬ 
rences  of  bimodal  directional  distributions.  Among  these,  there  were  many 
with  significant  amounts  of  energy  contained  in  the  secondary  mode.  These 
conditions  indicate  multiple  sources  of  energy,  and  account  for  much  of  the 
variation  in  other  bulk  parameters  discussed  previously. 
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Numerous  multimodal  conditions  also  occur  in  frequency  spectra,  where 
the  modal  percentage  curves  are  labelled  /-modal  energy  in  Appendix  A. 
Though  occurring  throughout  the  year,  multimodal  frequency  spectra  are  most 
common  in  spring,  summer,  and  fell.  There  is  a  high  coincidence  of  multi¬ 
modal  frequency  spectra  and  multimodal  direction  spectra  in  the  months  of 
June,  July,  August,  and  September.  This  behavior  is  consistent  with  the 
existence  of  southern  swell,  where  modal  peaks  occur  not  only  at  different 
directions,  but  also  at  different  frequencies.  Multimodal  frequency  spectra  in 
spring  and  fall  are  likely  due  to  storm  events,  which  build  new  seas  on  exist¬ 
ing  swell.  In  most  of  the  occurrences  of  multimodal  frequency  spectra,  a 
significant  fraction  of  total  wave  energy  is  contained  in  the  secondary  mode. 


Bulk  Parameter  Probability  Densities 

To  gain  insight  as  to  ranges  and  probable  values  of  the  set  of  bulk  parame¬ 
ters  discussed  in  the  last  section,  bar  charts  of  probability  density  functions  are 
presented  and  discussed  in  this  section.  Being  based  on  the  whole  data  set  of 
2,339  observed  frequency-direction  spectra,  the  distributions  are  somewhat 
coarse,  without  distinction  of  seasonal  or  finer  time  scale  variations.  Nonethe¬ 
less,  they  give  a  broad  view  of  conditions  found  at  Harvest  Platform. 


Wave  height,  peak  frequency,  and  peak  direction 

Figure  8  shows  the  probability  density  of  in  terms  of  percent  of  total 
number  of  observations  in  bins  of  0.2-m  (0.66-ft)  width.  The  most  frequently 
occurring  characteristic  wave  height  is  approximately  1.6  m  (5.2  ft).  There 
are  no  observations  of  wave  heights  less  than  0.5  m  (1.6  ft),  and  only  about 
2  percent  of  the  cases  had  heights  below  the  1-m  (3.3-ft)  bin.  Extreme  char¬ 
acteristic  wave  heights  were  approximately  5  m  (16.4  ft).  The  considerable 
amount  of  wave  energy  passing  Harvest  Platform  is  seen  by  noting  that  about 
15  percent  of  observations  had  wave  heights  in  excess  of  3  m  (9.8  ft),  as 
determined  by  adding  the  percentages  in  all  bins  at  and  above  3  m  (9.8  ft). 

Figure  9  illustrates  the  distribution  of  .  The  most  common  peak  frequen¬ 
cy  is  approximately  0.07  Hz,  which  corresponds  to  peak  periods  of  about 
14  sec.  Such  periods  are  characteristic  of  swell  or  well-developed  storm 
waves.  Because  of  the  bounds  on  analyzed  fi'equencies  discussed  in  Chap¬ 
ter  2,  analyses  were  not  done  for  frequencies  less  than  0.04  Hz  or  greater  than 
0.16  Hz.  Consequently,  Figure  9  is  based  on  frequencies  between  0.04  and 
0.16  Hz,  so  that  true  peak  frequencies  outside  this  range  were  not  detected 
and  are  not  shown.  Figure  9  shows  that  very  few  cases  had  peak  frequencies 
near  O.04  Hz,  indicating  dominance  by  very  low-frequency  swell,  and  fewer 
than  1  percent  had  peak  frequencies  near  0.16  Hz,  which  typically  occur 
briefly  in  the  initial  stages  of  storms. 

The  distribution  of  6^  is  shown  in  Figure  10.  It  is  immediately  obvious 
that  the  bulk  of  energetic  peaks  originate  to  the  northwest  and  west  of  Harvest 
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Figure  1 0.  Distribution  of  peak  direction 


small  percentage  of  cases  in  this  range  of  directions  does  not  mean  that  south¬ 
ern  swell  is  present  in  only  4  percent  of  observations,  but  rather  that  it  is  the 
most  energetic  direction  in  4  percent  of  observations.  There  is  also  a  very 
small  number  of  cases  where  peak  direction  is  east-southeast.  This  is  the 
direction  of  Santa  Barbara  Channel,  which  has  enough  fetch  for  local  storms 
occasionally  to  induce  peak  waves  in  this  quadrant. 


Circular  moment  parameters 

Parameters  6^^,  a,  7,  and  6  derived  from  circular  moments  of  D{0^) 
reveal  additional  spectral  properties  when  considered  as  statistical  distribu¬ 
tions.  Figure  11  shows  the  distribution  of  As  a  characteristic  wave  direc¬ 
tion,  has  most  of  its  values  in  the  northwest  quadrant  (e.g.,  0^  in  Fig¬ 
ure  10).  However,  unlike  0^,  the  most  common  direction  is  more  to  the  west, 
and  the  influence  of  waves  from  the  south  only  moves  the  mean  direction  into 
the  west  side  of  the  southwest  quadrant.  This  behavior  occurs  because  mean 
direction  is  a  centroid  of  the  distribution,  whereas  peak  direction  is  the  ex¬ 
treme  of  the  distribution.  Equal  modes  of  energy  from  the  northwest  and 
south  then  result  in  0^,  at  some  intermediate  angle,  as  evidenced  in  Figure  11. 

The  distribution  of  circular  width  a  is  shown  in  Figure  12.  This  parame¬ 
ter  ranges  from  values  near  0.35  to  about  1.20.  There  appear  to  be  two  peaks 
in  the  distribution:  one  near  0.55,  and  another,  lesser  one  between  0.90  and 
1.00.  The  low-valued  peak  represents  narrow  directional  spreads,  and  by 
being  the  most  populous,  indicates  the  most  common  width  in  this  data  set. 
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Figure  1 1 .  Distribution  of  mean  direction 

The  high  valued  peak  represents  broad  directional  spreads,  and  is  evidently 
the  result  of  multimodal  distributions  of  energy  from  widely  divergent  sourc 
es.  Because  a  reasonable  correlation  between  a  and  the  quartile  spread  Ad 


Figure  1 2.  Distribution  of  circular  width 
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was  seen  in  Figure  6,  it  is  possible  to  assign  some  angular  measure  of  spread 
to  the  two  modal  peaks  in  Figure  12.  Specifically,  Figure  6  indicates  that  for 
a  »  0.55 ,  a  spread  angle  of  20  to  40  deg  is  common.  For  a  between  0.90 
and  1.00,  a  spread  angle  of  roughly  90  to  120  deg  may  be  expected. 

Circular  skewness  7  has  a  distribution  shown  in  Figure  13,  which  indicates 
a  range  of  -2.5  to  3.8.  The  peak  probability  of  skewness  occurs  for  7  «  0.8, 
and  most  of  the  observations  have  positive  skewness.  Only  about  10  percent 
of  the  observations  have  7  <  0 .  As  defined  in  this  report,  positive  skewness 
means  that  distributions  have  tail  regions  or  side  lobes  at  positive  angles  (rela¬ 
tive  to  that  are  not  balanced  by  similar  features  at  negative  angles.  Be¬ 
cause  resides  mostly  in  the  northwest  quadrant,  positive  skewness  suggests 
additional  wave  energy  from  the  west  or  southwest.  This  behavior  is  consis¬ 
tent  with  secondary  source  regions  for  wave  energy  at  Harvest  Platform  in  the 
central  and  south  Pacific  Ocean,  the  most  likely  alternative  source  regions  to 
those  in  the  northwest. 

The  distribution  of  bulk  circular  kurtosis  6  is  shown  in  Figure  14.  There 
is  a  broad  region  of  nearly  equal  probability  densities  for  6  in  the  range  1.5 
to  5.5,  and  a  tail  region  of  higher  values  that  extends  to  about  13.5.  As 
discussed  previously,  kurtosis  values  greater  than  3  indicate  distributions  that 
are  more  peaked  than  Gaussian  curves,  i.e.,  they  have  more  pointed  peak 
regions  and  broader  bases  than  standard  normal  distributions.  For  wave 
energy  distributions,  this  behavior  suggests  that  distributions  with  large  b  are 
dominated  by  waves  from  a  single  direction  with  little  spread  near  the  distri- 


Figure  13.  Distribution  of  circular  skewness 
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6 


Figure  14.  Distribution  of  circular  kurtosis 


bution  peaks.  Kurtosis  values  less  than  3  indicate  broad  peak  regions  with 
abrupt  drops  in  energy  at  the  outer  edges.  This  behavior  is  consistent  with 
broad  directional  distributions  of  wave  energy,  as  occurs  when  several  source 
regions  are  active  at  the  same  time.  The  sudden  drop  in  probability  density  of 
8  for  values  near  1.5  indicates  a  characteristic  behavior  that  is  seldom  exceed¬ 
ed.  In  this  set  of  data,  it  is  probably  the  coexistence  of  southern  swell  and 
local  wind  sea  from  the  northwest,  which  occurs  frequently  and  in  a  rather 
fixed  pattern. 


Quartile  parameters 

The  probability  density  of  quartile  spread  A6  is  shown  in  Figure  15. 
Overall,  the  spread  parameter  covers  a  very  broad  range  from  about  15  deg  to 
values  near  155  deg.  The  distribution  has  two  clear  peaks:  one  at  about 
30  deg,  and  a  second  at  about  105  deg.  As  suggested  by  the  circular  width 
parameter  o  in  Figure  12,  but  more  clearly  in  Figure  15,  the  narrower 
spreads  are  likely  due  to  unimodal  energy  distributions,  and  the  broader 
spreads  arise  because  of  waves  travelling  from  directionally  divergent  sources. 
Figure  15  suggests  that  unimodal  directional  distributions  have  Ad  in  the 
range  15  deg  to  approximately  70  deg,  with  higher  A6  due  to  multimodal 
distributions.  This  anticipated  behavior  is  examined  further  in  Chapter  8. 

Probability  density  of  the  quartile  measure  of  distribution  asymmetry  A  is 
illustrated  in  Figure  16.  Like  circular  skewness  y  in  Figure  13,  A  shows  a 
propensity  for  slightly  positive  values,  with  peak  probability  density  for 
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Figure  15.  Distribution  of  quartile  spread 


A  ^  0.2,  and  about  75  percent  of  cases  with  A  greater  than  zero.  This  be¬ 
havior  supports  the  interpretation  of  Figure  13,  which  was  that  positive  asym¬ 
metry  occurs  because  of  dominant  waves  from  the  northwest  and  simultaneous 
secondary  waves  from  the  west  and  south.  As  with  Ad ,  the  verity  of  this 
interpretation  is  shown  in  Chapter  8. 
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Figure  16,  Distribution  of  quartile  asymmetry 
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Modal  analysis 


Distributions  of  modal  parameters  based  on  D(0^)  and  normalized  fre¬ 
quency  spectra  are  shown  in  Figures  17  and  18,  respectively.  In  Figure  17, 
the  probability  density  of  number  of  directional  modes  is  shown  on  the  left- 
hand  side.  A  maximum  of  3  modes  was  found  in  any  of  the  Di_6^)  analyzed, 
and  that  condition  occurred  in  only  about  1  percent  of  the  cases.  The  most 
common  number  of  modes  is  2,  which  occurred  in  approximately  60  percent 
of  the  observations.  The  remaining  cases,  about  39  percent,  were  unimodal  in 
direction.  Because  of  the  highly  constrained  definition  of  modal  separation  as 
described  in  Chapter  4,  the  bimodal  distributions  discerned  in  this  study  are 
quite  distinct,  and  suggest  clear,  directionally  deviate  sources  of  wave  energy. 

The  right-hand  side  of  Figure  17  shows  the  probability  density  of  percent 
of  wave  energy  in  the  primary,  or  most  energetic,  mode.  This  computation 
was  done  because  there  could  be  a  strong  indication  of  bimodal  spectra,  but 
the  secondary  modes  could  all  have  very  low  energy,  and  consequently  be  of 
little  interest.  Figure  17  indicates  that  approximately  45  percent  of  all  cases 
had  between  97.5  and  100  percent  in  the  primary  mode.  This  bin  has  more 
cases  (45  percent)  than  the  number  of  unimodal  cases  (39  percent)  because 
some  secondary  modes  (roughly,  6  percent  of  cases)  had  less  than  2.5  percent 
of  the  energy  in  a  distribution,  and  so  were  included  in  the  95-to-lOO-percent 
bin.  Slightly  less  than  10  percent  of  all  cases  had  92.5  to  97.5  percent  of  total 
energy  in  the  primary  mode,  indicating  2.5  to  7.5  percent  of  energy  was  in  a 
secondary  mode.  By  adding  adjacent  bin  percentages,  it  can  be  seen  that 
approximately  25  percent  of  observations  had  secondary  modes  with  17.5  to 
50  percent  of  total  energy.  This  result  means  the  amount  of  energy  in  second¬ 
ary  modes  is  often  of  considerable  in:q)ortance. 


Modes  Percent  0-Modal  Energy 


Figure  17.  Distribution  of  directional  modes 


Chapter  5  Bulk  Parametric  Behavior 


o 


Figure  1 8.  Distribution  of  frequency  modes 

Figure  18  illustrates  results  of  modal  analysis  for  normalized  frequency 

spectra  165(/„)///^.  A  maximum  of  three  modes  was  found,  with  approxi¬ 
mately  1  percent  of  cases  having  three  modes,  37  percent  of  cases  having  two 
modes,  and  62  percent  of  cases  being  unimodal.  Note  that  the  percentage  of 
unimodal  frequency  spectra  need  not  be  the  same  as  the  percentage  of  unimo¬ 
dal  direction  spectra  because  two  directional  modes  can  occur  at  the  same  fre¬ 
quency,  or  two  frequency  modes  can  occur  in  the  same  direction.  Figure  18 
shows  that  approximately  62  percent  of  all  cases  have  between  97.5  and  100 
percent  of  wave  energy  in  the  primary  mode.  About  25  percent  of  all  cases 
have  between  17.5  and  47.5  percent  of  wave  energy  in  secondary  and  tertiary 
modes,  again  making  secondary  modes  important.  It  is  likely  that  some  of  the 
secondary  frequency  modes  are  the  same  as  the  secondary  direction  modes  in 
cases  where  southern  swell  propagates  through  a  local  wind  sea. 


Parameter  Correlations  with  Wave  Height 

Of  considerable  interest  is  characterization  of  the  wave  field  at  Harvest 
Platform  in  high-energy  conditions.  High-energy  waves,  as  indicated  by  large 
,  passing  this  site  from  the  open  Pacific  Ocean  are  directed  toward  the 
California  coast,  and  it  is  typically  for  these  waves  that  a  coastal  engineer 
must  design  shore  protection  structures  or  beach  renourishment  projects.  A 
reasonable  description  of  wave  fields  in  high-energy  conditions  can  be 
achieved  by  correlating  bulk  spectral  parameters  with  to  determine  ranges 
of  these  parameters  that  apply  in  extreme  conditions. 
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Peak  frequency  and  peak  direction 

Figure  19  relates  and  .  There  are  a  total  of  13  smooth  frequency 
bands  in  all  frequency  related  spectra  discussed  in  this  report.  Figure  19 
shows  that  all  bands  contain  spectral  peaks  for  some  observations.  At  low 
energy,  for  which  «  1  m  (3.3  ft),  ranges  fi’om  0.05  to  0.16  Hz.  At 
intermediate  energy,  with  =  3  m  (9.8  ft),  peak  frequency  appears  in  the 
lowest  frequency  band  (0.04  Hz),  but  only  ranges  to  about  0.12  Hz,  probably 
because  it  is  normal  for  peak  frequency  to  move  to  lower  values  in  growing 
wind-driven  seas.  At  the  highest  energies  >  4  m  (13.1  ft)],  peak  fre¬ 
quency  typically  resides  in  a  narrow  range  from  0.05  to  0.08  Hz. 

Figure  20  shows  the  correlation  of  6^  with  .  Clearly,  the  bulk  of  cases 
have  peak  directions  in  the  northwest  to  west,  and  this  property  holds  even  for 
the  highest  waves.  Dominance  of  the  spectrum  by  southern  energy  sources 
occurs  only  for  lower  energy  conditions.  Rare  events  provide  intermediate 
energy  levels  with  peak  directions  from  the  southwest  or  east-southeast. 


Circular  moment  parameters 

Correlation  of  H  with  mean  direction  derived  from  circular  moment 

mo  u 

n,  is  shown  in  Figure  21.  The  gross  pattern  of  points  in  Figure  21  is  like  the 


Figure  19.  Correlation  of  wave  height  with  peak  frequency 
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Figure  21.  Correlation  of  wave  height  and  mean  direction 
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pattern  for  peak  direction  shown  in  Figure  20.  The  bulk  of  points  in  Figure 
21  are  in  the  northwest  to  west  directions  at  intermediate  to  high  .  For 
low  wave  heights,  there  is  a  deviation  in  mean  direction  to  the  southwest. 

This  behavior  is  likely  due  to  the  averaging  effect  on  wind  sea  from  the  north¬ 
west  quadrant  and  southern  wave  sources  in  the  computation  of  0^;  for  two 
distinct  modes  of  energy,  mean  direction  will  fall  between  the  modal  peak 
directions.  Of  note  is  the  fact  that  mean  and  peak  directions  are  approximate¬ 
ly  equal  (to  within  ±5  deg)  in  high-energy  conditions  for  which  >  4  m 
(13.1  ft).  This  behavior  suggests  that  directional  distributions  are  naturally 
narrow  and  symmetric  in  these  conditions. 

The  suggested  narrowness  of  directional  distributions  at  high  energy  is 
verified  in  Figure  22,  which  is  a  correlation  of  with  circular  width  a . 

The  broad  swath  of  points  indicates  generally  large  a  for  small  ,  and 
small  a  for  large  For  high  energy  >  4  m  (13.1  ft)],  a  is  in  the 
range  0.3  to  0.5,  which  clearly  includes  the  smallest  observed  circular  widths. 
Based  on  the  correlation  shown  in  Figure  6,  this  range  of  a  occurs  for  quar- 
tile  spreads  A0  in  the  range  10  to  40  deg,  which  are  rather  narrow. 

Evidently,  the  small  directional  spreads  at  high  energy  are  indicative  of 
moderately  low  values  of  circular  skewness  7,  as  shown  in  Figure  23. 

Though  there  is  a  large  scatter  of  points  in  Figure  23,  the  grouping  tightens 


Figure  22.  Correlation  of  wave  height  with  circular  width 
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Figure  23.  Correlation  of  wave  height  with  circular  skewness 


somewhat  for  high  .  For  >  4  m  (13. 1  ft),  7  is  in  the  range  -0.2  to 
2.2,  which  is  considerably  narrower  than  the  range  of  -2.5  to  3.8  indicated  for 
near  1  m  (3.3  ft).  Furthermore,  because  the  high-energy  spectra  have  small 
directional  spreads,  side  lobes  of  energy  that  cause  a  non-zero  7  may  not  be 
greatly  removed  from  the  main  energy  peaks.  Consequently,  high-energy 
spectra  may  appear  relatively  unskewed.  These  properties  will  be  examined 
farther  in  Chapter  8  where  characteristic  directional  distributions  are  dis¬ 
cussed. 

Figure  24  shows  how  circular  kurtosis  b  is  related  to  .  The  main 
cluster  of  points  indicates  a  gross  tendency  for  5  to  increase  with  increasing 
being  less  than  4  for  the  smallest  wave  heights,  but  ranging  approximate¬ 
ly  between  4  and  9  for  >  4  m  (13.1  ft).  Kurtosis  greater  than  3  indicates 
that  maxima  of  the  directional  distributions  are  more  peaked  than  a  Gaussian 
curve.  This  property  of  high-energy  directional  distributions  seems  consistent 
with  narrowly  spread  curves  of  small  to  moderate  skewness. 


Quartile  spread  and  asymmetry 

Figures  25  and  26  illustrate  the  correlation  of  with  quartile  spread  M 
and  quartile  asymmetry  A,  respectively,  as  derived  from  D(d^).  Quartile 
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Figure  24.  Correlation  of  wave  height  with  circular  kurtosis 
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spread  spans  a  very  large  range  for  low-energy  distributions,  but  narrows  in 
range  and  becomes  moderately  small  for  >4m(13.1ft).  There  appear 
to  be  two  clusters  of  A0  for  between  1  and  2  m  (3.3  and  6.6  ft).  The 
largest  group,  having  a  A(9  range  of  20  to  approximately  70  deg,  is  associated 
with  wind  seas  and  swell  from  a  single  direction.  The  smaller  group,  with 
A0  >  90  deg,  is  attributed  to  wave  trains  from  multiple  source  directions. 

For  high-energy  spectra,  A0  is  confined  to  the  range  15  to  45  deg.  For  the 
extreme  high  ,  A$  =  30  deg,  which  happens  to  be  the  most  common 
spread  as  shown  in  Figure  15. 

Quartile  asymmetry  also  shows  a  wide  range  for  low  and  a  narrow 
range  for  high  ,  as  shown  in  Figure  26.  The  reasons  for  this  behavior  are 
the  same  as  for  the  behavior  of  A6 ;  multiple  wave  sources  with  varying  inten¬ 
sities  induce  large  apparent  asymmetries  in  the  total  directional  distributions, 
whereas  intense  single  sources  lead  to  more  nearly  symmetric  distributions. 
For  the  extreme  high  in  Figure  26,  A  is  nearly  zero  (i.e.,  symmetric  in 
the  sense  of  the  definition  of  A).  For  >  4  m  (13.1  ft),  A  is  roughly  in 

the  range  -0.6  to  0.6,  and  so  is  nearly  symmetric. 


Figure  26.  Correlation  of  wave  height  with  quartile  asymmetry 
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Summary  of  Bulk  Parameters 


In  this  chapter,  parameters  characterizing  integrated  direction  spectra 

D{d^  and  normalized  frequency  spectra  16S(/„)//f^  have  been  examined  in 
several  ways  to  deduce  some  of  the  directional  behavior  of  waves  at  Harvest 
Platform.  Time  series  graphs  in  Appendix  A  show  timing,  duration,  and 
seasonal  variation  of  these  parameters  as  a  gross  characterization  of  the  wave 
climate.  Statistical  distributions  of  individual  parameters  in  probability  density 
form  indicate  ranges  and  most  common  values  of  these  parameters,  and  sug¬ 
gest  ways  that  groups  of  parameters  are  affected  in  various  wave  regimes. 
Correlations  of  bulk  parameters  with  characteristic  wave  height  reveal  ranges 
and  groupings  of  parameters  at  various  energy  levels,  especially  in  conditions 
of  high  energy. 

Table  1  contains  a  summary  of  some  of  the  indicator  statistics  described  in 
this  chapter.  A  distinction  is  made  between  all  observations  and  a  subset 
called  high  energy,  which  has  been  defined  arbitrarily  as  any  condition  in 
which  the  characteristic  wave  height  exceeds  4  m  (13.1  ft).  Minimum  and 
maximum  values  are  listed  for  each  of  the  nine  bulk  parameters  discussed  in 
this  chapter  for  each  of  the  two  data  groupings.  The  column  labeled  mode  in 
the  grouping  for  all  data  is  the  most  common  value  of  each  parameter  as 
deduced  from  the  maximum  probability  density  shown  in  the  set  of  figures  in 
the  second  part  of  this  chapter.  Bulk  parameters  discussed  in  this  chapter  are 
gross  indicators  of  frequency-direction  spectral  structure,  but  may  be  indica¬ 
tive  of  greater  detail.  Details  of  spectral  structure  are  examined  more  closely 
in  subsequent  chapters,  using  bulk  parameters  for  general  categorization. 


Table  1 

Summary  of  Bulk  Parameter  Statistics 

Parameter 

All  Data 

High  Energy  >  4  m)  | 

Minimum 

Mode 

Maximum 

Minimum 

Maximum 

(m) 

0.61 

1.60 

5.08 

4.00 

5.08 

f,  (Hz) 

0.04 

0.07 

0.16 

0.05 

0.08 

Bp  (deg) 

-180 

50 

180 

45 

100 

Oo  (deg) 

-180 

70 

177 

49 

105 

a 

0.33 

0.55 

1.21 

0.38 

0.58 

Y 

-2.35 

0.80 

3.76 

-0.10 

2.15 

6 

1.22 

4.00 

13.68 

4.90 

8.70 

M  (deg) 

11 

30 

155 

13 

43 

A 

-1.77 

0.20 

2.59 

-0.67 

0.49 
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6  Frequency-Dependent 
Parameters 


More  detail  about  frequency-direction  spectral  structure  is  obtained  by 
examining  parameters  derived  from  directional  distribution  functions 
on  a  frequency-by-frequency  basis.  At  this  level,  properties  of  directional 
distributions  D{d^)  based  on  integrated  frequency-direction  spectra  can  be 
explained  in  terms  of  contributions  from  the  various  frequency  bands  that  are 
the  basis  of  integration.  For  example,  large  directional  spreads  or  asymme¬ 
tries  may  be  due  to  influences  from  a  limited  range  of  frequencies,  which,  in 
turn,  may  be  of  local  or  more  general  importance  to  students  of  wave  genera¬ 
tion  and  propagation  if  the  influential  frequencies  are  low  or  high. 

In  this  chapter,  statistical  distributions  of  nine  frequency-dependent  proper¬ 
ties  are  examined.  Spectral  density  S(f^)  and  peak  direction  6^  are  two  of 
them,  and  indicate  energy  levels  and  dominant  direction  of  wave  propagation. 
The  subscript  n  remains  an  index  of  discrete  frequency,  and  is  useful  for 
distinguishing  frequency-dependent  properties  from  bulk  parameters.  The 
four  circular  moment  parameters  defined  by  Kuik,  van  Vledder,  and 
Holthuijsen  (1988)  are  included  to  indicate  moment  properties  of  distributions. 
These  four  parameters  are  mean  direction  ,  circular  width  ,  circular 
skewness  ,  and  circular  kurtosis  8^ .  Use  of  these  parameters  at  the  fre¬ 
quency  level  is  more  appropriate  to  the  intent  of  the  paper  by  Kuik,  van 
Vledder,  and  Holthuijsen  (1988),  though  these  parameters  in  bulk  form  were 
quite  useful  in  analyzing  D{6^)  in  the  previous  chapter  of  this  report.  Quar- 
tile  spread  Ad^  and  quartile  asymmetry  are  used  to  augment  interpretation 
of  circular  width  and  skewness,  respectively.  Finally,  a  modal  analysis  is 
done  for  the  set  of  directional  distributions  at  each  frequency. 


Spectral  Density  and  Peak  Direction 

Figure  27  shows  statistical  distributions  of  frequency  spectral  densities  at 
each  frequency  in  the  form  of  probability  density  functions.  The  intent  is  to 
determine  which  frequencies  have  the  highest  wave  energy,  ranges  of  wave 
energy  observed  at  each  frequency,  and  most  probable  wave  energies  at  each 
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frequency.  Note  that  the  abscissas  in  Figure  27  are  logarithmic  scales  so  that 
each  numeral  represents  a  power  of  10.  Because  resolution  bandwidths  are 
constant  at  approximately  0.01  Hz,  sea  surface  displacement  variances  repre¬ 
sented  on  the  abscissas  include  values  from  10^  to  1  m^  (lO"^  ^  to  10“”  ft^),  a 
rather  large  range.  Note  also  that  the  order  of  the  subplots  is  such  that  fre¬ 
quency  increases  in  raster  fashion,  across  the  page,  and  then  down. 

There  is  a  clear  distinction  of  spectral  content  among  groups  of  frequencies 
shown  in  Figure  27.  At  the  lowest  frequency,  energy  density  tends  to  be 
lower  than  at  any  of  the  other  frequencies,  with  a  most  common  variance 
density  of  10"^  m^/Hz  (lO^’^  ft^/Hz),  and  a  maximum  near  10  m^/Hz 
(102.03  jfj2/jj2).  For  the  next  three  frequencies,  /=  0.054,  0.064,  and 
0.074  Hz,  the  modal  variance  density  increases  in  orderly  fashion  to  about 
10^  *  mVHz  (10®  “  ft^/Hz),  and  then  to  1  mVHz  (10‘  ft^/Hz).  These  three 
frequencies  are  distinguished  by  having  variance  densities  up  to  10'  ®  mVHz 
(10^®  ft^/Hz),  higher  than  any  other  frequencies.  This  behavior  is  consistent 
with  spectral  peak  frequencies  equalling  one  of  these  frequencies  in  the  high- 
energy  conditions  shown  in  Figure  19.  At /=  0.084  Hz,  variance  density  is 
nearly  uniformly  distributed  between  10  '  and  10  mVHz  (10®  ®^  and 
102.03  ft2/jj2).  For  />  0.093  Hz,  the  distributions  are  rather  uniform  in 
shape,  ranging  roughly  from  10  '  to  10  m^/Hz  (10®  °^  to  10^  ®^  ft^/Hz),  with 
modal  values  from  1  to  KIP'*  m^/Hz  (10'  to  10'  ft^/Hz).  The  similarity  in 
the  distributions  of  the  high-frequency  wave  energy  is  interesting  because  it 
means  that  each  of  these  frequencies  has  a  given  energy  level  roughly  the 
same  fraction  of  the  time. 

Distributions  of  frequency-dependent  peak  direction  are  shown  in  Fig¬ 
ure  28.  It  is  readily  seen  that  all  frequencies  have  signifrcant  contributions  of 

from  the  northwest  quadrant,  with  such  conditions  being  dominant  for 
/  >  0.074  Hz.  Wave  fields  with  peak  directions  in  the  southwest  and  south 
have  comparable  or  dominant  percentages  for  the  three  lowest  frequencies. 

As  illustrated  in  Figure  20,  distributions  with  bulk  peak  directions  in  the 
southwest  and  south  are  of  primarily  low  and  occasionally  intermediate  energy 
(as  represented  by  H^).  The  same  range  of  energy  is  represented  in  Figure 
20  by  peak  waves  from  the  west  to  southwest,  which  appear  in  Figure  28  to 
have  larger,  though  not  dominant,  percentages  in  the  upper  half  of  the  overall 
frequency  band.  Though  difficult  to  see  in  Figure  28,  the  intermediate-energy 
waves  from  the  east-southeast  shown  in  Figure  20  appear  to  be  represented  in 
distributions  of  6^  at  the  three  highest  frequencies  in  Figure  28. 

Note  that  parameters  shown  in  Figures  20  and  28  differ  in  that  Figure  20 
shows  bulk  peak  direction  based  on  an  integration  of  the  frequency-direction 
spectra  over  all  frequencies.  Figure  28  shows  distributions  of  directional 
peaks  at  individual  frequencies.  Because  the  two  figures  agree  qualitatively  in 
content,  it  can  be  deduced  that  waves  at  frequencies  away  from  tend  to 
align  with  6^  (i.e.,  6^  «  0p),  in  a  given  frequency-direction  spectrum.  This 
behavior  is  consistent  with  the  idea  that  waves  at  different  frequencies  from  a 
given  source  all  tend  to  have  the  same  characteristic  direction. 
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Figure  28.  Distributions  of  peak  direction  for  all  frequencies 
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Circular  Moment  Parameters 


Distributions  of  frequency-dependent  mean  direction  6^^  are  shown  in 
Figure  29,  which  indicates  that  mean  directions  have  distributions  that  are 
qualitatively  similar  to  peak  directions  (Figure  28)  at  corresponding  frequen¬ 
cies.  However,  it  is  noted  that  distributions  of  mean  direction  are  smoother 
and  have  fewer  contributors  at  extreme  directions  than  those  of  peak  direc¬ 
tions.  This  relative  behavior  of  frequency-dependent  peak  and  mean  direc¬ 
tions  exists  for  the  same  reasons  given  for  the  relation  between  bulk  peak  and 
mean  directions  discussed  in  Chapter  5;  where  directional  distributions  are 
broad  or  have  multiple,  directionally  divergent  modes,  peak  direction  will  be 
the  maximum  of  the  dominant  mode,  whereas  mean  direction  will  be  some 
direction  intermediate  between  the  modes  because  of  the  smoothing  nature  of 
its  integral  definition  (Equation  40). 

Consequently,  when  there  are  directional  modes  from  both  southern  swell 
and  western  storms,  peak  direction  will  indicate  either  west  or  south,  depend¬ 
ing  on  the  peak  of  the  dominant  mode,  and  mean  direction  will  indicate  some¬ 
where  in  the  southwest,  depending  on  the  areas  under  the  two  modes.  This 
behavior  is  quite  evident  at  the  lowest  three  frequencies  in  Figures  28  and  29; 
peak  direction  has  a  significant  fraction  of  cases  from  due  south,  and  mean 
direction  has  very  few.  The  behavior  of  and  distinction  between  these  two 
parameters  is  important  in  modeling  efforts.  Neither  represents  bimodal 
spectra  very  well  (analysis  of  individual  modes  is  required  for  this).  Models 
of  coastally  propagating  waves,  which  are  extremely  direction-sensitive 
(O’Reilly  and  Guza  1991,  1993),  will  yield  different  results  for  unimodal 
directional  distributions  if  one  or  the  other  of  these  parameters  is  used  as 
characteristic  direction,  and  can  yield  quite  incorrect  results  in  bimodal  situa¬ 
tions  if  either  parameter  is  used. 

Figure  30  shows  distributions  of  direction-dependent  circular  width  . 
Extreme  widths  are  seen  at  the  lowest  frequency,  which  is  likely  due  to  a 
combination  of  waves  from  widely  divergent  sources  and  analysis  noise. 
Analysis  noise  is  more  noticeable  at  the  lowest  frequency  because  wave  ener¬ 
gy  levels  are  typically  small,  as  shown  in  Figure  27.  For  most  frequencies, 
the  most  common  a  is  near  0.5.  There  are  larger  percentages  of  high-width 
cases  for  /  in  the  range  0.064  to  0.093  Hz.  These  frequencies  are  common 
to  both  wind  seas  and  southern  swell,  and,  further,  are  often  energetic  so  that 
large  circular  widths  are  likely  due  to  multiple-source  waves.  The  narrowest 
distributions,  with  a  «  0.3,  occur  at  low  and  intermediate  frequencies  and  are 
likely  due  to  unimodal  swell  and  well-developed,  isolated  seas  created  by 
steady  winds. 

Distributions  of  frequency-dependent  circular  skewness  are  shown  in 
Figure  3 1 .  The  lowest-frequency  energy  distributions  have  a  narrow  range  of 
skewness  with  dominant  values  zero  or  slightly  negative,  suggesting  a  balance 
of  wave  energy  distributions  about  the  mean  directions  for  this  frequency. 

The  widest  ranges  of  skewness  occur  at  intermediate  and  lower  frequencies. 
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Figure  29.  Distributions  of  mean  direction  for  all  frequencies 
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The  modal,  or  most  common,  values  shift  from  slightly  negative  at  lower 
frequencies  to  slightly  positive  at  higher  frequencies.  Negative  skewness 
indicates  side  lobes  at  positive  angles  from  the  mean  direction.  In  the  coordi¬ 
nate  system  used  in  this  report,  negative  skewness  is  consistent  with  mixed 
western  wind  seas  and  southern  swell,  with  the  wind  sea  energy  dominating. 
When  southern  swell  energy  dominates,  skewness  becomes  positive.  Fig¬ 
ure  31  shows  this  latter  case  to  occur  often  at  the  lower  frequencies.  At  the 
highest  frequencies,  ranges  of  skewness  are  small,  with  most  common  values 
zero  or  slightly  positive,  suggesting  nearly  symmetric  energy  distributions. 

Circular  kurtosis  distributions  are  shown  in  Figure  32.  As  discussed  in 
Chapter  5  for  bulk  parameters,  kurtosis  is  a  measure  of  flatness  or  peakedness 
of  a  distribution,  with  6  =  3,  the  kurtosis  of  a  Gaussian  distribution,  being  a 
standard  comparator.  In  Figure  32,  the  distribution  of  8  at  the  lowest  fre¬ 
quency  is  predominantly  at  values  less  than  3,  indicating  that  most  of  the 
distributions  are  relatively  flat.  For  the  highest  seven  frequencies,  the  distri¬ 
butions  of  kurtosis  range  from  about  1  to  as  high  as  20,  with  modal  values 
typically  between  4  and  6.  These  high-frequency  distributions  indicate  that 
typical  directional  energy  distributions  are  more  peaked  than  Gaussian  curves. 
For  the  lower  frequencies  (except  the  lowest),  the  distributions  are  rather 
uniform  for  8  in  the  range  2  to  6,  with  a  significant  fraction  of  cases  having 
8  >  6.  These  low-frequency  results  indicate  a  mix  of  distributions  with  some 
very  peaked  curves,  suggesting  well-directed  waves  with  single  direction 
modes,  and  some  very  flat  curves,  indicating  waves  from  directionally  diverse 
sources. 


Quartile  Spread  and  Asymmetry 


Directional  spread  Ad^  and  asymmetry  parameters  based  on  quartile 
directions  of  energy  distributions  at  each  frequency  can  be  compared  to 
frequency-dependent  circular  width  and  skewness  ,  respectively,  to  aug¬ 
ment  interpretation.  Figure  33  shows  distributions  of  directional  spread  de¬ 
duced  from  the  set  of  observations.  At  the  lowest  frequency,  there  is  a  broad 
range  of  spread  parameter,  like  the  circular  width  at  this  frequency  in  Figure 
30,  but  there  is  a  most  common  value  of  about  25  deg,  suggesting  the  fre¬ 
quent  occurrence  of  very  narrow  directional  distributions.  For  frequencies 
higher  than  0.054  Hz,  there  appear  to  be  two  modes  in  Ad  distributions. 

With  a  boundary  between  these  two  modes  at  roughly  60  to  80  deg,  direction¬ 
al  spreads  in  the  high-spread  mode  can  be  as  high  as  100  to  140  deg,  which 
suggests  the  presence  of  multiple  wave  sources.  The  low-spread  modes  in 
Figure -33  indicate  Ad  generally  in  the  range  10  to  60  deg,  which  are  typical 
of  unimodal  seas  from  single  sources  in  various  stages  of  development  (Long 
and  Oltman-Shay  1991). 

Figure  34  shows  distributions  of  quartile  asymmetry.  The  shapes  of  sever¬ 
al  of  the  distributions  in  Figure  34  look  very  much  like  the  distributions  of 
circular  skewness  shown  in  Figure  31.  This  observation  applies  to  the  distri- 
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buttons  at /=  0.044,  0.084,  0.093,  and  0.103  Hz,  which  suggests  that  7  and 
A  are  measures  of  the  same  property.  At  other  frequencies,  there  are  subtle 
differences  in  the  distributions  shown  in  Figures  31  and  34  that  suggest  the 
two  parameters  are  more  remotely  related.  For  example,  at  frequencies  0.054 
and  0.064  Hz,  the  modes  in  distributions  of  7  are  slightly  negative,  whereas 
modes  in  distributions  of  A  are  slightly  positive.  For  these  two  frequencies, 

7  indicates  common  skewness  with  mild  secondary  lobes  at  directions  smaller 
than  the  mean  directions,  whereas  A,  if  interpreted  as  a  skewness,  indicates 
the  opposite.  Additionally,  at  the  highest  four  frequencies,  A  indicates  higher 
percentages  of  negative  values  than  does  7.  However,  for  all  />  0.074  Hz, 
the  most  common  values  of  both  parameters  are  slightly  positive.  The  general 
indication  is  that  the  two  parameters  are  closely  related  so  that  either  can 
indicate  distribution  asymmetry. 


Modal  Analysis 

Within  a  frequency-direction  spectrum,  the  directional  distribution  at  any 
given  frequency  can  have  several  modes,  or  peaks,  depending  on  the  direc¬ 
tional  separation  and  strength  of  sources  of  waves  at  that  frequency,  as  well  as 
low-level  peaks  that  are  due  to  fluctuations  in  background  noise  levels.  Modal 
analysis  discussed  in  Chapter  5  was  based  on  the  normalized,  integrated  direc¬ 
tion  spectrum  D(6^).  By  virtue  of  its  integration  over  all  frequencies,  such  a 
directional  distribution  is  smoother  (i.e.,  has  fewer  modes)  in  general  than  the 
frequency-dependent  directional  distributions  that  contributed  to  the  integral. 
Consequently,  it  is  useful  to  examine  the  statistical  distributions  of  numbers 
and  relative  importance  of  modes  at  each  frequency. 

Figure  35  is  the  result  of  such  an  analysis.  For  each  frequency,  there  is  a 
probability  density  function  for  the  number  of  modes,  and  another  probability 
density  for  the  percentage  of  wave  energy  in  the  primary,  or  largest,  mode. 

As  many  as  five  modes  were  found  in  some  frequency-dependent  distributions, 
whereas  a  maximum  of  three  modes  were  found  in  the  bulk  parameterization 
shown  in  Figure  17.  The  lowest  frequency  shown  in  Figure  35  is  the  easiest 
to  evaluate.  Approximately  96  percent  of  all  directional  distributions  were 
unimodal,  and  approximately  96  percent  of  all  distributions  had  97.5  to 
100  percent  of  energy  in  the  primary  mode. 

For  the  following  three  low  frequencies  (f=  0.054,  0.064,  and  0.074  Hz), 
there  were  four  or  five  modes  represented  by  significant  fractions  of  observa¬ 
tions.  At  these  frequencies,  fewer  than  40  percent  of  observations  had  97.5  to 
100  percent  of  energy  in  the  primary  modes,  and  fewer  than  60  percent  of 
observations  had  92.5  to  100  percent  of  energy  in  the  primary  modes.  This 
result  means  a  significant  fraction  of  cases  had  7.5  to  50  percent  of  energy 
distributed  among  secondary  modes.  This  is  important  because  two  of  these 
frequencies  are  peak  frequencies  in  conditions  of  highest  energy  (Figure  19), 
and  are  also  common  frequencies  for  southern  swell  (Figure  28). 
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For  /  >  0.093  Hz,  most  of  the  distributions  are  described  by  one  or  two 
modes,  with  bimodal  distributions  being  most  common.  Between  50  and 
74  percent  of  observations  had  97.5  to  100  percent  of  wave  energy  in  the 
primary  mode,  and  73  to  88  percent  of  observations  had  87.5  to  100  percent 
of  wave  energy  in  the  primary  mode.  Thus,  12  to  27  percent  of  cases  had 
12.5  to  50  percent  of  energy  in  secondary  modes.  This  indicates  that  although 
most  of  the  high-frequency  observations  were  bimodal,  not  a  lot  of  the  sec¬ 
ondary  modes  were  energetic. 
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This  chapter  contains  a  detailed  examination  of  statistical  distributions  of 
spectral  properties  and  parameters  associated  with  individual  frequency  bands. 
These  properties  include  spectral  density,  peak  direction,  mean  direction, 
higher  circular  moments  (width,  skewness,  and  kurtosis),  quartile  directional 
spread  and  asymmetry,  and  analysis  of  directional  modes.  As  illustrated  in 
Figures  27  to  35,  there  are  generally  significant  differences  in  these  properties 
over  the  full  range  of  frequencies  examined.  Furthermore,  parametric  behav¬ 
ior  at  the  frequency  level  does  not  necessarily  parallel  bulk  parametric  distri¬ 
butions.  Knowledge  of  parametric  behavior  at  the  frequency  level  enables 
classification  of  directional  distribution  functions  for  further  display  and  analy¬ 
sis,  which  is  done  in  subsequent  chapters  of  this  report. 
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7  Examples  of  Spectra 


Part  of  the  purpose  of  this  report  is  to  characterize  frequency-direction 
spectral  observations  from  Harvest  Platform  by  illustration.  This  is  a  some¬ 
what  daunting  task  because  there  are  more  than  2,000  observations  in  the  data 
set,  and,  as  indicated  by  the  bulk  and  frequency-dependent  parameters  dis¬ 
cussed  in  Chapters  5  and  6,  there  is  a  very  broad  variety  in  energy  distribu¬ 
tions.  Space  does  not  allow  picturing  all  the  spectra,  so  a  subset  of  represen¬ 
tative  spectra  was  identified  for  illustration. 

To  define  a  representative  subset,  each  member  of  the  whole  set  of  obser¬ 
vations  was  classified  in  terms  of  its  quartile  spread  Ad  and  asymmetry  A . 
These  two  bulk  parameters  give  indications  of  the  directional  arcs  over  which 
wave  energy  is  distributed,  and  the  existence  and  relative  importance  of  side 
lobes  and  modes  in  the  distributions.  Observations  were  grouped  in  classes  of 
Ad  that  were  10  deg  wide,  and  centered  at  10-deg  intervals,  starting  at 
10  deg.  The  narrowest  spread  class  thus  included  observations  with  Ad  in  the 
range  5  to  15  deg.  Sixteen  spread  classes  were  employed,  so  the  broadest 
spread  class  satisfied  155  deg  <  Ad  <  165  deg.  Within  each  spread  class, 
observations  were  further  grouped  by  /I  in  1 1  bins  of  width  0.4,  starting  at 
A  =  -l.S,  and  extending  to  y4  =  2.6 .  The  most  symmetric  distributions  were 
thus  in  the  asymmetry  class  for  which  -0.2  <  A  <  0.2. 

By  counting  the  number  of  observations  in  each  of  the  classification  cells, 
the  most  common  conditions  (within  the  definition  used  here)  could  be  identi¬ 
fied.  Example  spectra  were  chosen  from  24  of  the  most  populous  cells. 
Because  there  is  much  interest  in  high-energy  conditions,  24  examples  were 
chosen  based  on  the  case  with  the  highest  in  each  of  the  populous  cells. 

In  three  of  the  populous  cells  with  Ad  near  100  deg,  additional  cases  with  low 
or  intermediate  energy  were  included  for  intercomparison. 

A  total  of  28  examples  are  shown  in  this  report.  One  of  them  is  Figure  3, 
discussed  in  conjunction  with  primary  data  analysis  in  Chapter  3.  The  re¬ 
mainder  are  shown  in  Appendix  B.  Table  2  is  a  listing  of  the  sample  cases, 
and  includes  for  each  case  the  figure  where  it  is  graphed  as  well  as  its  set  of 
characterizing  bulk  parameters.  In  Appendix  B,  graphs  are  drawn  three  per 
page  in  order  of  low  to  high  Ad  for  each  A  ranging  from  most  negative  to 
most  positive.  The  most  populous  classification  cells  were  in  two  groups: 
low  to  moderate  Ad  with  low  to  moderate  A ,  and  high  Ad  with  moderate  to 
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Table  2 

Parameters  for  Example  Spectra 


Rgure 

H^o 

m 

Hz 

deg 

deg 

■ 

AS 

deg 

3 

1.41 

0.064 

44 

72 

0.895 

1.25 

2.16 

103 

1.50 

Bit’ 

4.37 

0.064 

78 

71 

0.391 

0.21 

6.42 

27 

-0.43 

Bim 

3.53 

0.083 

118 

97 

0.555 

-0.20 

3.75 

43 

-0.26 

B1b 

4.06 

0.064 

64 

66 

0.485 

2.10 

8.13 

20 

0.10 

B2t 

4.53 

0.083 

96 

101 

0.468 

0.01 

5.38 

28 

0.18 

B2m 

3.85 

0.083 

100 

99 

0.518 

-0.21 

4.16 

35 

-0.07 

B2b 

3.18 

0.083 

62 

87 

0.560 

0.46 

3.70 

45 

0.13 

B3t 

2.26 

0.113 

70 

95 

0.654 

0.09 

2.78 

56 

0.09 

B3m 

3.86 

0.064 

68 

74 

0.492 

1.65 

7.27 

21 

0.29 

B3b 

4.71 

0.083 

90 

99 

0.441 

0.29 

5.86 

28 

0.30 

B4t 

3.40 

0.083 

54 

60 

0.498 

0.46 

4.12 

36 

0.29 

B4m 

3.38 

0.074 

88 

102 

0.600 

0.30 

3.37 

46 

0.54 

B4b 

3.30 

0.083 

92 

96 

0.724 

-0.05 

3.09 

56 

0.23 

B5t 

3.96 

0.064 

66 

73 

0.518 

1.74 

5.94 

24 

0.64 

B5m 

3.47 

0.064 

88 

98 

0.639 

1.80 

4.91 

30 

0.80 

B5b 

3.54 

0.064 

84 

97 

0.625 

1.41 

4.27 

37 

0.92 

B6t 

3.01 

0.074 

60 

81 

0.832 

1.90 

3.38 

48 

0.64 

B6m 

3.40 

0.093 

52 

83 

0.695 

0.83 

2.88 

60 

0.60 

B6b 

1.67 

0.083 

52 

97 

0.847 

0.65 

1.91 

96 

0.93 

B7t 

2.24 

0.083 

52 

72 

0.677 

1.73 

3.85 

43 

1.12 

B7m 

2.52 

0.074 

64 

82 

0.796 

1.95 

3.40 

47 

1.13 

B7b 

1.45 

0.083 

58 

92 

0.858 

1.06 

2.21 

IB 

B8t 

0.81 

0.093 

58 

92 

0.896 

1.22 

2.31 

88 

m\ 

B8m 

2.00 

0.064 

58 

74 

0.927 

1.22 

2.29 

1.26 

B8b 

1.68 

0.113 

40 

83 

0.915 

1.21 

2.22 

99 

IB 

B9t 

1.58 

0.113 

44 

84 

0.927 

1.16 

2.04 

105 

1.26 

B9m 

0.61 

0.103 

52 

92 

0.978 

0.84 

1.77 

113 

1.22 

B9b* 

1.62 

0.132 

40 

74 

0.920 

1.53 

2.32 

99 

IB 

^The  letters  t,  m,  and  b  indica 
in  Appendix  B. 

te  top,  middle,  an 

i  bottom  graphs. 

espectiv* 

3ly,  in  figures 
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high  A .  Consequently,  the  graphs  in  Appendix  B  appear  somewhat  scattered. 
However,  the  figures  presented  and  discussed  in  this  chapter  are  for  illustra¬ 
tive  rather  than  analytic  purposes,  so  the  apparent  lack  of  system  in  their 
order  is  not  detrimental. 

The  example  spectra  illustrate  many  properties  described  parametrically  in 
Chapters  5  and  6.  For  instance,  all  examples  indicate  significant  amoimts  of 
wave  energy  arriving  from  the  northwest  (wave  directions  between  0  and 
90  deg  in  the  figures)  as  indicated  in  Figure  10  for  statistical  distributions  of 

and  in  Figure  20  for  correlations  of  6^  and  .  Example  cases  with 
>  4  m  (13.1  ft)  (Figures  Bit,  Bib,  B2t,  and  B3b,  where  t,  m,  and  b 
indicate  top,  middle,  and  bottom  graphs,  respectively,  in  the  figures  in  Ap¬ 
pendix  B)  all  have  A6  between  20  and  30  deg  as  indicated  in  Figure  25,  and 
thus  give  a  graphic  image  of  narrowly  spread,  high-energy  spectra. 

Directionally  broad  examples  (Figures  3,  B6b,  B7b,  and  all  graphs  in  Fig¬ 
ures  B8  and  B9)  are  characteristically  wind  seas  from  the  northwest  on  which 
are  superimposed  peaks  of  energy  from  the  south.  The  added  effect  of  these 
southern  energy  peaks  is  to  increase  dramatically  estimates  of  A6 .  This  can 
be  seen  qualitatively  by  comparing  Figures  B3b  and  B8b,  where  the  wind  sea 
(from  the  northwest)  parts  of  the  spectra  have  roughly  the  same  features,  but 
Figure  B8b  shows  other  peaks  of  energy  from  the  south  and  at  low  frequency. 
For  Figure  B3b,  Ad  =  28  deg,  which  is  a  reasonable  estimate  for  the  wind  sea 
part  of  Figure  B8b.  However,  Figure  B8b  has  a  bulk  A6  of  99  deg,  which 
must  be  mostly  due  to  the  secondary  energy  mode  fi-om  the  south. 

There  can  be  substantive  differences  in  detailed  spectral  structure  within 
classes  characterized  by  ranges  of  Ad  and  A .  For  example.  Figures  B9t 
(Ad  =  105  deg,  A  =  1.26)  and  B9m  (Ad  =  113  deg,  A  =  1.22)  are  from  the 
same  class  cell,  the  two  having  been  included  in  the  example  set  for  compari¬ 
son  of  contrasting  energy  levels;  the  spectrum  of  Figure  B9t  represents  more 
than  six  times  the  energy  of  the  spectrum  in  Figure  B9m.  Figure  B9m  indi¬ 
cates  greater  directional  energy  spread  at  high  frequencies  and  greater  fre¬ 
quency  distribution  of  energy  in  the  southern  source  mode  than  the  case 
shown  in  Figure  B9t.  However,  the  integrated  direction  spectra  D(d^)  shown 
in  the  left  rear  panels  of  these  graphs,  and  firom  which  the  classifying  bulk 
parameters  are  derived,  have  comparable  structural  shapes.  Subtleties  in 
spectral  structure  are  not  all  energy-related.  Figures  3  and  B9b  are  from  the 
same  class  cell,  and  have  nearly  the  same  set  of  parameters  throughout,  except 
for  peak  frequency.  Visually,  it  appears  that  the  structure  of  the  wind  sea 
mode  of  energy  fi'om  the  northwest  is  substantially  narrower  in  Figure  B9b 
than  in  Figure  3.  Many  of  these  subtleties  can  be  resolved  by  analyzing  and 
characterizing  separately  each  energy  mode  in  a  spectrum.  Though  beyond 
the  scope  of  the  initial  investigation  reported  here,  such  analysis  clearly  can 
yield  a  more  detailed  and  uniform  characterization  of  these  spectra. 

Disregarding  minor  differences,  the  example  spectra  provide  images  of 
energy  distributions  as  they  are  identified  through  the  classifying  parameters. 
Spectra  with  narrow  directional  spreads  tend  to  be  unimodal  with  small  to 
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moderate  asymmetries.  Cases  in  the  high-energy  regime  tend  to  have  the 
narrowest  directional  spreads.  Cases  with  large  directional  spreads  tend  to  be 
multimodal  with  varying  contributions  of  energy  from  southern  sources  to 
extant  wind  seas.  Though  there  exist  many  additional  spectra,  some  with 
more  deviate  energy  distributions,  the  examples  used  in  this  chapter  are  from 
populous  directional  spread  and  asymmetry  classes,  and  so  are  representative 
of  much  of  the  data  set. 
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8  Characteristic  Directionai 
Distribution  Functions 


A  way  to  gain  further  insight  in  the  way  wave  energy  is  directionally  dis¬ 
tributed  in  observations  from  Harvest  Platform  is  to  examine  the  distribution 
functions  from  all  frequencies  from  all  observed  frequency-direction  spectra. 
Each  of  the  2,339  frequency-direction  spectra  considered  in  this  report  has 
13  frequency  bands  so  that  from  the  normalized  frequency-direction  spectra 
D(f^,0m)  defined  by  Equation  7  there  are  a  total  of  30,407  individual  direc¬ 
tional  distribution  functions.  These  functions  can  be  classified,  grouped,  and 
averaged  within  groups  to  define  characteristic  directional  distribution  func¬ 
tions.  Because  normalized  distributions  are  used,  any  effect  of  energy  content 
is  removed  in  this  type  of  analysis.  Because  all  distribution  functions  are 
averaged  without  regard  to  source  wave  frequency,  any  knowledge  of  frequen¬ 
cy  effects,  such  as  high-frequency  versus  low-frequency  processes,  is  lost. 

This  is  acceptable  at  the  level  of  analysis  done  here  because  the  objective  is  to 
define  general  characteristic  shapes,  along  with  some  knowledge  of  frequency 
of  occurrence  of  these  shapes. 


Classification  Algorithms 

In  this  report,  directional  distribution  functions  are  classified  and  grouped 
based  on  three  firequency-dependent  parameters  ,  A0„ ,  and  defined  in 
Chapter  4.  For  each  frequency  f„  in  a  normalized  frequency-direction  spec¬ 
trum  D(f^,  d^) ,  the  directional  distribution  function  is  translated  on  the  direc¬ 
tion  axis  by  its  mean  direction  6^^ ,  classified  as  being  in  a  particular  range  of 
directional  spread  Ad^  and  asymmetry  ,  and  stored  in  a  classification  buffer 
function  D.j^(d^  -^q)-  Mathematically,  this  process  is  expressed 

-  «»)  =  -ej  m  .  1, 2 . M  (50: 


Index  i  in  Equation  50  is  an  index  designating  a  range  of  directional  spread, 
and  is  defined  as  the  nearest  integer  to  A6J 10  so  that  the  i'*  range  includes 
directional  spreads  from  lOi  -  5  deg  to  lOi  +  5  deg.  Most  of  the 
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observations  were  classifiable  when  twelve  spread  ranges  were  used  so  that  i 
assumes  values  from  1  to  12.  Index  j  in  Equation  50  designates  a  range  of 
quartile  as5mimetry,  and  is  defined  as  the  nearest  integer  to  AJQA  so  that  the  j* 
range  includes  asymmetries  from  0.4y  -  0.2  to  0.4y  +  0.2 .  The  bulk  of  ob¬ 
servations  had  asymmetries  in  the  overall  range  of  -2.2  to  2.2  so  that  j  as¬ 
sumes  values  from  -5  to  5,  with  the  class  for  which  y  =  0  being  the  most 
symmetric.  Index  k  in  Equation  50  is  a  case  index  to  distinguish  among  the 
contributing  distribution  functions  in  the  accumulation  buffer.  The  upper  limit 
for  k  in  each  spread-asymmetry  classification  cell  is  K.j ,  which  indicates  the 
total  number  of  contributors  to  cell  ij.  The  argument  6^  on  the  left-hand  side 
of  Equation  50  is  a  generic  notation  that  indicates  all  contributors  to  the  accu¬ 
mulation  buffer  are  centered  on  their  mean  directions. 

In  this  classification  scheme,  30,044  directional  distributions  out  of  a  possi¬ 
ble  30,407  distributions  fit,  with  only  363  distributions  (1.2  percent)  exceed¬ 
ing  the  indicated  ranges  of  spread  or  asymmetry.  Furthermore,  as  will  be 
seen,  not  all  of  the  classification  cells  had  contributions.  Nevertheless, 

98.8  percent  of  all  observations  are  accounted  for,  and  the  classification  cells 
defined  here  form  a  regular  grid  in  the  spread-asymmetry  plane. 

Characteristic  distribution  functions  are  defined  as  averages  (at  each  dis¬ 
crete  direction)  of  contributing  distributions  in  each  classification  cell.  A 
measure  of  variability  of  the  members  within  a  cell  is  found  by  computing 
standard  deviations  of  contributing  distributions.  Such  statistics  require  a 
minimum  number  of  contributions  in  each  cell.  In  this  study,  no  statistics 
were  computed  for  any  cell  in  which  <  3.  In  cells  for  which  K..'>3, 
sums  were  computed  of  the  value  and  the  square  of  each  contributing  member 
in  the  accumulation  buffer  for  each  discrete  direction.  The  mathematical 
expressions  for  these  processes  are 


*=1 


and 


w  =  1, 2, ... ,  M 


/n  =  1, 2, ... ,  M 


(51) 


(52) 


A  characteristic  mean  distribution  function  -  0„)  is  then  defined  as 


=  m=l,2,...,M 


(53) 
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and  a  standard  deviation  -d^)  of  contributing  distributions  about  the 

mean  is  computed  from 


-  1 


w  =  1,  2, ... ,  M 


(54) 


Equation  53  provides  the  essential  definition  for  the  set  of  characterizing 
distribution  fonctions.  Equation  54  provides  a  measure  of  how  well  contribut¬ 
ing  distributions  agree  in  structural  form.  If  standard  deviation  values  are 
small  relative  to  means,  it  is  an  indication  that  contributing  shapes  are  quite 
similar.  If  standard  deviations  are  large,  there  is  considerable  variability  in 
the  forms  of  contributing  distributions. 


Resulting  Distributions 

Characteristic  distributions  resulting  from  these  computations  are  illustrated 
in  Appendix  C.  There  is  a  graph  representing  each  spread-asymmetry  class, 
with  a  page  of  graphs  for  each  class  of  directional  spread,  and  each  page 
shows  all  the  asymmetry  classes  for  the  indicated  directional  spread  class. 

Each  graph  shows  D^(6^-Bq)  as  a  solid  line.  Dashed  lines  show  standard 

deviations  in  the  forms  +  D^.  and  D..  -  D^j .  Each  graph  also  indicates  the 
number  of  contributing  distributions  used  in  the  computations.  If  a  cell  had 
too  few  contributions  (Ky  <  3),  the  number  of  cell  members  is  shown,  but  no 
curves  are  drawn. 

Shapes  of  the  mean  curves,  relative  values  of  standard  deviation  curves, 
and  numbers  of  contributing  distributions  are  all  relevant  to  interpreting  these 
graphs.  Characteristic  distributions  with  small  to  moderate  spreads  and  small 
to  neutral  asymmetry  shown  as  graphs  near  page  centers  in  Figures  Cl  to  C4 
tend  to  have  simple  shapes  with  a  strong  primary  peak  and  either  no  or  a  very 
slight  side  lobe  of  energy.  Characteristic  curves  with  broad  spreads  and  small 
to  neutral  asymmetries,  as  in  Figures  CIO  and  Cll,  appear  to  be  composed  of 
two  or  more  main  lobes  of  comparable  energy  content.  Curves  with  the 
largest  asymmetries,  which  do  not  appear  in  the  narrowest  classes,  but  begin 
with  Figure  C5,  all  have  well-displaced  (in  direction)  secondary  lobes  of 
moderate  relative  intensity.  The  classification  scheme  appears  to  preserve 
distribution  symmetry;  within  a  spread  class,  curves  in  as3mimetry  classes  of 
the  same  magnitude  but  of  opposite  sign  appear  qualitatively  to  be  mirror 
images  of  each  other.  This  property  holds  true  even  though  there  tends  to  be 
a  great  many  more  members  in  classes  with  positive  asymmetry  than  with 
negative. 

Quality  of  the  characterizing  nature  of  these  curves  appears  to  be  higher 
for  distributions  in  classes  with  small  to  moderate  spreads  than  for  those  with 
large  spreads.  In  Figures  Cl  to  C3,  standard  deviations  are  small  relative  to 
mean  values  near  distribution  peaks.  This  condition  suggests  that  all 
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contributing  distributions  have  nearly  the  same  shape.  In  contrast,  Figures 
CIO  to  C 12  show  distributions  with  standard  deviations  that  are  larger  per¬ 
centages  of  peak  values,  suggesting  considerable  variability  in  contributing 
shapes. 

However,  a  pattern  seems  evident  in  the  largest  spread  classes  shown  in 
Figures  CIO  to  C12.  Most  of  the  mean  distributions  appear  to  be  bimodal, 
and,  with  careful  examination,  it  can  be  seen  that  modal  separation  increases 
as  spread  classification  parameter  increases.  This  result  suggests  that  distribu¬ 
tions  assigned  to  the  higher  spread  classes  are  products  of  wave  conditions 
having  multiple  sources  (as  indicated  by  examples  in  Chapter  7)  so  that  spread 
classification  is  controlled  by  directional  separation  of  source  processes  rather 
than  energy  distributions  within  individual  modes.  In  fact,  if  modes  are  sepa¬ 
rated  at  minimum  points  in  the  high-asymmetry  distributions  in  Figures  CIO 
to  C12,  it  appears  that  individual  modes  would  be  classifiable  within  the  low- 
to  moderate-spread  classes  of  Figures  Cl  to  C5.  Though  clearly  a  valuable 
topic  for  future  research,  such  modal  separation  is  not  done  in  this  study 
because  it  is  important  to  retain  complete  distributions  for  overall  data  charac¬ 
terization,  and  for  estimation  accuracy  presented  in  Chapter  9. 

Not  so  obvious  in  the  figures  in  Appendix  C  is  the  relative  importance  of 
some  classes  over  others.  All  classes  are  shown,  which  suggests  equal  weigh¬ 
ting,  but  some  classes  are  much  more  populous  than  others,  and  so  are  more 
representative  of  common  or  t)q)ical  conditions  at  Harvest  Platform.  The 
single  most  populous  class  is  the  neutrally  symmetric  cell  with  spreads  near 
30  deg  (Figure  C3).  With  3,295  contributing  distributions,  this  one  class 
represents  approximately  10.8  percent  of  all  distributions.  In  contrast,  the 
classes  with  spreads  of  100  deg  or  greater  (Figures  CIO  to  C12)  have  only 
between  5  and  124  cases,  or  0.02  to  0.41  percent  of  all  observations. 

Clearly,  the  classes  with  small  to  intermediate  spreads  (Figures  Cl  to  C5) 
and  near-neutral  asymmetry  are  dominant.  For  the  group  of  classes  with 
Ad  <  55  deg  and  -0.2  <  A  <  0.2,  9,582  distributions  or  31.5  percent  of  all 
distributions  are  represented.  If  the  first  asymmetry  cells  on  either  side  of  the 
neutral  asymmetry  cells  are  included  (i.e.,  -0.6  <  A  <  0.6)  in  the  group  with 
Ad  <  55  deg,  the  number  of  distributions  rises  to  20,312,  or  66.8  percent  of 
all  cases,  clearly  a  dominant  group  of  classes.  Much  less  dominant  are  the 
high-spread  classes.  The  whole  set  of  classes  with  Ad  >  65  deg  (all  cases  in 
Figures  C7  through  C12)  contains  only  4,198  cases  or  13.8  percent  of  obser¬ 
vations.  Note  that  if  individual  modes  are  analyzed  separately,  as  mentioned 
previously,  many  members  of  the  high-spread  classes  would  add  to  the  low- 
spread,  low-asymmetry  classes. 

Some  contributions  to  the  high-spread  classes  may  be  less  due  to  wave 
signal  than  to  background  system  noise.  Included  in  the  classification  set  are 
all  distributions  from  the  lowest  analyzed  frequency  (0.044  Hz).  As  indicated 
in  Figure  27,  the  lowest  wave  frequency  is  dominated  by  very  low  wave 
energy.  In  low-energy  conditions,  collection  and  analysis  noise  become  more 
important.  A  signature  of  noisy  conditions  in  directional  analysis  is  very 
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broad  distributions.  Because  only  13  resolution  frequency  bands  were  ana¬ 
lyzed,  a  significant  part  of  one-thirteenth  (7.7  percent)  of  contributing  distri¬ 
butions  could  be  affected  by  noise.  This  condition  could  be  eliminated  in 
future  analysis  by  constraining  contributing  distributions  to  be  from  spectral 
densities  above  a  certain  minimum  level.  The  main  results  of  the  present 
work  are  not  seriously  affected  by  a  small  percentage  of  potentially  noisy 
results,  however,  because  such  results  reside  in  classes  with  low  populations. 

The  characteristic  directional  distribution  functions  shown  in  Appendix  C 
and  discussed  in  this  chapter  are  elucidative  in  that  they  reveal  some  of  the 
variety  of  common  shapes  of  energy  distributions  observed  at  Harvest  Plat¬ 
form.  The  more  common  shapes  are  readily  identified  by  populations  within 
the  various  classification  cells.  These  characteristic  shapes  are  also  extremely 
useful  in  evaluating  the  resolving  ability  of  the  Harvest  Platform  directional 
wave  gauge,  as  presented  in  the  following  chapter. 
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9  Accuracy  of  Directional 
Distribution  Estimates 


Central  to  any  gauging  device  is  the  accuracy  with  which  it  resolves  an 
input  signal.  This  is  especially  true  in  high-resolution  directional  wave  gaug¬ 
ing  because,  in  the  application  used  here,  an  array  of  separate  pressure  gauges 
is  employed,  and  the  computational  algorithm  is  complicated.  The  IMLE 
method  is  nonlinear,  which  complicates  error  estimates  relative  to  linear  sys¬ 
tems,  where  errors  can  often  be  described  analytically.  The  IMLE  method  is 
also  data  adaptive,  which  means  results  are  quite  sensitive  to  input  conditions. 
At  present,  the  only  known  way  to  estimate  directional  resolving  ability  with 
this  method  is  to  make  computations,  basing  input  on  known  directional  distri¬ 
butions  and  comparing  IMLE  output  estimates  to  these  known  conditions. 
Clearly,  the  choice  of  input  conditions  is  important,  and  should  be  closely 
related  to  conditions  at  a  given  observation  site.  For  the  Harvest  Platform 
data  set  discussed  in  this  report,  reasonable  test  conditions  are  the  characteris¬ 
tic  distribution  functions  described  in  Chapter  8.  The  present  chapter  de¬ 
scribes  test  methodology  and  resulting  estimates  of  uncertainty  for  these  test 
conditions. 


Test  Method 

If  a  discrete  directional  distribution  function  at  a  given  frequency 

is  known,  an  exact  (to  discretization  accuracy)  normalized  cross-spectral 
matrix  can  be  computed  as  a  basis  for  input  for  IMLE  estimation. 

The  mathematical  expression  is  like  Equation  12,  taking  the  form 

=  E  (55) 


where  element  indexes  i  and  j  range  from  1  to  7 ,  the  number  of  gauges  in 
the  test  array.  If  the  left-hand  side  of  Equation  55  is  used  in  place  of  the 
data-derived  cross-spectral  matrix  used  in  Equation  9,  then  the  set  of  Equa¬ 
tions  9  to  15  can  be  used  to  make  an  IMLE  estimate  of  D(f^,  6^) ,  and  the  two 
distribution  functions  can  be  compared. 
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This  procedure  is  correct,  and  is  often  used,  but  is  somewhat  unrealistic  in 
that  it  overlooks  the  natural  variability  of  a  sea  surface.  Displacement  of  a 
natural  sea  surface  is  commonly  treated  as  a  Gaussian  random  process,  which 
means  there  is  expected  to  be  noise  in  any  cross  spectrum  that  represents  real 
wave  signals.  Consequently,  noise  must  be  added  to  the  exact  cross  spectrum 
of  Equation  55  to  simulate  natural  variability. 

In  tests  reported  here,  noise  was  added  to  exact  cross  spectra  using  the 
method  of  Brennan  and  Mallett  (1976).  In  their  method,  the  of  Equa¬ 

tion  55  are  elements  of  the  7x7  Hermitian  matrix  M,  which  is  decomposed 
using  standard  matrix  procedures  as 

M  =  UA  I/-' 


where  1/  is  an  7  x  7  matrix  with  each  column  holding  one  eigenvector  of  M, 
A  is  a  diagonal  matrix  containing  the  7  eigenvalues  of  M,  and  the  superscript 
-1  indicates  the  matrix  inverse.  In  the  following  expressions,  use  is  made  of 
the  facts  that  the  eigenvector  matrix  satisfies 


where  the  asterisk  indicates  complex  conjugation,  and  the  subscript  T  denotes 
matrix  transpose,  and,  for  a  Hermitian  matrix,  eigenvalues  of  M  are  all  real. 
Using  these  relationships,  the  right-hand  side  of  Equation  56  can  be  rewritten 
as 


M  =  17A‘°  A*^  U-' 

=  17  A*'^  A*''  U;  (58) 

=  V  Vj. 


where 

V  =  U*  A*'" 


(59) 


The  advantage  of  this  manipulation  is  that  noise  can  be  introduced  by  multi¬ 
plying  the  matrix  V  by  a  column  vector  Z  of  7  complex  random  deviates 
from  a  Gaussian  distribution  with  zero  mean  and  variance  i  to  form  a  matrix 

with  noise  V  in  the  form 


V  =  vz 

=  77*  A^^Z 


(60) 
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Equation  60  can  be  substituted  into  the  right-hand  side  of  the  last  equality  of 
Equation  58  to  form  a  cross-spectral  matrix  with  noise  as 


M  =  V*  f  j. 

=  UA}'^Z*  Zj.A^'^Ut 


(61) 


where  it  will  be  remembered  that  matrix  product  order  is  reversed  on  transpo¬ 
sition  for  the  second  triplet  of  terms  in  the  second  equality  of  Equation  61. 

Equation  61  has  two  desired  properties  associated  with  adding  noise  to  a 

cross-spectral  matrix.  One  property  is  that,  on  average,  M  is  equal  to  M, 
which  can  be  shown  by  finding  the  expected  value  (the  average  over  a  large 
number  of  observations)  of  Equation  61.  Because  the  eigenvalue  and  eigen¬ 
vector  matrices  are  constant  for  a  given  simulation,  the  expectation  operation 
can  be  written 

E[M]  =  E[U\^'^Z*  Z^A.^'^U;] 

=  UA^'^E[Z*  Zj.]A^'^Ur 
=  UA''UA^^^Ur 
=  UA^'^A^'^U; 

=  M 

where  £■[  ]  is  the  expectation  operator,  and  I  is  the  identity  matrix.  The 
expected  value  of  Z*  Zj.  becomes  identity  because,  as  written,  each  of  the 

complex  random  numbers  in  Z  is  multiplied  by  its  own  conjugate,  forming 
the  sum  of  the  squares  of  two  Gaussian  random  deviates.  Each  of  the  squared 

values  has  an  expected  value  of  1 ,  so  they  sum  to  form  an  expected  value  of 
unity.  The  other  desired  property  is  that  the  noise  is  distributed  throughout 
the  matrix  elements  in  exactly  the  same  way  that  cross-spectral  estimates  are 
distributed  when  computed  from  raw  Fourier  coefficients  of  data  from  the  I 
gauges  of  the  test  array. 

The  cross  spectrum  with  noise  computed  from  Equation  61  is  appropriate 
for  raw  cross-spectral  estimates,  and  has  two  degrees  of  freedom.  To  simu¬ 
late  results  from  a  data  collection,  results  from  Equation  61  must  be  comput¬ 
ed,  accumulated,  and  averaged  a  number  of  times  corresponding  to  the  prod¬ 
uct  of  the  number  of  ensembles  used  and  the  number  of  frequency  bands 
averaged  in  an  actual  data  collection.  In  collections  analyzed  for  this  report, 
data  were  analyzed  using  eight  contiguous  ensembles  with  fi-equency  smooth¬ 
ing  over  ten  frequency  bands.  Consequently,  an  average  of  80  computations 
of  Equation  61,  each  with  new  random  numbers,  was  required  to  simulate  the 
160  degrees  of  freedom  obtained  in  prototype  data  collections. 
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If  the  elements  Mg(f„)  of  the  resulting  cross  spectrum  with  noise  are  used 
as  the  source  cross  spectrum  at  a  given  test  frequency  in  Equations  9  to  15, 
the  resulting  IMLE  directional  distribution  function  estimate  is  but 

one  of  a  large  number  (indexed  by  subscript  of  possible  results  as  a  conse¬ 
quence  of  the  added  noise.  To  find  an  expected  result  from  these  possible 
results,  the  procedure  outlined  above  for  finding  one  estimate  is  repeated  a 
large  number  of  times,  and  the  results  are  accumulated  and  averaged.  For  the 
tests  reported  here,  experimentation  revealed  that  reasonably  stable  statistics 
were  obtained  if  500  estimates  were  computed.  An  expected  distribution 

D(fn,0^)  and  standard  deviation  from  the  expected  distribution  for 

a  given  test  case  were  found  by  first  computing  the  sums 


500 

'"-i.z . M 

k=\ 


and 


500 


k=\ 


(63) 


(64) 


and  then  computing 


D(fM  = 


WnA) 

500 


m  =  1,2, ...  ,M 


(65) 


and 


= 


s,(fM-500D^(fn^ej 

499 


ffi  =  1, 2, ...  ,M 


(66) 


It  is  noted  that  the  statistics  computed  using  Equations  63  to  66  are  strictly 
legitimate  only  for  MLE  results,  where  variance  estimates  in  each  direction  6^ 
are  independent  of  each  other.  In  IMLE,  the  constraints  of  Equations  8  and 
16  result  in  some  coupling  of  results  across  all  directions.  Consequently, 
mean  and  standard  deviation  distributions  from  Equations  65  and  66,  respec¬ 
tively,  are  good  indicators  of  estimator  accuracy  but  are  not  completely  defen¬ 
sible  from  first  principles. 


Test  Distributions 

To  examine  resolvability  of  some  of  the  variety  of  directional  distributions 
that  occur  at  Harvest  Platform,  the  mean  characteristic  distributions  discussed 
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in  Chapter  8  and  illustrated  in  Appendix  C  were  used  as  test  distributions. 
These  distributions  have  the  advantages  that  they  are  representative  of  proto¬ 
type  conditions,  being  derived  as  averages  of  such  conditions,  and  they  are  of 
reasonable  number.  This  last  feature  is  important  because  it  takes  a  signifi¬ 
cant  amount  of  computer  time  to  run  a  test  condition  through  the  test  proce¬ 
dure  outlined  in  the  previous  section.  Testing  a  great  number  of  cases  then 
becomes  somewhat  impractical.  Within  the  directional  spread-asymmetry  cells 
shown  in  Appendix  C,  there  are  110  cells  for  which  mean  distributions  are 
shown.  These  1 10  mean  distributions  constitute  the  basic  test  set. 

The  1 10  test  distributions  were  defined  by  averaging  individual  distribu¬ 
tions  without  regard  to  the  frequencies  at  which  the  individual  distributions 
were  observed.  However,  resolving  ability  of  a  fixed  array  of  gauges  is  a 
function  of  the  wavelengths  of  incident  waves,  which  are  related  to  frequen¬ 
cies  by  virtue  of  the  dispersion  relation  (Equation  11).  To  conduct  a  more 
complete  set  of  tests,  each  of  the  test  distributions  was  run  for  three  frequen¬ 
cies.  These  frequencies  were  0.05,  0.10,  and  0.15  Hz  to  represent  low-, 
medium-,  and  high-frequency  regions,  respectively,  in  observed  spectra.  With 
1 10  test  shapes  run  at  each  of  three  frequencies,  there  were  a  total  of  330  test 
conditions. 


Test  Results 

Test  results  discussed  in  the  present  chapter  are  illustrated  in  Appendix  D. 
Results  varied  reasonably  smoothly  through  the  set  of  test  conditions  so  that 
some  space  is  saved  by  illustrating  a  subset  of  test  cases.  Shown  in  Appen¬ 
dix  D  are  19  figures,  each  representing  one  of  the  110  test  distributions.  Each 
figure  shows  results  for  all  three  of  the  test  frequencies.  Represented  are  tests 
based  on  characteristic  distributions  from  directional  spread  classes  centered  at 
10  deg  (the  narrowest  class),  20,  30,  and  40  deg  (the  most  populous  classes), 
60  and  100  deg  (intermediate  classes),  and  120  deg  (the  broadest  class). 

Within  each  directional  spread  class  except  the  narrowest,  illustrated  cases 
include  neutral  asymmetry  as  well  as  an  extreme  positive  asymmetry  and  an 
extreme  negative  asymmetry.  Only  the  neutral  asymmetry  case  is  shown  for 
the  narrowest  spread  class. 

Each  figure  contains  three  graphs,  one  for  each  frequency.  The  test  fre¬ 
quency  is  shown  in  the  upper  right-hand  comer  of  each  graph,  along  with  the 
wavelength  for  waves  of  that  frequency  in  water  depth  of  200  m.  All  graphs 
within  a  figure  are  drawn  to  the  same  scale,  but  scales  change  from  figure  to 
figure.  Each  graph  contains  three  curves.  A  solid  line,  called  the  tme  distri¬ 
bution,  is  the  test  distribution  obtained  from  the  characteristic  distribution 
illustrated  in  Appendix  C  for  the  ranges  of  and  A  indicated  at  the  tops  of 
the  figures  in  Appendix  D.  Dotted  lines  are  expected,  or  mean,  distributions 
based  on  500  estimates  as  expressed  by  Equation  65.  Ideally,  expected  and 
true  distributions  are  identical.  When  they  are  not  identical,  differences  indi¬ 
cate  the  amount  of  estimator  bias.  Dashed  lines  show  standard  deviation 
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distributions  as  computed  using  Equation  66,  and  provide  a  measure  of  the 
variability  of  the  set  of  estimates  from  the  expected  distribution. 

The  set  of  figures  in  Appendix  D  indicate  that  the  best  results  are  obtained 
for  distributions  with  narrow  spread  and  near-neutral  asymmetry  in  medium- 
to  high-frequency  wave  conditions.  Figures  Dl,  D3,  and  D6  all  illustrate 
very  small  differences  between  true  and  expected  distributions  for  test  fre¬ 
quencies  of  0. 10  and  0.15  Hz.  For  low-frequency  waves  at  test  frequency 
0.05  Hz,  there  are  differences  of  approximately  5  to  10  percent  between  the 
peaks  of  the  true  and  expected  distributions.  Low-frequency  waves  are  typi¬ 
cally  more  difficult  to  resolve  in  IMLE  because  their  wavelengths  (approxi¬ 
mately  600  m  in  these  examples)  are  much  larger  than  the  maximum  size  of 
the  array  (approximately  60  m,  as  shown  in  Figure  2)  so  that  phase  differenc¬ 
es  due  to  small  changes  in  wave  direction  are  more  difficult  to  resolve  than 
for  shorter  waves.  Note  that  because  the  total  areas  under  directional  distri¬ 
bution  functions  are  constrained  to  be  unity  by  Equation  8,  regions  where 
expected  distributions  are  lower  than  true  distributions  must  be  compensated 
by  regions  at  other  directions  where  expected  distributions  exceed  true  distri¬ 
butions.  Thus,  where  the  peak  expected  curve  is  lower  than  the  true  curve  for 
the  lowest  frequency  in  Figure  Dl,  the  shoulders  of  the  expected  curve  have 
higher  values  than  the  true  distribution. 

Standard  deviation  curves  appear  uniformly  maximum  in  the  regions  of  the 
maxima  of  the  true  curves.  In  the  regions  of  these  maxima,  standard  devia¬ 
tions  are  approximately  10  to  20  percent  of  the  true  distribution  values.  This 
behavior  indicates  that  in  a  single  estimate  using  data  from  real  waves,  a 
variation  of  the  maximum  by  10  to  20  percent  is  not  uncommon.  At  direc¬ 
tions  different  from  true  spectral  maxima,  the  standard  deviation  curves  be¬ 
come  smaller,  though  it  appears  there  are  increases  in  percentages  of  true 
spectral  values  represented  by  the  standard  deviation  curves.  For  example,  at 
the  lowest  frequency  illustrated  in  Figure  Dl,  the  standard  deviations  in  the 
directions  near  20  deg  and  -20  deg  appear  to  be  about  half,  or  50  percent,  of 
the  true  distribution  values.  The  figures  suggest  that  these  higher  percentages 
occur  in  regions  away  from  spectral  peaks,  and  so  may  be  less  dire  errors. 

Test  distributions  with  moderate  to  large  asymmetry  and  large  spread  are 
slightly  less  well-resolved  than  are  the  narrow,  neutrally  symmetric  cases. 

The  magnitude  of  asymmetry  does  not  necessarily  indicate  poorer  results, 
however.  For  example,  expected  distributions  for  low  and  intermediate  fre¬ 
quencies  in  the  narrow,  moderately  asymmetric  example  of  Figure  D2  show  a 
5-  to  10-deg  bias  of  peak  values  toward  the  high-energy  tail  of  the  distribu¬ 
tion,  along  with  small  variations  in  the  resolution  of  the  high-energy  tail.  In 
contrast.  Figure  D16  illustrates  a  case  with  spread  of  approximately  100  deg 
and  quartile  asymmetry  of  about  2.0,  yet  has  about  the  same  magnitude  bias 
toward  its  high-energy  tail,  and  likewise  small  variations  in  resolution  of  the 
tail  region. 

The  poorest  representation  of  true  distribution  by  expected  distribution  is 
shown  in  Figure  D18,  which  is  a  case  having  neutral  symmetry,  but  extreme 
spread  with  Ad  near  120  deg.  Both  true  and  expected  distributions  have 
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curiously  undulating  values  in  the  regions  of  their  broad  peaks.  At  none  of 
the  test  frequencies  are  the  agreements  exceptional;  expected  peak  directions 
deviate  from  true  peak  directions  by  15  to  30  deg.  A  reason  for  this  behavior 
is  that  the  IMLE  algorithm  was  derived  to  try  to  isolate  narrowly  distributed 
wave  signals  from  background  noise  (Davis  and  Regier  1977),  which  it  does 
quite  well  as  shown  in  the  example  of  Figure  Dl.  In  broadly  distributed  wave 
signals  like  those  of  Figure  D 18,  the  algorithm  tends  to  produce  spurious 
peaks,  which  can  be  seen  as  central  maxima  for  the  two  higher  frequencies  in 
Figure  D18.  Qualitatively,  the  agreement  of  true  and  expected  distributions  in 
Figure  D18  is  quite  good.  The  main  features  of  the  true  distribution  (broad, 
rather  uniform  high-energy  region,  and  nearly  constant  background  region) 
are  reasonably  well-represented,  but  not  as  well  as  in  narrower  examples. 

It  is  noted  that  all  the  test  distributions  taken  from  the  characteristic  shapes 
of  Appendix  C  have  mean  directions  near  0  deg.  This  means  that  the  centroid 
directions  are  those  of  waves  propagating  from  the  north,  and  so  only  propa¬ 
gate  through  one  orientation  of  the  array  geometry  shown  in  Figure  2.  In 
more  complete  tests,  true  distributions  would  be  shifted  so  their  centroids 
were  at  several  places  around  the  horizon,  and  thus  test  array  response  to 
varied  wave  angle  of  attack.  This  was  impractical  because  it  would  add  an¬ 
other  dimension  to  an  already  computationally  intensive  test  program. 

However,  there  are  indications  in  the  present  test  set  that  varying  wave 
angle  of  attack  does  not  substantially  change  array  response.  In  all  of  the 
asymmetric  examples  in  Appendix  D  for  Ad  >  35  deg  (Figures  D8,  DIO, 

Dll,  D13,  D14,  D16,  D17,  and  D19),  there  are  clear  indications  of  side 
lobes  having  substantial  fractions  of  total  energy.  Existence  of  these  side 
lobes  causes  distribution  mean  directions  to  shift  away  from  distribution  peaks 
so  that  peak  energy  conditions  are  not  at  0  deg,  but  deviate  by  as  much  as 
35  deg  (Figures  D17  and  D19).  Secondary  lobes  deviate  from  0  deg  by  as 
much  as  110  deg  (Figure  D14).  Because  distributions  with  both  positive  and 
negative  asymmetries  were  tested,  and  test  results  were  nearly  symmetric 
(compare  Figures  D14  and  D16,  for  example),  one  can  conclude  that  direc¬ 
tional  response  is  generally  uniform  for  all  waves  coming  from  the  northeast 
and  northwest  quadrants.  By  symmetry,  responses  to  waves  coming  from  the 
southeast  and  southwest  quadrants  will  be  the  same  because  their  orientation  to 
the  array  geometry  is  the  same  as  for  waves  from  the  two  northern  quadrants. 
Thus,  these  tests  suggest  that  array  response  is  reasonably  uniform  for  all 
wave  approach  directions. 

It  can  be  argued  that  test  conditions  were  somewhat  biased  because  obser¬ 
vations  from  the  Harvest  Platform  array  were  used  to  test  array  response,  and 
it  is  reasonable  to  expect  such  a  device  to  be  able  to  resolve  its  own  output 
quite  well.  However,  in  no  test  case  was  an  exact  array  result  used  as  a  test 
condition;  all  test  conditions  were  derived  from  averages  of  several  array 
results.  Because  the  standard  deviations  of  the  characteristic  distributions  of 
Appendix  C  were  finite,  it  is  unlikely  that  any  average  distribution  was  an 
exact  replication  of  any  single  result.  Furthermore,  the  addition  of  synthetic 
statistical  noise  ensured  that  no  two  input  conditions  were  identical.  Finally, 
because  the  expected  distributions  were  very  much  like  the  true  conditions  for 
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all  of  the  most  populous  characteristic  distributions  tested  (Figures  D2  to 
DIO),  it  is  likely  that  the  characteristic  distributions  used  for  these  tests  were 
largely  representative  of  actual  wave  conditions  at  Harvest  Platform. 

These  tests  indicate  that  the  directional  wave  response  of  the  Harvest  Plat¬ 
form  array  is  quite  good  for  any  application  that  can  tolerate  slight  biases  in 
peak  energy  levels  and  occasional  errors  in  peak  direction  of  up  to  10  deg. 
Results  from  the  array  may  only  be  qualitatively  useful  for  highly  direction- 
sensitive  models  like  the  one  described  by  O’Reilly  and  Guza  (1993).  Howev¬ 
er,  these  results  should  be  excellent  ground  truth  for  such  applications  as  deep 
ocean  wave  modeling,  remote  sensing  by  air-  and  spacecraft,  and  field  com¬ 
parisons  of  low-resolution  directional  wave  gauges.  Further  indications  of  the 
quality  of  array  response  are  provided  by  examination  of  the  behavior  of 
directional  parameters  that  characterize  these  data.  Such  an  examination  is 
done  in  Chapter  10. 
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10  Accuracy  of  Directional 
Parameters 


Characterizing  parameters  of  directional  distributions  of  wave  energy  rather 
than  the  distributions  themselves  is  a  common  practice  in  engineering  and 
modeling  efforts.  A  set  of  seven  frequency-dependent  parameters  that  de¬ 
scribe  various  features  of  directional  distributions  were  introduced  in  Chap¬ 
ter  4  of  this  report,  and  discussed  in  statistical  form  in  Chapter  6.  These 
parameters  are  distribution  peak  direction  6^^ ,  the  circular  moment  parameters 
describing  mean  direction  ,  circular  width  ,  skewness  7^ ,  and  kurtosis 
6^ ,  and  the  quartile  parameters  describing  directional  spread  ,  and  asym¬ 
metry  A^.  To  determine  how  well  these  parameters  can  be  resolved  in  typical 
wave  conditions  with  the  Harvest  Platform  array,  statistics  were  accumulated 
during  the  directional  distribution  tests  described  in  Chapter  9.  Summary 
statistics  of  these  parameters  are  displayed  and  discussed  here. 

Statistical  descriptors  were  simply  the  means  and  standard  deviations  of 
each  of  the  parameters  during  the  500  iterative  runs  of  the  110  test  distribu¬ 
tions  for  each  of  the  three  test  frequencies.  The  parameter  mean  value  indi¬ 
cates  its  expected  value  from  any  arbitrary  test  sample,  and  can  be  compared 
to  the  parameter  true  value  to  investigate  biases.  Standard  deviations  give  a 
measure  of  variability  of  test  results,  and  indicate  a  nominal  confidence  range 
for  a  parameter  in  a  given  test  condition. 


Peak  Direction 

Figure  36  illustrates  results  of  these  tests  for  peak  direction.  There  is  one 
graph  for  each  of  the  three  test  frequencies.  Each  graph  is  in  the  form  of  a 
correlation  of  true  with  mean  test  B^  (the  subscript  n ,  which  indicates  a 
frequency-dependent  parameter,  is  dropped  in  subsequent  discussion  because 
frequency  is  shown  explicitly  in  each  graph).  Superimposed  on  each  mean  B^ 
point  is  the  standard  deviation  of  the  test  set  in  the  form  of  vertical  error  bars 
that  extend  above  and  below  the  mean  point  by  one  standard  deviation.  Sym¬ 
bols,  identified  at  the  top  of  Figure  36,  indicate  the  directional  spread  class  of 
the  test  distributions,  and  correspond  with  the  spread  classes  of  the  character¬ 
istic  distributions  shown  in  Appendix  C.  For  each  symbol,  there  is  a  data 
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point  for  each  of  the  asymmetry  classes  associated  with  a  given  spread  class  in 
Appendix  C.  For  example,  circles  represent  the  narrowest  spread  class  in  the 
test  set,  and  there  are  four  circles  in  each  graph  in  Figure  36  because  there 
are  only  four  asymmetry  classes  shown  in  Figure  Cl  for  this  spread  class. 

Distinction  of  symbols  is  important  because  the  eye  is  drawn  to  the  large 
error  bars  in  Figure  36,  yet  the  large  error  bars  do  not  represent  the  most 
common  cases  in  the  observed  data.  Close  examination  of  Figure  36  indicates 
that  the  cases  with  greatest  standard  deviations  are  mostly  from  classes  with 
large  spreads  and  large  6^ .  Mean  results  for  these  cases  are  also  biased  to¬ 
ward  small  peak  directions,  in  some  cases  by  more  than  60  deg.  Curiously, 
this  is  a  desired  result.  Examination  of  the  test  distributions  for  cases  with 
extreme  spread  (Figure  C12)  shows  that  peak  directions  are  extreme  for  the 
cases  with  near-neutral  asymmetry.  Because  such  distributions  have  broad, 
rather  flat  peaks,  minor  variations  in  distributions  from  the  test  cases  can  be 
detected  as  distribution  peaks  scattered  over  a  wide  range  of  directions.  On 
average,  one  would  expect  such  a  distribution  to  fall  near  the  centers  of  the 
broad  flat  regions,  yielding  peak  directions  near  0  deg.  In  the  extreme  cases 
of  Figure  36,  this  is  exactly  what  happens,  yielding  what  appear  to  be  large 
biases,  and  very  large  standard  deviations.  The  three  cases  of  extreme  6^ 
from  each  of  the  four  broadest  cases  (diamond-plus,  circle-plus,  double¬ 
triangle,  and  box-plus  symbols)  all  exhibit  this  behavior  for  all  test  frequencies 
shown  in  Figure  36.  Another  large-error  group  of  cases  centered  near  0  deg 
is  from  classes  of  intermediate  spread  with  near-neutral  asymmetry  that  also 
have  considerable  uncertainty  in  location  of  peak  direction,  but  for  which 
peaks  of  test  spectra  are  near  0  deg. 

If  these  extreme  cases  are  eliminated  from  Figure  36,  much  tighter  group¬ 
ings  of  points  emerge.  The  group  with  the  smallest  variations  is  the  set  of 
narrowest  cases  (open-circle  symbols),  for  which  the  standard  deviation  indi¬ 
cations  are  about  the  same  size  as  the  plotted  symbol  (approximately  1  or 
2  deg).  The  next  broader  group  (base-down-triangle  symbols)  has  slightly 
larger  scatter  (approximately  3  to  4  deg),  and  an  occasional  bias  of  about  the 
same  size.  Note  that  none  of  these  narrow  cases  extends  to  very  large  and 
so  are  clustered  within  10  to  20  deg  of  0  deg. 

A  general  interpretation  of  Figure  36  is  that  errors  in  peak  direction  are 
small  for  narrow  distributions  or  broader  distributions  with  well-defined, 
narrow  modal  peaks.  When  modal  peaks  become  broad  and  flat,  peak  direc¬ 
tion  varies  dramatically  due  to  minor  output  fluctuations  that  vary  with  differ¬ 
ent  combinations  of  statistical  noise.  Comparison  of  the  three  graphs  in  Fig¬ 
ure  36  indicates  somewhat  smaller  ranges  of  error  bar  lengths  for  cases  with 
narrow  and  intermediate  spreads  at  the  highest  frequency.  This  is  a  conse¬ 
quence  of  improved  IMLE  resolution  of  high-frequency  distributions  as  dis¬ 
cussed  in  Chapter  9. 
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Circular  Moment  Parameters 


Peak  direction  might  be  expected  to  be  noisy  because  it  is  a  single-point 
parameter,  relying  on  directional  estimates  in  a  particular  direction  to  be 
stable  in  the  presence  of  naturally  varying  input  conditions.  Moment-based 
parameters  are  expected  to  be  less  noisy  because  moments  are  weighted  sums 
of  contributions  from  all  directions  in  a  distribution,  a  process  that  acts  as  a 
stabilizing  filter.  In  this  section,  the  behaviors  of  the  four  frequency- 
dependent  parameters  based  on  distribution  circular  moments  are  examined. 

Mean  direction 

Tests  of  mean  direction  reported  here  are  somewhat  limited  because  all  the 
characteristic  distributions  of  Appendix  C  that  were  used  for  testing  were 
centered  on  their  mean  directions.  Thus,  results  from  all  statistical  tests 
should  ideally  yield  the  same  result.  Figure  37  shows  the  means  and  standard 
deviations  of  all  test  cases  of  plotted  not  as  functions  of  true  0^,  which 
should  be  constant,  but  rather  as  functions  of  true  asymmetry  parameter  A , 
which  may  reveal  sensitivity  to  irregular  energy  distributions.  As  with  peak 
direction.  Figure  37  contains  three  graphs,  one  for  each  test  frequency.  Sym¬ 
bols,  defined  at  the  top  of  the  figure,  denote  directional  spread  classes  of  test 
distributions  shown  in  Appendix  C.  Horizontal  lines  through  each  graph  show 
the  value  of  true  0g  for  all  the  test  cases.  Due  to  an  indexing  error  in  an 
analysis  program,  all  test  distributions  were  centered  one  direction  bin  too 
high,  so  that  true  6^  =  -2  deg  instead  of  0  deg.  Note  that  this  error  has  no 
effect  on  any  other  tests  or  discussion  in  this  report,  and  only  shifts  the  refer¬ 
ence  lines  in  Figure  37  to  -2  deg  instead  of  0  deg. 

For  a  given  spread  and  asymmetry  class  shown  in  Figure  37,  the  smallest 
errors  occur  at  the  highest  test  frequency,  and  increase  modestly  for  the  two 
lower  test  frequencies.  For  a  given  asymmetry  class,  error  bars  are  smallest 
for  the  narrowest  test  cases  and  increase  with  increasing  spread  range.  For 
example,  for  the  cases  at  which /=  0.15  Hz,  and  A  «  0,  the  standard  devia¬ 
tion  is  less  than  1  deg  for  the  narrowest  spread  class  (open-circle  symbol),  but 
is  more  than  6  deg  for  the  broadest  spread  class  (box-plus  symbol).  Error 
bars  are  roughly  the  same  size  for  all  asymmetries  in  a  given  spread  class  and 
at  a  particular  fi'equency. 

At  the  highest  frequency  shown  in  Figure  37,  there  is  very  little  bias  in  the 
mean  test  0g.  The  s3mibols  tend  to  be  within  + 1  deg  of  the  reference  line  at 
this  frequency.  For  the  middle  frequency,  there  is  a  slight  bias  that  appears 
correlated  with  asymmetry.  This  bias  is  about  -3  deg  for  the  largest  negative 
asymmetry,  is  near  0  deg  at  neutral  asymmetry,  and  becomes  positive  by 
about  2  deg  for  the  largest  positive  asymmetry.  A  similar,  but  stronger  pat¬ 
tern  is  seen  in  the  graph  for  the  lowest  frequency,  where  the  bias  ranges  from 
about  -5  deg  at  extreme  negative  asymmetry  to  about  3  deg  at  extreme  posi¬ 
tive  asymmetry.  The  bias  is  due  to  the  fact  that,  at  low  frequencies,  there  is  a 
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tendency  for  IMLE  to  yield  too  little  energy  under  a  true  distribution  peak, 
and  consequently,  too  much  energy  in  the  tails.  When  weighted  by  the  trigo¬ 
nometric  functions  in  its  definition  (Equation  40),  the  resulting  6^  is  shifted 
away  from  the  peak  and  toward  the  high-energy  tail.  This  behavior  makes  it 
correlate  with  asymmetry,  for  which  the  high-energy  tail  changes  sides  as 
asymmetry  changes  sign. 

An  overall  characterization  of  mean  direction  accuracy  is  represented  by 
the  graph  of  the  center  frequency  in  Figure  37.  A  nominal  standard  deviation 
is  +5  deg.  A  nominal  bias  is  2  to  3  deg,  depending  on  asymmetry.  If  more 
exact  characterization  of  error  is  necessary,  the  information  plotted  in  Fig¬ 
ure  37  can  be  used. 


Circular  width 

Means  and  standard  deviations  of  the  sets  of  a  estimates  are  shown  as 
functions  of  true  circular  width  in  Figure  38,  following  the  same  class  identifi¬ 
cation  patterns  used  in  Figures  36  and  37.  Figure  38  shows  a  to  be  relatively 
robust,  with  no  strong  biases  evident,  and  reasonably  uniform  standard  devia¬ 
tions  within  each  test  frequency  group.  Though  errors  are  modestly  smaller 
for  the  highest  test  frequency,  it  is  relatively  straightforward  to  assign  a  nomi¬ 
nal  uncertainty  of  ±0.05  to  estimates  of  a . 


Circular  skewness 

Like  the  circular  width  parameter,  circular  skewness  7  is  reasonably  well- 
behaved.  Figure  39  illustrates  the  correlation  of  estimated  mean  skewness 
with  true  skewness,  along  with  vertical  error  bars  to  show  standard  deviations 
of  estimates.  For  each  test  frequency,  the  standard  deviations  are  nominally 
the  same,  with  slight  reductions  in  the  vicinity  of  7  =  0.  As  usual,  errors  are 
least  at  high  frequencies,  and  increase  with  decreasing  frequency.  Nominal 
standard  deviations  are  ±0.4  for  /=  0.15  Hz,  and  ±0.5  for  the  other  two  test 
frequencies.  Biases  are  evident  in  the  results  at  extreme  values  of  7,  with  the 
bias  directions  being  toward  smaller  7.  For  7  »  4,  a  nominal  bias  of  -0.3 
can  be  associated  with  the  lowest  test  frequency,  -0.2  with  the  center  test  fre¬ 
quency  and  less  than  -0. 1  for  the  high-frequency  test.  Biases  are  very  small 
near  7  =  0,  and  are  largest  for  the  narrowest  classes  of  test  cases  (open-circle, 
base-down-triangle,  ±,  and  x  symbols),  suggesting  sensitivity  to  minor  varia¬ 
tions  in  skewed,  but  closely  grouped  distributions  of  high-energy  points.  This 
parameter  is  better  resolved  for  directionally  broad,  highly  asymmetric  distri¬ 
butions  than  for  any  of  the  narrow,  moderately  asymmetric  distributions. 


Circular  kurtosis 

Test  results  for  circular  kurtosis  b  are  shown  in  Figure  40.  Standard 
deviations  of  results  tend  to  increase  with  increasing  b .  A  nominal  estimate 
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Figure  38.  Test  results  for  circular  width 
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Figure  40.  Test  results  for  circular  kurtosis 
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of  standard  deviation  based  on  Figure  40  is  15  percent  of  8 .  Errors  are 
somewhat  smaller  at  the  highest  test  frequency.  Biases  are  also  evident  in 
Figure  40  for  the  intermediate  and  low  test  frequencies.  Biases  are  low  by 
approximately  8  percent  for  /=  0.05  Hz,  and  5  percent  for  /=  0.10  Hz.  As 
with  skewness,  the  largest  errors  are  associated  with  the  narrowest  classes  of 
test  distributions.  These  distributions  have  high  kurtosis,  meaning  they  are 
more  peaked  than  Gaussian  curves,  and  have  increased  scatter  in  their  kurtosis 
estimation  because  kurtosis  is  a  higher  moment  of  a  distribution,  and  the 
energetic  parts  of  these  distributions  are  defined  by  a  small  number  of  points. 
Higher  moments  are  typically  more  difficult  to  resolve  than  lower  moments 
like  mean  direction  and  circular  width.  Absolute  errors  are  smallest  for  low 
kurtosis  distributions,  which  are  from  the  very  broad  spread  classes. 


Quartile  Parameters 


Two  additional  parameters  tested  are  quartile  spread  Ad  and  asymmetry  A , 
which  are  based  on  quartile  directions  of  distribution  fimctions.  Because 
quartile  points  are  determined  through  integration  of  a  distribution  function 
(Equation  44),  these  parameters  are  expected  to  have  stabilities  comparable  to 
those  of  circular  moment  parameters. 


Quartile  spread 

Figure  41  shows  correlations  of  mean  test  Ad  with  distribution  true  Ad, 
along  with  vertical  error  bars  showing  standard  deviations  of  test  spread  pa¬ 
rameters.  Because  these  graphs  are  of  Ad ,  it  is  pointless  to  classify  results  by 
spread  grouping,  so  symbols  are  employed  to  represent  asymmetry  classes, 
and  are  defined  at  the  top  of  Figure  41 .  Error  bars  tend  to  be  small  for  small 
spreads,  being  less  than  3  deg  for  the  narrowest  spread  class  at  all  test  fre¬ 
quencies.  Error  bars  increase  with  increasing  Ad,  reaching  nominal  values  of 
5  to  15  deg  for  Ad  near  60  deg,  and  then  are  approximately  constant  for 
higher  Ad .  High-frequency  distributions  tend  to  have  smaller  errors,  with 
nominal  values  of  ±5  deg  for  large  Ad ,  than  the  two  lower-frequency  groups. 

There  is  a  slight  bias  toward  smaller  spreads  evident  in  the  intermediate- 
and  low-frequency  test  sets,  being  approximately  5  deg  at  the  highest  Ad. 

This  bias  is  due  to  the  tendency  for  IMLE  to  shift  estimates  of  modal  peaks  in 
bimodal  distributions  such  that  the  peaks  move  closer  together.  This  has  the 
effect  of  compressing  the  distribution  along  the  direction  axis,  resulting  in 
lower  estimates  of  A6 . 

A  nominal  value  for  imcertainty  in  A0  is  ±  10  deg  based  on  the  data  set  as 
a  whole,  though  this  is  clearly  an  overestimate  of  the  errors  for  small  Ad . 
Note  that  this  nominal  uncertainty  is  qualitatively  consistent  with  uncertainty 
in  circular  width  a,  which  is  also  a  measure  of  directional  spread.  Based  on 
the  general  slope  of  the  correlation  curve  of  Figure  6  relating  bulk  estimates 
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of  these  parameters,  a  variation  of  +0.05  in  a  is  approximately  equivalent  to 
a  variation  of  ±  10  deg  in  Ad . 


Quartile  asymmetry 

The  final  parameter  examined  in  these  tests  is  quartile  asymmetry  A ,  for 
which  test  results  are  illustrated  in  Figure  42.  Correlations  of  test  A  with  true 
A  are  shown  for  three  test  frequencies,  and  symbols  denote  spread  classes. 
Error  bars,  denoting  standard  deviations  of  test  A,  seem  of  reasonably  uni¬ 
form  size  within  each  graph.  Nominal  standard  deviation  values  are  0.3,  0.4, 
and  0.5  for  /=  0.15,  0.10,  and  0.05  Hz,  respectively.  The  middle-  and  low- 
frequency  results  reveal  biases  at  large  A  toward  smaller  values.  The  reason 
for  these  biases  is  that  the  IMLE  tendency  to  underestimate  peak  values  in 
wide,  bimodal  distributions  has  the  effect  of  making  estimated  distributions 
more  uniform  along  the  direction  axis,  thereby  reducing  their  asymmetry. 


Summary  of  Parameter  Accuracy 


Table  3  is  a  siunmary  listing  of  nominal  parameter  accuracies  and  bias 
ranges  for  the  circular  and  quartile  parameters  based  on  information  shown  in 
Figures  37  to  42.  The  intent  of  Table  3  is  to  provide  citable  accuracies  that 
represent  the  bulk  of  typical  observations  at  Harvest  Platform.  Clearly,  the 
table  is  not  complete  because  many  of  the  parameters  have  greater  or  lesser 
variability  than  the  cited  numbers  as  can  be  seen  by  examining  the  test  result 
figures.  Of  note  is  the  fact  that  tests  of  mean  direction  were  extremely  limit¬ 
ed,  having  been  conducted  for  only  a  single  test  direction.  Results  cited  in 
Table  3  are  based  on  that  single  test  direction,  but  may  not  be  representative 
of  the  broader  suite  of  mean  directions  from  which  waves  arrive  at  this  gauge 


Table  3 

Accuracies  and  Biases  of  Characterizing  Parameter! 

5 

Parameter 

Test  Range 

Nominal  Accuracy 

Nominal  Bias  Range 

So 

-2  deg’ 

±5  deg 

-5  deg  to  3  deg 

a 

0.38  to  1.17 

±0.05 

None 

Y 

-1.6  to  3.8 

±0.5 

-0.3  to  0.1 

6 

1.4  to  13 

±  1 5  percent 

-8  percent  to  None 

Ad 

1 2  deg  to  1 20  deg 

±  1 0  deg 

-5  deg  to  None 

A 

-2.0  to  2.0 

±0.4 

-0.3  to  0.3 

I  ^All  tests  were  run  with  a  constant  mean  direction  || 
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Figure  42.  Test  results  for  quartile  asymmetry 
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Missing  from  Table  3  are  characterizing  accuracy  and  biases  for  peak 
direction.  This  parameter  was  omitted  because  test  results  were  highly  vari¬ 
able,  and  could  not  be  represented  well  by  single  numbers.  Figure  36  must 
be  used  to  estimate  accuracy  and  bias  for  peak  direction. 

It  should  be  remembered  that  the  test  results  presented  here  are  based  on 
the  addition  of  synthetic  Gaussian  noise  to  IMLE  input  cross-spectra  derived 
from  characteristic  directional  distributions.  The  noise  was  theoretically  ap¬ 
propriate  for  160  degrees  of  freedom,  which  matches  the  degrees  of  freedom 
claimed  in  the  basic  directional  analysis  described  in  Chapter  3.  Results  using 
longer  or  shorter  data  records  than  those  used  in  this  report  would  have  great¬ 
er  or  fewer  degrees  of  freedom,  and  consequent  differences  in  parameter 
resolution. 

An  advantage  of  evaluating  accuracies  and  biases  of  directional  parameters 
is  that  it  allows  a  quantification  of  the  ability  of  an  array  to  resolve  directional 
distributions  of  wave  energy,  which  are  difficult  otherwise  to  evaluate.  In 
combination  with  distribution  test  results  discussed  in  Chapter  9,  the  paramet¬ 
ric  evaluation  in  the  present  chapter  indicates  the  extent  to  which  the  Harvest 
Platform  array  can  resolve  directional  distributions,  and  reasons  for  inaccura¬ 
cies  and  biases.  With  this  knowledge,  a  user  can  then  decide  if  reported 
results  are  adequate  for  a  given  purpose.  Because  there  are  no  other  long¬ 
term,  deepwater,  high-resolution  directional  wave  gauges  (as  of  this  writing), 
the  present  results  are  the  best  in  existence.  They  should  add  immensely  to 
imderstanding  deep  ocean  wave  conditions,  which  will  allow  better  definitions 
of  seaward  boundary  conditions  for  applied  models  of  shoreward  propagating 
wave  energy.  Furthermore,  the  existence  of  the  Harvest  Platform  gauge  (or, 
at  least,  the  use  at  other  sites  of  the  technical  principles  described  in  this 
report)  should  allow  better  understanding  of  results  from  commonly  used  low- 
resolution  directional  gauges,  and  improved  ground  truth  for  tests  of  applied 
remote  sensing  techniques. 
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1 1  Summary 


A  high-resolution  directional  wave  gauge  has  been  installed  on  the  Texaco 
Oil  Company  Harvest  Platform  to  initiate  long-term  monitoring  of  the  direc¬ 
tional  wind  wave  climate  at  a  deepwater  site  that  can  be  used  to  represent 
open  ocean  conditions  for  waves  approaching  the  coast  of  southern  California. 
In  addition  to  wind  wave  climatology,  results  from  this  gauge  can  be  used  for 
seaward  boundary  conditions  in  coastal  wave  propagation  models,  studies  of 
ocean  wave  evolution,  comparison  studies  with  locally  deployed  low- 
resolution  directional  wave  gauges,  and  ground  truth  in  remote  sensing  re¬ 
search  and  tests.  This  report  describes  various  attributes  of  the  directional 
wave  gauge,  and,  in  some  detail,  a  summary  of  results  from  one  year  of 
operation. 


Gauge  and  Basic  Data  Set 

Harvest  Platform  is  located  approximately  20  km  (10.8  n.m.)  west  of  Point 
Conception,  California,  at  about  the  200-m  (656-ft)  water  depth.  The  direc¬ 
tional  gauge  is  an  array  of  six  pressure  gauges  mounted  on  the  Harvest  Plat¬ 
form  structure  at  a  nominal  depth  of  16  m  (53  ft).  Utilizing  the  maximum 
available  horizontal  extent  of  the  platform  structure,  the  array  is  rectangular  in 
shape  with  outside  dimensions  approximately  23  by  59  m  (75.5  by  193.6  ft). 
Data  are  acquired  eight  times  daily  in  collections  of  time  series  having  dura¬ 
tions  equal  to  2  hr  16  min,  with  all  six  gauges  sampled  simultaneously  at 
1  Hz.  Data  quality  is  assured  by  intercomparing  frequency  spectra  from  all 
six  gauges,  and  further  by  intercomparing  frequency  cross  spectra  over  redun¬ 
dant  lag  spacings  within  the  array.  Data  that  passed  these  quality  checks  were 
further  processed  into  frequency-direction  spectra  using  the  iterative  maximum 
likelihood  estimation  method  derived  and  described  by  Pawka  (1983).  This 
report  is  based  on  all  results  from  calendar  year  1993,  from  which  2,339 
satisfactory  collections  were  processed  from  a  total  of  2,920  possible  collec¬ 
tions. 
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Bulk  Parameters 


A  set  of  characterizing  parameters  were  defined  and  used  to  reduce  the 
substantial  information  in  the  frequency-direction  spectra  to  manageable  form. 
In  addition  to  characteristic  wave  height  and  peak  frequency  normally  defined 
for  wave  data,  a  set  of  direction-dependent  parameters  were  computed  both 
for  bulk  integrated  direction  spectra  (directional  analogs  of  frequency  spectra), 
and  for  directional  distributions  at  individual  frequencies.  The  parameter  set 
included  peak  direction,  four  moment  parameters  defined  by  Kuik,  van 
Vledder,  and  Holthuijsen  (1988),  and  two  quartile  parameters  defined  by  Long 
and  Oltman-Shay  (1991).  The  four  moment  parameters  were  mean  direction, 
circular  width,  skewness,  and  kurtosis.  The  two  quartile  parameters  were 
directional  spread  and  asymmetry.  Additionally,  directional  and  frequency 
modes  were  analyzed  to  obtain  number  and  relative  importance  of  peaks  in 
uni-  and  multimodal  distributions. 

Bulk  parameters  were  used  to  characterize  the  general  behavior  of 
frequency-direction  spectra  as  whole  units.  Time  series  of  these  parameters, 
illustrated  in  Appendix  A,  indicated  wave  climatology  for  the  collection  year. 
Statistical  distributions  of  parameters  indicated  common  and  extreme  values 
that  were  observed.  Characteristic  wave  heights  ranged  from  0.6  m  (2.0  ft)  to 
more  than  5  m  (16.4  ft)  with  common  values  between  1  and  2  m  (3.3  and 
6.6  ft).  Most  common  peak  frequency  was  0.07  Hz,  suggesting  a  predomi¬ 
nance  of  swell,  but  all  frequency  bands  were  occasionally  most  energetic, 
indicating  the  occurrence  of  initial  and  well-developed  stages  of  local  storms 
as  well  as  long-travelled  waves  from  distant  sources.  Peak  directions  were 
mostly  due  to  waves  from  the  northwest,  but  about  5  percent  of  observations 
were  from  the  south,  suggesting  occasional  dominance  of  the  spectra  by  south¬ 
ern  swell.  Mean  direction  demonstrated  a  similar  pattern,  but,  being  an  inte¬ 
gral  parameter,  was  only  shifted  slightly  toward  the  south,  indicating  a  com¬ 
mon  coexistence  of  local  wind  seas  and  southern  swell. 

Circular  width  and  directional  spread  both  showed  two  main  groupings  of 
cases.  A  narrow  group  with  circular  width  commonly  0.55  and  directional 
spread  30  deg  was  associated  with  unimodal  spectra  from  single  sources.  A 
broad  group  with  circular  width  peaking  between  0.9  and  1.0,  and  directional 
spread  having  a  maximum  near  100  deg  was  attributed  to  multimodal  seas 
from  directionally  divergent  sources.  Circular  skewness  and  asymmetry  pa¬ 
rameters  both  showed  common  positive  values,  indicating  that  energy  distri¬ 
butions  were  not  symmetric,  but  had  energetic  tail  regions  preferentially  coun¬ 
terclockwise  from  main  energy  peaks.  Circular  kurtosis  was  dominated  by 
values  greater  than  3,  indicating  distribution  shapes  sharper  or  more  peaked 
than  Gaussian  curves.  Approximately  24  percent  of  kurtosis  values  were  less 
than  3,  suggestive  of  flattening  that  would  occur  either  in  broad  unimodal 
distributions  or  multimodal  sources  widely  separated  in  direction. 

Most  of  the  bulk  distributions  were  either  uni-  or  bimodal.  Frequency 
spectra  were  approximately  61  percent  unimodal  and  38  percent  bimodal. 
Directional  spectra  were  only  39  percent  unimodal,  and  60  percent  bimodal. 
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Secondary  modes  were  not  inconsequential.  Approximately  25  percent  of 
directional  distributions  had  25  percent  of  total  energy  in  second  modes. 

Correlations  of  bulk  parameters  with  characteristic  wave  height  gave  indi¬ 
cations  of  spectral  behavior  in  various  wave  conditions,  especially  in  energetic 
conditions  that  may  be  of  concern  to  coastal  interests.  High-energy  waves 
typically  had  peak  frequencies  in  the  range  0.07  to  0.09  Hz,  and  both  peak 
and  mean  directions  in  the  west-northwest.  Directional  spectra  of  energetic 
waves  were  narrow  with  circular  widths  between  0.35  and  0.50,  and  direc¬ 
tional  spreads  from  15  to  45  deg.  Circular  skewness  and  asymmetry  were 
both  near  zero,  with  a  tendency  toward  slightly  positive  values.  Circular 
kurtosis  achieved  some  of  its  highest  values  in  high-energy  conditions,  indicat¬ 
ing  very  peaked  energy  distributions. 


Frequency-Dependent  Parameters 

More  detailed  information  is  obtained  by  examining  statistical  distributions 
of  characterizing  parameters  on  a  frequency-by-frequency  basis.  This  type  of 
analysis  avoids  the  smoothing  effect  of  integration  across  all  frequencies  re¬ 
quired  for  definition  of  the  directional  spectrum  from  which  the  bulk  parame¬ 
ters  were  derived.  At  the  frequency  level,  there  appear  three  groupings  with 
common  features:  very  low  frequencies,  identified  by  the  single  band  at 
0.044  Hz,  low  frequencies  in  the  range  0.054  to  0.093  Hz,  and  middle  to  high 
frequencies,  which  extend  to  0.16  Hz,  the  highest  frequency  examined. 

Very  low-frequency  distributions  typically  have  very  low  energy,  with  peak 
and  mean  directions  in  the  west-northwest,  and  across  much  of  the  southwest 
quadrant.  Circular  width  and  directional  spread  extend  over  broad  ranges, 
with  a  having  common  values  from  0.5  to  1.2,  and  M  ranging  from  less 
than  20  deg  to  more  than  100  deg.  Asymmetry  and  circular  skewness  param¬ 
eters  indicate  such  distributions  to  be  highly  symmetric,  and  kurtosis  is  com¬ 
monly  equal  to  2.0,  suggesting  broad,  flat  distributions.  Approximately 
96  percent  of  all  such  distributions  are  unimodal.  Because  energy  levels  are 
low,  distributions  at  this  frequency  are  likely  to  be  mostly  noise,  which  tends 
to  be  directionally  white. 

Low-frequency  distributions  contain  broad  ranges  of  energy,  including  the 
highest  energies  of  any  group  of  frequencies.  Mean  and  peak  directions  have 
two  common  groupings:  one  from  the  northwest  quadrant,  and  another,  with 
about  equal  commonality,  from  the  south.  This  behavior  is  indicative  of  the 
interplay  between  well-developed  western  storm  waves  and  swell  from  distant 
southern  sources.  Circular  width  shows  broad  ranges  of  common  values  from 
0.4  to  p.9,  and  directional  spread  is  not  only  widely  distributed,  but  also 
shows  two  main  groupings  for  this  range  of  frequencies.  One  group,  from  10 
to  60  deg  in  spread,  suggests  unimodal,  locally  generated  seas.  The  other 
group,  with  spreads  from  60  to  120  deg,  indicates  multimodal  distributions 
with  energy  from  widely  separated  directions.  Skewness  and  asymmetry  both 
show  a  gradual  shift  of  common  values  from  near  zero  to  slightly  positive  as 


96 


Chapter  1 1  Summary 


frequency  increases  through  this  band.  Kurtosis  has  a  wide  range  of  common 
values  from  2  to  6  for  these  frequencies,  again  suggesting  a  variety  of  distri¬ 
bution  shapes  from  broad  to  peaky.  This  range  of  frequencies  is  noted  for  its 
high  number  of  directional  modes,  as  many  as  six,  with  secondary  modes 
commonly  energetic  because  90  percent  or  more  of  total  energy  is  contained 
in  the  primary  mode  for  only  50  to  70  percent  of  the  observations. 

High-frequency  distributions  tend  to  be  more  regular,  with  intermediate 
levels  of  energy,  peak  and  mean  directions  primarily  centered  in  the  northwest 
quadrant,  and  modest  angular  distributions,  with  a  typically  0.4  to  0.7,  and 
A0  in  the  range  10  to  60  deg,  peaking  near  30  to  40  deg.  These  distributions 
tend  to  be  neutrally  to  slightly  positively  skewed,  though  there  are  modest 
ranges  of  these  parameters  with  y  between  -1  and  3,  and  A  between  -1.0 
and  1.5.  Kurtosis  ranges  from  2  to  18  at  these  frequencies,  but  has  common 
values  of  4  to  5,  indicating  relatively  peaked  shapes.  The  most  common 
number  of  directional  modes  for  these  distributions  is  two,  but  the  effect  of 
the  secondary  modes  is  less  than  for  the  low  frequencies  because  approximate¬ 
ly  80  percent  of  all  observations  have  90  percent  of  energy  in  primary  modes. 


Representative  Results 

To  illustrate  the  appearance  of  observed  frequency-direction  spectra,  27 
examples  are  shown  in  Appendix  B.  These  examples  are  primarily  high- 
energy  cases  because  such  cases  are  typically  of  most  interest,  though  a  few 
low-  and  intermediate-energy  cases  are  shown  for  comparison.  The  examples 
were  chosen  as  representative  of  various  bulk  directional  spreads  and  as5mime- 
tries  by  which  the  entire  observation  set  was  classified.  Examples  shown  are 
mostly  the  highest-energy  case  from  each  of  the  most  populous  classification 
cells. 

Of  the  30,407  individual  directional  distribution  functions  observed  during 
1993,  110  characteristic  distributions  were  defined  to  illustrate  variety  and 
frequency  of  occurrence  of  typical  shapes.  All  individual  distributions  were 
classified  and  grouped  in  cells  defined  by  ranges  of  quartile  spread  and  asym¬ 
metry.  All  members  of  each  cell  were  then  averaged,  without  regard  to 
source  wave  frequency,  to  determine  characteristic  shapes.  Standard  devia¬ 
tions  of  each  cell  population  were  also  computed  to  indicate  shape  variability 
within  a  class.  Illustrated  in  Appendix  C,  the  most  common  classes  by  far 
were  those  with  narrow  to  moderate  spreads  and  small  to  neutral  asymmetries. 
Broad  distributions  with  multiple  modes  were  not  common.  This  behavior 
suggests  that  where  a  frequency-direction  spectrum  can  be  considered  multi¬ 
modal,  the  cause  is  a  series  of  mostly  unimodal  distributions  at  different  fre¬ 
quencies,  such  as  local  wind  sea  from  the  northwest  at  intermediate  frequen¬ 
cies  coupled  with  southern  swell  at  low  frequencies. 
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Gauge  Accuracy 


To  investigate  the  directional  resolving  ability  of  the  Harvest  Platform 
array  and  its  analysis  algorithm,  a  number  of  tests  were  run  using  characteris¬ 
tic  distributions  as  basic  test  elements.  Ideal  cross  spectra  were  computed  for 
these  test  cases,  and  then  synthetic,  Gaussian  noise  appropriate  to  the  compu¬ 
tational  degrees  of  freedom  was  added  to  the  cross  spectra.  When  used  as 
input  to  the  IMLE  method,  a  cross  spectrum  with  noise  resulted  in  an  output 
distribution  that  represented  one  possible  result  from  a  natural  sea.  This 
process  was  repeated  500  times  for  each  test  case  so  that  statistics  could  be 
computed  from  a  population  of  500  results.  All  110  of  the  characteristic 
distributions  were  treated  this  way,  and,  because  results  are  frequency- 
sensitive,  all  distributions  were  tested  for  three  representative  frequencies: 
0.05,  0.10,  and  0.15  Hz. 


From  the  500-member  population  of  test  results  for  each  case  and  each  test 
frequency,  mean  and  standard  deviation  distributions  were  computed  to  deter¬ 
mine  an  expected  result  and  a  measure  of  variability  of  results.  Some  of  the 
test  results  are  illustrated  in  Appendix  D,  where  it  is  seen  that  results  are 
excellent  for  high-frequency  waves  having  distributions  with  low  to  moderate 
spreads  with  any  as)unmetry.  Results  degraded  somewhat  with  decreasing 
frequency  for  the  low-  to  moderate-spread  cases.  Directionally  broad  distribu¬ 
tions  with  large  asymmetries  were  well-resolved,  but  had  a  tendency  to  bias 
distribution  peaks  toward  high-energy  distribution  tails.  Results  were  worst, 
but  still  qualitatively  good,  for  broad  distributions  with  small  asymmetries. 
Evaluations  of  test  distributions  are  necessarily  qualitative  because  there  are 
many  distribution  shapes  that  represent  the  full  suite  of  observations,  and  no 
single  quantitative  statistic  can  represent  directional  resolution  in  much  detail. 

However,  quantitative  estimates  of  error  can  be  ascribed  to  the  set  of  direc¬ 
tional  parameters  that  represent  the  test  distributions.  For  the  500  sample 
runs  of  the  110  test  cases  at  each  of  the  three  test  frequencies,  means  and 
standard  deviations  were  computed  for  peak  direction,  mean  direction,  circular 
width,  skewness,  kurtosis,  directional  spread,  and  asymmetry.  When  ordered 
by  directional  spread  and  asymmetry  classes,  these  means  and  standard  devia¬ 
tions  give  quantitative  measures  of  directional  resolving  ability  of  the  Harvest 
Platform  array.  Results  varied  widely  from  parameter  to  parameter,  and  the 
reader  is  referred  to  Table  3  and  Figures  36  to  42  in  Chapter  10  for  specific 
numerical  values.  Essentially,  these  tests  quantify  and  substantiate  observa¬ 
tions  deduced  from  the  directional  distribution  tests  discussed  in  Chapter  9. 


Closing  Comment 

It  should  be  noted  that  the  Harvest  Platform  array  is  the  only  long-term, 
deepwater,  high-resolution  directional  wave  gauge  in  existence.  Tests  dis¬ 
cussed  in  this  report  indicate  that  the  gauge  is  reliable,  and  gives  generally 
excellent  results  for  the  most  common  frequency-direction  spectra  observed  at 
this  site.  Though  the  array  makes  optimum  use  of  the  space  available  in  the 
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platfomi  structure,  it  still  has  quite  limited  spatial  coverage,  and  so  is  less  able 
to  resolve  detailed  directional  distributions  of  low-frequency  wind  waves,  i.e., 
those  with  long  wavelengths.  This  property  may  make  the  array  less  than 
ideal  for  use  in  numerical  models  having  extreme  directional  sensitivity. 
Nonetheless,  observations  made  with  this  array  are  more  than  adequate  for  a 
multitude  of  other  purposes,  among  them  general  directional  wind  wave  cli¬ 
mate,  deepwater  directional  spectral  evolution  during  storms,  intercomparison 
tests  with  and  development  of  commonly  used  low-resolution  directional  gaug¬ 
es,  and  ground  truth  in  remote  sensing  applications.  Resulting  improved 
knowledge  of  deep  ocean  wind  wave  conditions  can  only  enhance  abilities  to 
model  the  physics  of  open  ocean  wave  processes,  and  the  consequent  near¬ 
shore  wave  climate  required  in  coastal  engineering  computations  as  such 
waves  propagate  landward. 
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Figure  A3.  Bulk  data  for  March  1993  (Continued) 
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Figure  A4.  (Concluded) 


Appendix  A  Time  Series  Graphs  of  Bulk  Parameters 


A9 


100  ^  "  primary 

^ "  secondary 
A  =  tertiary 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31 

May  1993 


Figure  A5.  Bulk  data  for  May  1993  (Continued) 


Appendix  A  Time  Series  Graphs  of  Bulk  Parameters 


Appendix  A  Time  Series  Graphs  of  Bulk  Parameters 


All 


L9  (deg)  B  (deg) 

90  180  -180  0  180 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  16  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30 


*S< 


I 


++ 

VHL#-+-iJ' 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30 

June  1993 


Figure  A6.  Bulk  data  for  June  1993  (Continued) 
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Figure  A11.  (Concluded) 
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Figure  A12.  Bulk  data  for  December  1993  (Continued) 
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Figure  B2.  Spectra,  top  to  bottom:  A6  =  28°,  A  =  0.18; 
A  =  -0.07;  and  A9  =  45°,  A  =  0A3 
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Figure  B9.  Spectra,  top  to  bottom:  A6  =  105°,  >1  =  1.26; 

Ae  =  113°,  A  =  1.22;  and  A0  =  99°,  >4  =  1.50 
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Figure  C2.  Distribution  functions  with  LQ  in  the  range  15  to  25  deg 
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Figure  C7.  Distribution  functions  with  L9  in  the  range  65  to  75  deg 
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Figure  C8.  Distribution  functions  with  A0  in  the  range  75  to  85  deg 
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Figure  CIO.  Distribution  functions  with  A©  in  the  range  95  to  105  deg 
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Figure  C11.  Distribution  functions  with  LQ  in  the  range  105  to  115  deg 
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Figure  C12.  Distribution  functions  with  A0  in  the  range  115  to  125  deg 
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Figure  D1.  Test  results  for  5  <  A0  (deg)  <  15,  and  -0.2  <  A  <  0.2 
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Figure  D2.  Test  results  for  15  <  L9  (deg)  <  25,  and  -1.0  <  A  <  -0.6 
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Figure  D8.  Test  results  for  35  <  A0  (deg)  <  45,  and  -1.4  <  A  <  -1.0 
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Figure  Dll.  Test  results  for  55  <  A8  (deg)  <  65,  and  -1.8  <  A  <  -1.4 
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Figure  D12.  Test  results  for  55  <  AS  (deg)  <  65.  and  -0.2  <  A  <  0.2 
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Figure  D14.  Test  results  for  95  <  A0  (deg)  <  105,  and  -2.2  <  A  <  -1.8 
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Figure  D16.  Test  results  for  95  <  A6  (deg)  <  105,  and  1.8  <  A  •<  2.2 
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Figure  D17.  Test  results  for  115  <  A0  (deg)  <  125.  and  -2.2  <  A  <  -1.8 
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Figure  D18.  Test  results  for  115  <  A9  (deg)  <  125,  and  -0.2  <  A  <  0.2 
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Appendix  E 
Notation 


*  Complex  conjugation  operator 

Og  Normalizing  coefficient  in  MLE  to  ensure  unit  area  in 

a  directional  distribution  estimate 

Normalizing  coefficient  for  r'*  iteration  of  IMLE  to 
ensure  unit  area  in  a  directional  distribution  estimate 

A  Bulk  quartile  asymmetry  parameter 

Frequency-dependent  quartile  asymmetry  parameter  at 
n"'  discrete  frequency 

C.j(f^)  Coincident  spectral  density  between  gauges  i  and  j  at 

frequency  /„ 

d  Water  depth 

df  Discrete  frequency  increment 

dd  Discrete  direction  increment 

D(d^)  Bulk  directional  distribution  function  formed  by  inte¬ 

grating  frequency-direction  spectrum  with  respect  to 
frequency,  and  normalizing  by  total  variance 

D(fn,0m)  Directional  distribution  function  at  frequency 

Mean  directional  distribution  function  at  test  frequency 
,  derived  as  the  average  at  each  6^  of  a  number  of 
synthetic  estimates 


Appendix  E  Notation 


DMM 


MLE  estimate  of  directional  distribution  function  at 
frequency  /„ 

Mean  directional  distribution  function  displaced  by  its 
mean  direction  6^  for  the  i"‘  asymmetry  class,  and  they* 
directional  spread  class 

k*  directional  distribution  function  displaced  by  its 

mean  direction  in  the  i*  asymmetry  class  and  y* 
directional  spread  class 

Jt*  directional  distribution  function  at  test  frequency  /„ 
of  a  set  of  estimates  subject  to  synthetic  noise 

IMLE  estimate  of  directional  distribution  function  at 
frequency  after  r*  iteration 


E[] 

f 

fn 


Intermediate,  uncorrected  IMLE  estimate  of  directional 
distribution  function  at  frequency  during  r*  itera¬ 
tion 

Standard  deviation  directional  distribution  function  at 
test  frequency  ,  derived  as  the  square  root  of  the 
variance  at  each  of  a  number  of  synthetic  estimates 

Standard  deviation  directional  distribution  function 
displaced  by  its  mean  direction  6^  for  the  f*  asymme¬ 
try  class,  and  the  y*  directional  spread  class 

Unit  vector  in  the  x -direction 

Unit  vector  in  the  y -direction 

Expectation  operator 

Frequency 

71*  discrete  wave  frequency 


y^  Spectral  peak  frequency 

g  Gravitational  acceleration 


H  Characteristic  wave  height 

mo 


E2 
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i  Complex  notation  ^-1  [in  exponent  or  on  main  equa¬ 

tion  line] 

Gauge  index  [as  subscript  in  all  chapters  except  Chap¬ 
ter  8] 

Asymmetry  class  index  [as  subscript  Chapter  8] 

I  Number  of  gauges  in  an  array 

)  Integral  of  D(d^)  beginning  at  direction  6^^ 

)  Integral  at  the  n*  discrete  frequency  of 

beginning  at  direction  6^ 

I  Identity  matrix 

j  Gauge  index  [as  subscript  in  all  chapters  except  Chap¬ 

ter  8] 

Directional  spread  class  index  [as  subscript  in  Chap¬ 
ter  8] 

k  Index  of  contributing  distribution  function  [as  subscript 

in  Chapter  8] 

Directional  distribution  estimate  counter  [as  subscript 
in  Chapter  9] 

Magnitude  of  wave  number  vector  associated  with  n"' 
discrete  frequency 

"k  (6  )  Wave  number  vector  for  wave  direction  associated 

with  n*  discrete  frequency 

/  Summation  index 

m  Discrete  direction  index  [as  subscript] 

ffi.  First  cosine  moment  of  D(ff^) 

ntj  Second  cosine  moment  of  D {6^) 

m,„  First  cosine  moment  of 

Second  cosine  moment  of  D(/„,0„) 
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E3 


'"min 

Discrete  index  of  direction  of  minimum  D{f^,6^) 

M 

Upper  limit  of  direction  index  m 

M,(fn) 

Element  of  cross-spectral  matrix  between  gauges  i 
and  j  at  frequency 

Element  of  cross-spectral  matrix  with  added  synthetic 
noise  between  guages  i  and  j  at  test  frequency 

Estimate  of  element  of  cross-spectral  matrix  between 

gauges  i  and  j  at  frequency  during  r*  iteration  of 
IMLE 

M 

Cross-spectral  matrix 

M 

Cross-spectral  matrix  with  added  synthetic,  Gaussian 
noise 

n 

Discrete  frequency  index  [as  subscript] 

First  sine  moment  of  D(0^) 

Second  sine  moment  of  D{d^) 

First  sine  moment  of  D{f^,6^) 

«2„ 

Second  sine  moment  of 

N 

Upper  limit  of  discrete  frequency  index 

Qijifn) 

Quadrature  spectral  density  between  gauges  i  and  j 

r 

Iteration  count  for  IMLE 

R 

Upper  limit  for  IMLE  iterations 

(fnA) 

Sum  of  synthetic  directional  distribution  functions  at 
each  for  test  frequency 

(fn>V 

Sum  of  squares  of  synthetic  directional  distribution 
functions  at  each  0^  for  test  frequency 
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s,  ..(0  -  Sum  of  contributing  directional  distribution  fimctions 

ly  ^  m  0 ' 

shifted  by  their  mean  directions  in  quartile  asymmetry 
class  i  and  directional  spread  class  j 

s  ..(0  -0^)  Sum  of  squares  of  contributing  directional  distribution 

2y '  m  0 '  *  ,  .  - 

functions  shifted  by  their  mean  directions  m  quartile 
asymmetry  class  i  and  directional  spread  class  j 

S(f^)  Frequency  spectrum 

S(0^)  Direction  spectrum 

S(fn,  0^)  Frequency-direction  spectrum 

S(f^,  0^)  Directional  spectrum  at  spectral  peak  frequency 

T  Matrix  transpose  operator  [as  subscript  in  Chapter  9] 

U  Eigenvector  matrix  of  M 

V  Matrix  product  of  conjugate  of  eigenvector  matrix  and 
square  root  of  eigenvalue  matrix 

V  Matrix  product  of  V  and  column  vector  Z  of  random 
Gaussian  deviates 

Position  vector  of  gauge 

Position  vector  of  j"'  gauge 

jc  Horizontal  coordinate  increasing  northward 

y  Horizontal  coordinate  increasing  westward 

Z  Column  vector  of  random  Gaussian  deviates 

/3  Exponential  convergence  rate  parameter  in  IMLE 

7  Convergence  coefficient  in  IMLE 

Circular  skewness 

7^  Frequency-dependent  circular  skewness  at  frequency  /„ 

r,  (/,)  Coherence  of  signals  from  gauges  i  and  j  at  frequen¬ 

cy /n 
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8 


Circular  kurtosis 


8^  Frequency-dependent  circular  kurtosis  at  frequency 

A6  Quartile  directional  spread 

Ad^  Frequency-dependent  quartile  directional  spread  at 

frequency  /„ 

Convergence  check  parameter  at  r""  IMLE  iteration 
Mean  direction 

Frequency-dependent  mean  direction  at  frequency 

02555  First  quartile  direction  of  D (6 J 

0^^  Second  quartile  direction  of  D(0^) 

0^5^  Third  quartile  direction  of  D(0^) 

^25%  It  frequency-dependent  quartile  direction  of 

^(fit’V  at  frequency 

050%  n  Second  frequency-dependent  quartile  direction  of 

at  frequency /„ 

0^5  „  Third  frequency-dependent  quartile  direction  of 

at  frequency /„ 

0,  I"'  discrete  direction 

0  m"'  discrete  direction 

m 

0  Direction  of  minimum  energy 

^min 

0  Peak  direction  of  D(0) 

0^  Peak  direction  of  D(f^,0^)  at  frequency 

X  Wavelength 

\(f  ,0  )  IMLE  correction  function  at  the  r"’  iteration 

A  Diagonal  matrix  containing  eigenvalues  of  M 
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a  Circular  width 

Frequency-dependent  circular  width  at  frequency 

Cross-spectral  phase  function  of  signals  between  gaug¬ 
es  i  and  j  at  frequency  /„ 
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